﻿ UNIVERSITATEA “POLITEHNICA” TIMIŞOARA FACULTATEA DE ELECTRONICĂ ŞI TELECOMUNICAŢII Teza de doctorat Contribuţii la modulaţia codată trellis multidimensională folosită în realizarea modemurilor analogice de bandă vocală Coordonator ştiinţific: prof dr ing Miranda NAFORNIŢĂ Doctorand: ing Florin DĂRĂBAN -2003- -Teza de doctorat- CUPRINS CAPITOLUL I 3 INTRODUCERE 3 VOCALĂ 1 1 MODEMURI ANALOGICE DE BANDĂ 3 RIMI CARACTERISTICE 4 1 2 MĂ 1 3 CANALUL TELEFONIC 6 I MODEMURI DIGITALE UTILIZATE ÎN PSTN 6 1 4 MODEMURI ANALOGICE Ş 1 5 SCURT ISTORIC 8 CAPITOLUL II 10 SISTEME DE TRANSMISII DE DATE CU MODULAŢIE CODATĂ TRELLIS (TCM) 10 2 1 INTRODUCERE 10 2 1 1 Detecţia simbol cu simbol (SSD) 11 2 1 2 Detecţia de probabilitate condiţionată maximă (MLD) 13 2 2 DIAGRAMA TRELLIS 14 IONAREA CONSTELAŢIEI DE SEMNALE 14 2 3 PARTIŢ TORUL TCM 16 2 4 TRANSMIŢĂ 2 4 1 Reprezentarea Ungerboeck a transmiţătorului TCM 16 2 4 2 Reprezentarea Calderbank-Mazo a transmiţătorului TCM 18 2 5 RECEPTORUL TCM 23 2 5 1 Receptorul TCM cu detecţie de probabilitate condiţionată maximă (MLD) 23 2 5 2 Receptorul TCM cu detecţie simbol cu simbol (SSD) 25 2 5 3 Algoritmul Viterbi 26 CAPITOLUL III 29 PERFORMANŢELE SISTEMELOR DE TRANSMISII DE DATE CU MODULAŢIE CODATĂ TRELLIS 29 A PROBABILITĂŢII ERORII DE SECVENŢĂ 3 1 LIMITA SUPERIOARĂ 29 3 1 1 Diagrama de stare a erorii 31 3 1 2 Consideraţii de simetrie a funcţiei de transfer matriciale a diagramei de stare a erorii 32 3 1 3 Consideraţii asimptotice a funcţiei de transfer matriciale a diagramei de stare a erorii 34 3 1 4 Limita superioară a probabilităţii erorii de bit 35 3 1 5 Consideraţii de convergenţă a funcţiei de transfer scalare a diagramei de stare a erorii 36 3 1 6 Cazul unui canal de comunicaţie general 36 A PROBABILITĂŢII ERORII DE SECVENŢĂ 3 2 LIMITA INFERIOARĂ 38 3 2 1 Limita inferioară a probabilităţii erorii de bit 39 IEI DE TRANSFER SCALARĂ A DIAGRAMEI DE STARE A ERORII 39 3 3 CALCULUL FUNCŢ EI EUCLIDIENE MINIME 40 3 4 CALCULUL DISTANŢ 3 4 1 Calculul distanţei euclidiene minime utilizînd diagrama de stare a erorii 40 3 4 2 Algoritmul Saxona-Mulligan-Wilson 41 3 4 3 Algoritmul produs trellis 44 3 4 4 Limita inferioară a distanţei euclidiene minime 45 3 4 5 Limita superioară a distanţei euclidiene minime 46 DE PUTERE A SEMNALULUI DE LINIE TRANSMIS 51 3 5 DENSITATEA SPECTRALĂ 1 -Teza de doctorat- CAPITOLUL IV 54 CODURI TCM MULTIDIMENSIONALE 54 4 1 INTRODUCERE 54 IA WEI PENTRU CODURI TCM MULTIDIMENSIONALE 55 4 2 CONSTRUCŢ 4 2 1 Codul TCM 4-D, cu rata 2/3 şi cu 16 stări 56 4 2 2 Codul TCM 8-D, cu rata 3/4 şi cu 64 stări 64 IA STERIAN PENTRU CODURI TCM MULTIDIMENSIONALE 71 4 3 CONSTRUCŢ 4 3 1 Codul TCM 6-D, cu rata 3/4 şi cu 64 stări 74 4 3 2 Codul TCM 12-D, cu rata 4/5 şi cu 256 stări 82 IEI WEI PENTRU CODURI TCM MULTIDIMENSIONALE 92 4 4 GENERALIZAREA CONSTRUCŢ 4 4 1 Construcţia constelaţiei de semnale 2N-D 92 4 4 1 1 Metode de construcţie a unei constelaţii de semnale 2-D optimală 92 4 4 1 1 1 Modelarea constelaţiei 92 4 4 1 1 2 Corespondenţă în inele 97 4 4 1 2 Extinderea constelaţiei de semnale 2-D constituentă a constelaţiei de semnale 2N-D 100 4 4 1 2 1 Extinderea simplă a constelaţiei de semnale 2-D 101 4 4 1 2 2 Extinderea optimală a constelaţiei de semnale 2-D 104 4 4 2 Partiţionarea constelaţiei de semnale 2N-D în subseturi 2N-D 112 4 4 3 Asignarea subseturilor 2N-D la tranziţiile de stare ale codului TCM 2N-D 113 4 4 4 Realizarea corespondenţei dintre cei NQ+1 biţi şi constelaţia de semnale 2N-D 114 IEI WEI A CODORULUI TCM 2N-D 114 4 5 GENERALIZAREA CONSTRUCŢ IEI WEI A DECODORULUI TCM 2N-D 120 4 6 GENERALIZAREA CONSTRUCŢ IA WEI GENERALIZATĂ A CODURILOR TCM 2N-D 127 4 7 EXEMPLE PENTRU CONSTRUCŢ 4 7 1 Codul TCM 4-D, cu rata 2/3 şi cu 16 stări 127 4 7 2 Codul TCM 6-D, cu rata 3/4 şi cu 64 stări 134 4 7 3 Codul TCM 8-D, cu rata 3/4 şi cu 64 stări 143 4 7 4 Codul TCM 12-D, cu rata 4/5 şi cu 256 stări 150 CAPITOLUL V 160 MODEMURI ANALOGICE DE BANDĂ VOCALĂ CU MODULAŢIE CODATĂ TRELLIS MULTIDIMENSIONALĂ 160 IONALĂ 5 1 SCHEMA BLOC FUNCŢ 160 5 1 1 Cifratorul şi decifratorul 160 5 1 2 Convertorul serie/paralel şi convertorul paralel/serie 162 5 1 3 Codorul neliniar 163 5 1 4 Filtrul de preaccentuare 164 5 1 5 Egalizarea adaptivă 167 5 1 6 Modulatorul QAM Demodulatorul QAM Blocul de refacere a semnalului purtător 171 5 1 7 Filtrul de interpolare/decimare Blocul de refacere a tactului de simbol 175 5 1 8 Blocul de control automat a cîştigului 178 5 1 9 Supresorul de ecou 181 5 2 IMPLEMENTARE HARD 183 5 2 1 Procesorul digital de semnal (DSP) 184 5 2 2 Memoria SRAM 184 5 2 3 Memoria ROM 184 5 2 4 Interfaţa cu calculatorul gazdă 185 5 2 5 Interfaţa analogică 185 5 2 6 Interfaţa cu linia telefonică 186 5 3 IMPLEMENTARE SOFT 188 CAPITOLUL VI 190 CONCLUZII ŞI CONTRIBUŢII PERSONALE LA MODULAŢIA CODATĂ TRELLIS MULTIDIMENSIONALĂ 190 IA CODATĂ TRELLIS MULTIDIMENSIONALĂ 6 1 CONCLUZII PERSONALE CU PRIVIRE LA MODULAŢ 190 II PERSONALE LA MODULAŢIA CODATĂ TRELLIS MULTIDIMENSIONALĂ 6 2 CONTRIBUŢ 193 BIBLIOGRAFIE CITATĂ ÎN TEZA DE DOCTORAT 195 BIBLIOGRAFIE CONSULTATĂ PENTRU DOCUMENTARE 199 ANEXA 1 Tabel cu parametrii calculaţi ai codurilor TCM 2N-D pentru N=2≤8 204 ANEXA 2 Modelele MATLAB-Simulink pentru sistemele de transmisii de date cu TCM 4-D (construcţia Wei), TCM 6-D (construcţia Sterian) şi TCM 6-D (construcţia propusă) -205 ANEXA 3 Rezultatele simulărilor - Graficele ratei erorilor de bit în funcţie de raportul semnal/zgomot al canalului de comunicaţie AWGN -206 2 -Teza de doctorat- CAPITOLUL I Introducere 1 1 Modemuri analogice de bandă vocală În domeniul dinamic al telecomunicaţiilor apare adesea tendinţa de a prezice declinul rapid al unei tehnologii existente în momentul în care o tehnologie alternativă mai puternică apare la orizont Acesta este şi cazul modemurilor analogice de bandă vocală, despre care în anii 1985 se credea că vor dispare complet odată cu apariţia reţelei digitale cu integrarea serviciilor (ISDN=Integrated Services Digital Network) Cu toate acestea, în prezent, în condiţiile dezvoltării reţelei Internet şi a utilizării calculatoarelor personale (PC=Personal Computer), modemurile analogice sunt cele mai răspîndite vehicule pentru comunicaţiile de date În anul 1998, conform unui studiu realizat de Georgia Tech, în SUA, 69,6 % din toţi utilizatorii reţelei Internet (persoane fizice şi persoane juridice) au accesat reţeaua Internet prin modemuri analogice În acest studiu se prezintă repartizarea vitezei de conectare în funcţie de locaţie, sex, grupe de vârstă, grupe de experienţă În Tabelul 1 1 se prezintă repartizarea vitezei de conectare în funcţie de locaţia geografică a utilizatorilor reţelei Internet Tabelul 1 1 Repartizarea vitezei de conectare în funcţie de locaţia geografică a utilizatorilor persoane fizice şi persoane juridice a reţelei Internet Tip acces Viteza de conectare SUA Europa Altele (răspunsul subiecţilor [%] [%] [%] chestionaţi) Modemuri analogice 45 Mb/s (FDDI**) 0,6 4 0,8 Total 26,4 53,8 34 Nu ştiu ce modem Nu ştiu ce modem 3,9 6,7 3,3 folosesc folosesc *T1,T3=Flux PCM (Pulse Code Modulation) primar, ternar în SUA/Japonia **FDDI=Fiber Distributed Data Interface 3 -Teza de doctorat- În anul 2000, conform unui studiu realizat de Forrester Research, în SUA, 87,9 % din utilizatorii reţelei Internet (persoane fizice) au accesat reţeaua Internet prin modemuri analogice, iar 80 % din aceştia cu un modem analogic de 56 kb/s (conform recomandării ITU V 90) (Tabelul 1 2) Tabelul 1 2 Repartizarea pe modemuri a utilizatorilor reţelei Internet (persoane fizice) Tip acces 1998 1999 2000 [%] [%] [%] Modemuri analogice 96,5 92,1 87,9 Modemuri digitale Cablu coaxial 2,4 5,7 9,5 ISDN, DSL*, Wireless** 1,1 2,3 2,6 *=Digital Subscriber Line **=Fără fir La succesul modemurilor analogice au contribuit : • progresele în tehnicile de codare, modulaţie şi egalizare adaptivă, • progresele în tehnologiile de prelucrare digitală a semnalelor, atît în ceea ce priveşte algoritmii de prelucrare digitală a semnalelor, cît si a procesoarelor digitale de semnal (DSP=Digital Signal Processor), care sunt din ce în ce mai performante şi mai ieftine, • dezvoltarea reţelei telefonice publice comutate (PSTN=Public Switched Telephone Network) şi creşterea calităţii conexiunilor telefonice, prin introducerea pe scară largă a centralelor telefonice digitale şi a echipamentelor de transmisiuni digitale, • standardizarea realizată de Uniunea Internaţională a Telecomunicaţiilor- sectorul de standardizare în Telecomunicaţii (ITU-T=International Telecommunications Union-Telecommunication Standardization Sector) (fostă CCITT = Consultative Committee Intenational Telegraph and Telephone) care începînd cu anii 1970 a adoptat succesiv standardele din Seria V pentru modemurile analogice 1 2 Mărimi caracteristice În cazul unei transmisii sincrone, definim debitul binar (rata de bit) Db [bit/s] ca numărul de biţi ce pot fi transmişi în unitatea de timp prin canalul de comunicaţie : 1 [bit/s] (1 1) bD = bT unde Tb = durata unui bit [s] În cazul unei transmisii asincrone, definim viteza de modulaţie (rata de simbol) vm [simboluri/s] sau [Baud] ca viteza cu care se schimbă stările semnalului modulat în canalul de comunicaţie : 4 -Teza de doctorat- 1 [simboluri/s] sau [Baud] (1 2) mv= sT unde Ts = durata unui simbol [s] m Dacă semnalul modulat are M=2 nivele atunci avem relaţ ia: vMlogvD= m [bit/s] (1 3) m2mb= unde m = numărul de biţi necesari codării În cazul unei transmisii oarecare, definim debitul binar ca raportul dintre numărul total de biţi transmişi (biţi utili + biţi suplimentari ceruţi de protocol) şi durata necesară transmisiei lor Eficienţa spectrală η [bit/s/Hz] este definită ca numărul de biţi per secundă pe care îi putem transmite într-o lăţime de bandă de 1 Hz a canalului de comunicaţie : D b = [bit/s/Hz] (1 4) η B unde B = lăţimea de bandă a canalului de comunicaţie care este disponibilă pentru transmisie [Hz] iar Db = debitul binar [bit/s] Rata de transmisie Q [bit/utilizare canal] sau [bit/interval de semnalizare] este definită ca numărul de biţi per utilizarea canalului sau ca numărul de biţi per interval de semnalizare : [bit/utilizare canal] sau [bit/interval de semnalizare] (1 5) TQ= bsD unde prin utilizarea canalului de comunicaţie se înţelege transmisia unui simbol de informaţie prin canalul de comunicaţie Viteza maximă de transmisie a datelor printr-un canal de comunicaţie ideal (stabil, omogen, invariant în timp, fără distorsiuni) este limitată Această limită, numită capacitatea canalului de comunicaţie C [bit/s], a fost stabilită de Shannon pentru un canal de comunicaţie ideal cu zgomot aditiv, alb, gaussian, cu o putere medie a semnalului util mărginită şi presupunînd o codare infinită ca fiind: S BC+= )(1log [bit/s] (1 6) 2 N unde C = capacitatea canalului de comunicaţie [bit/s], S = puterea medie a semnalului util [W] iar N = puterea medie a zgomotului [W] S Dacă raportul semnal/zgomot se consideră în dB deci 10SNR= )(log 10 N atunci: [bit/s] (1 7) 0,33C≅ B SNR 5 -Teza de doctorat- 1 3 Canalul telefonic Canalul telefonic poate fi considerat ca un canal de comunicaţie cu zgomot redus, de bandă de frecvenţe limitată şi de putere limitată Canalul telefonic poate fi presupus cu o bună aproximaţie ca un canal de comunicaţie liniar Canalul telefonic este un canal de comunicaţie invariant în timp şi are o caracteristică de frecvenţă relativ uşor egalizabilă Un canal telefonic obişnuit are banda de frecvenţe cuprinsă între 300 şi 3400 Hz deci B=3100 Hz În toate cazurile de conexiune telefonică avem SNR=28 dB, iar în majoritatea cazurilor de conexiune telefonică avem SNR=32-38 dB În general, canalul telefonic este considerat ca un canal de comunicaţie cu zgomot aditiv, alb, gaussian (AWGN=Aditive White Gaussian Noise) şi cu interferenţă intersimbol (ISI= InterSymbol Interference) 1 4 Modemuri analogice şi modemuri digitale utilizate în PSTN În Tabelul 1 3 se prezintă caracteristicile modemurilor analogice şi digitale utilizate în PSTN Modemurile sunt conectate la PSTN prin perechi de fire de cupru, torsadate, neecranate, cu diametrul cuprins între φ =0,3-0,4 mm (cabluri urbane) şφ=0,8-0,9 mm i (cabluri interurbane) Se observă că modemurile analogice oferă debite binare mici în comparaţie cu modemurile digitale Din acest motiv modemurile analogice sunt aproape imposibil de utilizat în aplicaţii cum ar fi: ¾ acces de foarte mare viteză la reţeaua Internet, ¾ acces la servicii interactive de bandă largă, ¾ videoconferinţe Modemurile digitale permit debite binare mari deoarece au fost proiectate pentru o banda de frecvenţe a canalului de comunicaţie de ordinul a cîţiva MHz deci mult mai mare decît B=3000-3500 Hz ca în cazul modemurilor analogice În cazul modemurilor digitale nu se mai folosesc ca şi la modemurile analogice filtrele trece bandă de 4 kHz (antialiasing şi interpolare) în CODEC-urile interfeţelor de linie din unităţile de racordare locale/distante ale centralei telefonice sau din unităţile de reţea optică (ONU=Optical Network Unit) ale reţelei de acces (AN=Acces Network) 6 -Teza de doctorat- Tabelul 1 3 Caracteristicile modemurilor analogice şi digitale utilizate în PSTN Tip modem Db [kbit/s] ηModulaţia Linii de transmisie Lungimea [bit/s/Hz] (cod de linie) maximă a liniei de transmisie [Km] (φ=0,3-0,4 mm) MODEMURI ANALOGICE V 21 0,3 1 FSK comutate/închiriate 1,8 V 22 0,6 2 PSK comutate/închiriate 1,8 V 22 bis 2,4 4 QPSK comutate 1,8 V 32 9,6 4 TCM comutate/închiriate 1,8 V 32 bis 14,4 6 TCM comutate/închiriate 1,8 V 33 14,4 6 TCM închiriate 1,8 V 34 28,8 8,4 TCM comutate/închiriate 1,8 V 34 bis 33,6 9,8 TCM comutate/închiriate 1,8 MODEMURI DIGITALE ISDN BRI 144 2 PAM(4B3T) comutate/închiriate 4 ISDN PRI 2048 2 PCM(HDB3) comutate/închiriate 0,3 HDSL 2048 2 PAM(2B1Q) închiriate 4 sau CAP SDSL 768 2 PAM(2B1Q) închiriate 3,7 IDSL 128 2 PAM(2B/1Q) închiriate 4 ADSL Downstream= 8 CAP sau închiriate 3-5,5 =8448-1554 DMT Upstream= =640-16 VDSL Downstream= 4 CAP sau închiriate 0,3-1 =52000-13000 DMT Upstream= =2300-1500 ISDN BRI = Integrated Services Digital Network Base Rate Interface ISDN PRI = Integrated Services Digital Network Primary Rate Interface HDSL = High bit rate Digital Subscriber Line SDSL = Symmetric Digital Subscriber Line IDSL = ISDN Digital Subscriber Line ADSL = Asymmetric Digital Subscriber Line VDSL = Very high bit rate Digital Subscriber Line FSK = Frequency Shift Keying PSK = Phase Shift Keying QPSK = Quadrature Phase Shift Keying TCM = Trellis Coded Modulation PAM = Pulse Amplitude Modulation 2B1Q = 2 Binary 1 Quaternary 4B3T = 4 Binary 3 Ternary PCM = Pulse Code Modulation HDB3 = High Density Bipolar 3 CAP = Carrierless Amplitude and Phase Modulation DMT = Discrete MultiTone modulation 7 -Teza de doctorat- În plus modemurile digitale folosesc tehnici digitale de modulaţie ca: • modulaţia impulsurilor în amplitudine (PAM=Pulse Amplitude Modulation) cu un cod de linie cuaternar 2B1Q (2 Binary 1 Quaternary) • modulaţia impulsurilor în cod (PCM=Pulse Code Modulation) cu un cod de linie ternar HDB3 (High Density Bipolar 3) • modulaţia de amplitudine şi de fază fără purtătoare (CAP=Carrierless Amplitude and Phase Modulation) în care se realizează 2 modulaţii digitale cu ajutorul a două filtre digitale transversale, care au funcţiile de transfer cu aceeaşi amplitudine dar defazate cu 90° Semnalele de la ieşirile filtrelor digitale sunt însumate iar semnalul obţinut este convertit cu un convertor digital-analogic (DAC=Digital-to-Analog Converter) pentru a fi transmis prin canalul de comunicaţie • modulaţia multiton discretă (DMT=Discrete MultiTone modulation) în care se realizează o multiplexare prin diviziune în frecvenţă (FDM=Frecquency Division Multiplexing), a mai multor canale de 4 kHz, cu modulaţie codată trellis (TCM=Trellis Coded Modulation) Pentru fluxul de date spre utilizator (Downstream) sunt multiplexate 249 de canale, în banda de frecvenţe (133,8- 1100) kHz iar pentru fluxul de date dinspre utilizator (Upstream) sunt multiplexate 25 de canale, în banda de frecvenţe (26-133,8) kHz Majoritatea canalelor sunt folosite pentru transmisii de date ( 15η0≤≤ bit/s/Hz) dar există şi canale folosite pentru sincronizare Semnalele purtătoare sunt distanţate la intervale de 4,3125 kHz Cu toate acestea modemurile analogice fiind mult mai ieftine şi fiind primele care au apărut pe piaţă sunt mult mai răspîndite decît modemurile digitale 1 5 Scurt istoric În 1974 Massey a sugerat formal că performanţele sistemelor digitale de telecomunicaţii ar putea fi îmbunătăţite prin folosirea codării şi modulaţiei ca o entitate combinată şi nu ca două operaţii separate În 1976 Ungerboeck şi Csajka prezintă ideea de modulaţie codată trellis (TCM=Trellis Coded Modulation) la simpozionul de teoria informaţiei din Suedia În 1977 Imai şi Hirakawa au publicat o lucrare în care codarea bloc sau codarea convoluţională au fost combinate cu modulaţia În 1982 Ungerboeck a descris principiile de baza ale TCM În 1984 Forney a arătat că utilizarea TCM la construcţia modemurilor analogice permite obţinerea unor debite binare mai apropiate de capacitatea canalului de comunicaţii (canalul telefonic) în comparaţie cu alte medii de transmisii fizice În 1984 Wei a prezentat o schemă TCM 2-D cu 8 stări care a fost adoptată în acelaşi an de către ITU-T ca recomandarea V 32 pentru modemurile analogice care operează full-duplex la 9600 bit/s pe linii comutate sau închiriate cu 2 fire În 1987 Wei a prezentat lucrarea „TCM cu constelaţii de semnale multidimensionale” în care a prezentat trei scheme TCM pentru Q=7 biţi/interval de semnalizare şi anume: 8 -Teza de doctorat- ¾ o schemă TCM 4-D cu 16 stări, ¾ o schemă TCM 4-D cu 64 stări, ¾ o schemă TCM 8-D cu 64 stări În 1994 ITU-T a adoptat recomandarea V 34 bis pentru modemurile analogice care operează full-duplex la 33600 bit/s pe linii comutate sau închiriate cu 2 fire Se folosesc trei scheme TCM si anume: ¾ o schemă TCM 4-D cu 16 stări propusă de Wei ¾ o schemă TCM 4-D cu 32 stări propusă de Williams ¾ o schemă TCM 4-D cu 64 stări propusă şi corectată de Wei după ce Rossin şi Heegard au descoperit şi au semnalat în 1993 o greşeală în schema iniţială În 1998 ITU-T a adoptat recomandarea V 90 pentru comunicaţia de date dintre un modem analogic şi un modem digital care operează full-duplex la 56000 bit/s în sensul modem digital → modem analogic şi la 33600 bit/s în sensul modem analogic → modem digital pe linii comutate cu 2 fire Modemul analogic se conectează la PSTN printr-o interfaţă analogică obişnuită Modemul analogic emite semnale V 34 şi recepţionează semnale G 711 după trecerea lor printr-un decodor G 711 Modemul digital se conectează la PSTN printr-o interfaţă ISDN de bază (ISDN BRI=Integrated Services Digital Network Basic Rate Interface) sau printr-o interfaţă ISDN primară (ISDN PRI= Integrated Services Digital Network Primary Rate Interface) Modemul digital emite semnale G 711 şi recepţionează semnale V 34 după trecerea lor printr-un codor G 711 În prezent TCM rămîne un domeniu de cercetare activ în care continuă să apară idei noi ce deschid calea spre noi scheme practice Progresul este relativ lent pentru că problema găsirii unor metode noi de codare şi modulaţie este în general foarte complicată 9 -Teza de doctorat- CAPITOLUL II Sisteme de transmisii de date cu modulaţie codată trellis (TCM) 2 1 Introducere Sistemele de transmisii de date cu modulaţie codată trellis (TCM) prezintă o foarte bună protecţie la zgomot şi se folosesc la transmisii de date de foarte mare viteză În cazul unui canal de comunicaţie de putere limitată, creşterea eficienţei energetice se poate realiza prin utilizarea unor coduri corectoare de erori Codurile corectoare de erori adaugă biţi suplimentari secvenţei codificate transmise, deci impune modulatorului să lucreze la un debit binar mai mare necesitînd o bandă de frecvenţe mai mare În cazul unui canal de comunicaţie de bandă de frecvenţe limitată, creşterea eficienţei spectrale se poate realiza prin utilizarea unor modulatoare de ordin mare Modulatoarele de ordin mare necesită o putere a semnalului de linie mai mare pentru a menţine aceeaşi distanţă între semnale, cu alte cuvinte aceeaşi probabilitate de eroare Modulaţia codată trellis (TCM) este o tehnică ce combină utilizarea unor coduri corectoare de erori (în general coduri convoluţionale) cu utilizarea unor modulatoare de ordin mare Ideea modulării multidimensionale a simbolurilor codate convoluţional a apărut cu mult înaintea apariţiei modulaţiei codate trellis, operaţiile de codare–decodare şi modulare-demodulare fiind independente Performanţele acestor sisteme nu sunt satisfăcătoare deoarece: ♦ în receptor deciziile se iau simbol cu simbol înainte de decodare ceea ce poate conduce la pierderi de informaţie, ♦ codurile convoluţionale optimizate după criteriul distanţei Hamming nu asigură o separare maximă între semnale O mai bună protecţie la zgomot se asigură dacă semnalele de linie emise diferă cît mai mult unul de altul şi o măsură a distanţei dintre ele este distanţa euclidiană Pentru a mări distanţa euclidiană este necesară creşterea numărului de semnale, astfel încît să apară o redundanţă la codare, al cărei scop este maximizarea distanţei euclidiene minime Aspectul nou al modulaţiei codate trellis constă în faptul că operaţiile de codare şi modulaţie nu sunt tratate ca operaţii independente ci ca o operaţie unică Semnalul de linie recepţionat, în loc să fie mai întîi demodulat şi apoi decodat, este procesat de un receptor care combină demodularea şi decodarea într-o singură operaţie Procesul de detecţie implică mai mult decizii soft decît decizii hard (semnalul de linie recepţionat este procesat înainte de a lua decizia cărui simbol al sursei de date îi corespunde) Considerăm un canal de comunicaţie cu zgomot aditiv alb gaussian Semnalele de linie transmise pot fi considerate vectori într-un spaţiu euclidian N-dimensional N (R) numit spaţ iu de semnale 10 -Teza de doctorat- Pentru fiecare simbol al sursei de date, vectorul semnal transmis s este ales dintr-un set Ω′, format din M′ vectori semnal, numit constelaţie de semnale Vectorul semnal recepţionat este: r=s+n (2 1) unde n este vectorul zgomot a cărui componente sunt variabile aleatoare gaussiene N 0 independente, cu media 0 şi dispersia 2 Energia medie a vectorilor semnal din setul Ω′ este egală cu: 21 E′=∑ s (2 2) M′ ′ ∈s Ω 1K− Fie o secvenţă de K vectori semnal transmişi {si} Receptorul TCM care 0i= minimizează probabilitatea erorii de secvenţă operează astfel: 1K− ♦ observă secvenţa de K vectori semnal recepţionaţi {ri}, 0i= 1K− pătratul ♦ decide că s-a transmis secvenţa de vectori semnal {Si}, dacă 0i= −1K 22 distanţei euclidiene d=− ∑ iisr este minim pentru si=Si, i=0,…,K-1 Cu =0i 1K− de secvenţa alte cuvinte secvenţa de vectori semnal {Si} este mai apropiată 0i= 1K− secvenţă de vectori de vectori semnal recepţionaţi {ri} decît oricare altă 0i= semnal Probabilitatea erorii de secvenţă (ca şi probabilitatea erorii de simbol) este 2 ′ mind i este bine limitată superior de o funcţie descrescătoare a raportului ş 0N aproximată de această expresie cînd raportul semnal/zgomot este mare În expresia de 2 mai sus d′ tratul distanţei euclidiene minime dintre două secvenţe posibile min este pă de vectori semnal din Ω′ 2 1 1 Detecţia simbol cu simbol (SSD) La detecţia simbol cu simbol (SSD=Symbol by Symbol Detection) secvenţa 1K− de vectori semnal {si} este o secvenţă de vectori independenţi Astfel secvenţele 0i= K posibile de vectori semnal aparţin lui Ω′ Pătratul distanţei euclidiene este minimizat 2 dacă se minimizează separat termenii sr − ii pentru si∈Ω′ Probabilitatea erorii de secvenţă (ca şi probabilitatea erorii de simbol) este limitată superior de expresia: 11 -Teza de doctorat- 2 ′ d11M′ min− ))erfc(( (2 3) P(e)≤ N22 0 fiind bine aproximată de această expresie cînd raportul semnal/zgomot este mare În relaţia (2 3): 22 ′′′ ′ (2 4) mind=ssminii− ′′′ Ωs,s′ ii∈ ∞+ 22 t− erfc(x)=dte (2 5) ∫ π x Concluzia 2 1-Problema proiectării unui sistem de transmisii de date cu modulaţie codată trellis cînd se dau N, E′, M′ se reduce la alegerea unui set Ω′ de vectori semnal astfel încît distanţa euclidiană minimă dintre oricare doi vectori semnal să fie maximă Pentru compararea a două constelaţii de semnale se definesc următoarele două mărimi: Mlog′ 2 ⇒ eficienţa spectrală: η= [bit/s/Hz] (2 6) N 2 2mind′ ⇒ eficienţa energetică: δ = [bit/s/W] (2 7) bE E′ unde Eb= = energia medie per bit de informaţ ie [W/bit] Mlog′ 2 Probabilitatea erorii de secvenţă (ca şi probabilitatea erorii de simbol) este limitată superior de expresia: Eδ 1M− b′ ))erfc(( (2 8) P(e)≤ N22 0 fiind bine aproximată de această expresie cînd raportul semnal/zgomot este mare Concluzia 2 2-Aceeaşi probabilitate a erorii de secvenţă (a erorii de simbol) poate fi bE2 obţinută cu un raport semnal/zgomot mai mic dacă este eficienţa energetică δ 0N mai mare 12 -Teza de doctorat- 2 1 2 Detecţia de probabilitate condiţionată maximă (MLD) La detecţia de probabilitate condiţionată maximă (MLD=Maximum 1K− Likelyhood Detection) secvenţa de vectori semnal {si} este o secvenţă de vectori 0i= interdependenţi Astfel secvenţele posibile de vectori semnal aparţin unui subset din K ceea ce implică scăderea eficienţei spectrale η Ω′ Pentru a evita scăderea eficienţei spectrale se înlocuieşte Ω′cu Ω ⊃ Ω′ şi M′ cu M > M′, deci se extinde constelaţia de semnale Astfel rezultă o secvenţă de vectori semnal mai puţin asemănători între ei şi se obţine o creştere a pătratului distanţei euclidiene minime dintre două secvenţe posibile de vectori semnal Se definesc două tipuri de cîştig: 2 mind ⇒ cîştigul de distanţă: ε= (2 9) 2 ′ mind 2 d min E ⇒ cîştigul asimptotic de codare: γ= (2 10) 2 d′ min E′ 2 ei euclidiene minime dintre două secvenţe posibile de unde dmi=pă ntratul distanţ vectori semnal din setul Ω, iar E=energia medie a vectorilor semnal din setul Ω Vectorul semnal si transmis la momentul de timp discret i depinde nu numai de simbolul sursei de date ai transmis la acelaş i moment de timp discret i, ci şi de un număr finit L de simboluri anterioare ale sursei de date: si=f(ai, ai-1, ai-2,…,ai-L) (2 11) Se defineşte starea transmiţătorului TCM la momentul de timp discret i astfel : σ i=(ai-1, ai-2,…,ai-L) (2 12) Astfel se obţin relaţiile : si=f(ai, σ i) (2 13) σ i+1=g(ai, σi) (2 14) unde se presupune că funcţiile f şi g sunt invariante în timp Funcţia f arată că fiecare vector semnal depinde nu numai de simbolul corespondent al sursei de date ci şi de starea transmiţătorului TCM La fiecare moment de timp discret i, vectorul semnal si este selectat dintr-o subconstelaţ ie de semnale 13 -Teza de doctorat- selectată de valoarea stării transmiţătorului σ i Funcţia g arată evoluţia stărilor transmiţătorului TCM şi descrie memoria transmiţătorului TCM Modelul general al transmiţătorului TCM este prezentat în Fig 2 1 Fig 2 1 Modelul general al transmiţătorului TCM 2 2 Diagrama trellis Pentru reprezentarea grafică a funcţiilor f şi g se foloseşte o diagramă trellis Nodurile diagramei trellis sunt reprezentate prin valorile pe care le poate lua σ i (stările transmiţătorului TCM) Fiecărui simbol ai al sursei de date îi asociem o ramură care porneşte de la un nod al diagramei trellis la un moment de timp discret i şi care ajunge la un nod al diagramei trellis la momentul de timp discret i+1 Funcţia f determină care vector semnal si este asociat cu fiecare ramură de-a lungul diagramei trellis Funcţia g determină nodurile diagramei trellis m′ biţi/simbol sursă) Dacă avem o sursă de date care are M′=2 simboluri (m′ atunci din fiecare nod al diagramei trellis pornesc M′ ramuri (1 ramură/simbol sursă) Dacă două sau mai multe ramuri conectează aceleaşi perechi de noduri atunci apar tranziţii paralele Dacă două sau mai multe ramuri pornesc din acelaşi nod şi ajung în acelaşi nod atunci avem tranziţii adiacente 2 3 Partiţionarea constelaţiei de semnale În cazul detecţiei simbol cu simbol (SSD) diagrama trellis are la un moment de timp discret un singur nod şi toate tranziţiile sunt paralele În cazul detecţiei de probabilitate condiţionată maximă (MLD) odată observată secvenţa de vectori semnal recepţionaţi se caută traiectoria cea mai probabilă prin diagrama trellis Datorită zgomotului aditiv din canalul de comunicaţie traiectoria aleasă poate să difere de traiectoria corectă Dacă traiectoria aleasă şi traiectoria corectă diverg la 14 -Teza de doctorat- momentul de timp discret i şi converg la momentul de timp discret i+L atunci a avut loc o eroare de lungime L Atunci dmin este egală cu distanţa euclidiană minimă dinte vectorii semnal asociaţi cu o pereche de traiectorii prin diagrama trellis care formează o eroare de lungime L (Fig 2 2) 22222 )D,C(d )B,A(dd++= )B,A(dd= minmin unde X,Y sunt subseturi de vectori semnal asociaţi cu ramurile diagramei trellis şi 2 tratul distanţei euclidiene minime dintre vectorii semnal care aparţin lui )Y,X(d=pă X şi vectorii semnal care aparţin lui Y Fig 2 2 Pereche de traiectorii care formează o eroare de lungime L Concluzia 2 3-Subsetul de vectori semnal asociaţi cu tranziţii paralele trebuie să aibă o distanţă euclidiană minimă cît mai mare Concluzia 2 4-Subsetul de vectori semnal asociaţi cu tranziţii adiacente trebuie să aibă o distanţă euclidiană minimă cît mai mare Pentru partiţionarea optimă a unei constelaţii de semnale de mărime M se utilizează tehnica lui Ungerboeck (Fig 2 3) : ♦ se partiţionează succesiv constelaţia de semnale în 2, 4, 8,… subconstelaţii de MMM semnale de mărimi , , , ş i cu distanţele euclidiene minime 842 (3)(2)(1) ddd este necesar ca funcţia f să fie neliniară De fapt o funcţie f liniară corespunde efectului interferenţei intersimbol (ISI) şi determină 1γ= De exemplu, dacă simbolurile sursei de date sunt ai∈ {-1, +1} şi funcţia f este liniară de forma: f(a1, a2,…,an)=h1a1+h2a2+…+hnan (2 19) unde h1, h2,…,hn=constante reale, atunci, alegînd o traiectorie prin diagrama trellis corespunzătoare secvenţei de simboluri {+1} şi o altă traiectorie corespunzătoare aceleaşi secvenţe de simboluri dar cu un singur simbol {+1} schimbat în {-1}, obţinem: 22  )h hh()h h(hd ++++−−+++= n21n21min  2  )h h(h)h h(h ++++−−+++ n21n21  2  )h h(h)h h(h =−++−+++ n21n21  (2 20)  222 )h h4(h +++= n21   222 )h h(hE +++= n21   2 d min  4 =  E 18 -Teza de doctorat- În absenţa codării avem h1=1 şi hi=0, 1i>:  222 ′  4h)h()(hd =−−= 111min   2 ′ hE (2 21) =  1  2 ′ d  min 4 =  ′  E deci cîştigul asimptotic de codare este 1γ= Fie a=(a1, a2,…,an) un vector n-dimensional a că rui componente sunt variabile aleatoare independente, identic distribuite, care iau K valori reale dintr-un set A În general funcţia f(a) este neliniară şi poate fi scrisă ca o serie ortogonală Volterra: nnn )2()2()1()1()0( )(Qk)(Qkk)(f+ += aaa (2 22) ∑∑∑ ijijii+ == 1j1i1i= )l( iile de unde )(Qa=polinom de grad l în variabila a aa care satisface condiţ li 2i1ili2i1i (m)(l) ortogonalitate 0(a)](a)QE[Q= dacă l≠m sau dacă l=m dar j1j2…jm nu jjj iii m21l21 n şi m≤n este o permutare a lui i1i2…il unde l≤ )l( inem o serie ortonormală Dacă normalizăm polinoamele )(Qa obţ li 2i1i Volterra n Teorema 2 1-Seria ortonormală Volterra f(a) are K termeni l indicii obţinuţi din indicii i1i2…il după Demonstraţie-Se notează j1j2…js cu s≤ eliminarea repetiţiilor Se notează vk=numă rul de indici i egali cu jk, k=1, 2,…,s Atunci se obţine: )l( )a(P) a(P)a(P)(Q= sv22v1vi iia (2 23) s1l21 i )a(P=polinom ortonormal de grad vk în variabila ak deci unde v1+v2+…+vs=l ş kv k E[ )a(P)a(P]=0 dacă vr, 1≤t,r≤n tvrv vt≠ tr Se consideră monoamele liniar independente f0(a), f1(a),… unde a=a1, a2, …,an Printr-o ortonormalizare Gram-Schmidt se obţine: 19 -Teza de doctorat-      = (a)f(a)P 00     (a)f (a)f(a)f   − 01kk     (a)](a)fE[f (a)](a)fE[f(a)](a)fE[f −−−− 01k1k1kk1k   = det(a)P  k       (a)](a)fE[f (a)](a)fE[f(a)](a)fE[f − 001k0k0    (2 24) unde k=1, 2,… Numărul maxim de polinoame Pk(a) este egal cu K-1 )l( in ca produse de polinoame Pk(a) deci combinînd Polinoamele )(Qa se obţ li 2i1i cîte un element de pe următoarele linii: 1   )(aP )(aP )(aP 11-K1211   )(aP )(aP )(aP (2 25)  21-K2221  …   )(aP )(aP )(aP n1-Kn2n1  11 Numărul de produse de cîte 1 polinom Pk(a) este egal cu C)1K(1− n+ 22 Numărul de produse de cîte 2 polinoame Pk(a) este egal cu C)1K(− n … nn Numărul de produse de cîte n polinoame Pk(a) este egal cu C)1K(− n )l(n Numărul total de polinoame )(Qa este egal cu K (conform formulei lui Newton) li 2i1i i un număr finit Corolarul 2 1-Dacă fi(a)=a atunci calculul polinoamelor Pk(a) implică de medii statistice ale variabilelor aleatoare a i Corolarul 2 2-Dacă fi(a)=a atunci toate polinoamele Pk(a) au media statistică 0 cu excepţia polinomului P0(a) care are media statistică 1 =0k,1  E[Pk(a)]=E[Pk(a)P0(a)]= (2 26)  −≤≤1Kk1,0  )l( i cu formula: Corolarul 2 3-Coeficienţii k pot fi calculaţ i ii l21 )l()l( k=E[f(a))(Qa] (2 27) i iili 2i1i l21 20 -Teza de doctorat- Forma matricială a seriei ortonormale Voterra este: f=Qk (2 28) unde: n f=vector coloană K dimensional format din valorile funcţ iei f(a): 1)f(1,1, ,     − 1,1, ,1)f(  (2 29) f=     −1)f(1,1, ,   n)l( k=vector coloană K dimensional format din valorile coeficienţ ilor k: i ii l21 (0)   k   (1) k   1   k=(1) (2 30) k   2     (n) k   123 n   nn)l( formată din valorile polinoamelor )(Qa Q=matrice K x K dimensională li 2i1i )l( Datorită ortonormalităţii polinoamelor )(Qa există relaţia: li 2i1i Tn QQ=KI (2 31) nn unde I este matricea identitate K x K dimensională T ine : Dacă se înmulţeşte relaţia (2 28) la stînga cu Q se obţ -nT k=KQf (2 32) În concluzie vectorul k este o transformată a vectorului f Tipul transformatei depinde de matricea Q adică de distribuţia variabilelor aleatoare a În general sn, vectorul semnal la momentul de timp discret n, depinde de m biţ i de intrare curenţi (a1, a2,…, am) ş i v biţi de intrare anteriori (am+, am+,…, am+) ai 12v transmiţătorului TCM (Fig 2 6): 21 -Teza de doctorat- = n21n)a, ,a,f(as= nnnnnn (2 33) ++++= ∑∑ ∑∑∑∑ n2112 nljiijljiijii0 aaac aaacaacacc+ 1i1j1l1i1j1i====== ljiji 0 este parametrul Chernoff care nu depinde de i Dacă se consideră un canal de comunicaţie fără memorie (memoryless communications channel) atunci: L −⊕ ∑ L)]]f(,m[)]f(,[m[λ iiiiicrecr −⊕ )]]f(,m[)]f(,λ[m[crecr = iiiii1i (3 27) = ∏ ii}E{e}E{ecc = 1i Se defineşte: ⊕ )]]f(,m[)]f(,λ[m[crecr− (3 28) }E{e)]}f(),[f(∆exp{iiiiicecc= ⊕− iiiiλ unde parametrul λ ar trebui ales astfel încît să minimizeze limita Chernoff (3 26) Astfel limita pentru P[EL] probabilitatea medie a unei erori de lungime L cauzată de secvenţa de vectori eroare EL devine: ∆]exp{P[]W[]P[ECCCEE⊕ )]}f(),[f( =≤ (3 29) ∑ LLLλLLL− LC Dacă se presupune că simbolurile sursei de date ai sunt echiprobabile 1 -m′ rul de biţi folosiţi la codarea simbolurilor sursei) se (P(ai)==2 unde m' este numă M′ definesc G(ei) matricile ponderilor vectorilor eroare ei (matrici NxN dimensionale) i coloana q este : astfel încît [G(ei)]p,q elementul de pe lina p ş 1 [G(ei)]p,q=⊕ − ∆exp{ecc )]}f(),[f( (3 30) ∑ →→iqpqpλ ′ M qpc → unde D este o nedeterminată care poate fi considerată ca operaţia de întîrzîiere cu o iei din starea p unitate de timp şi cp este vectorul cuvînt de cod corespunză →qtor tranziţ în stare q (suma apare din cauza tranziţiilor paralele posibile între starea p şi starea q) Prezenţa biţilor necodaţi convoluţional în transmiţătorul TCM determină apariţia tranziţiilor paralele în diagrama trellis 37 -Teza de doctorat- este Concluzia 3 12-Pentru modelul canalului de comunicaţie general, parametrul ∆ λ 2 echivalent cu d, pă tratul distanţei euclidiene minime dintre două secvenţe posibile min de vectori semnal din setul Ω (constelaţia de semnale extinsă), pentru modelul canalului de comunicaţie cu zgomot aditiv, alb, gaussian 3 2 Limita inferioară a probabilităţii erorii de secvenţă Se consideră cazul unui canal de comunicaţie cu zgomot aditiv, alb, gaussian Calculul limitei superioare a probabilităţii erorii de secvenţă (de simbol) se bazează pe faptul că P(e) a oricărui receptor TCM real este mai mare decît P(e) a oricărui receptor TCM ideal care utilizează informaţiile furnizate de un observator Receptorul TCM ideal operează astfel: observatorul analizează o secvenţă lungă de vectori semnal, sau echivalent o secvenţă lungă de cuvinte de cod transmise CK=(ci, ci+1,…,ci+K-1) ş i îi comunică receptorului TCM ideal că secvenţa de vectori cuvinte de C′ are cea mai mica distanţă cod transmise a fost C′ K=(c′i, c′i+1,…,c′i+K-1), dacăK euclidiană fată de CK (nu neapă rat egală cu dmi = distanţa euclidiană minimă dintre n două secvenţe posibile de vectori semnal din setul Ω, constelaţia de semnale extinsă) Probabilitatea erorii de secvenţă (de simbol) pentru receptorul TCM ideal este : CC d1)(f)(f1′ KK− min )(erfc][P][I)(erfc][P)e(P≥= ∑∑ KKKidealCCC (3 31) N22N22 00CC LL unde: ′ −=)(f)(fmind pt,1CC  KKmin  I[CK]=KC (3 32)   rest,0  Concluzia 3 13-Probabilitatea erorii de secvenţă (de simbol) pentru receptorul TCM real este limitată inferior astfel : d1 min (3 33) Ψ )erfc( P(e)≥ N22 0 unde: 1 ][P][I=ψ ∑ KKCC este probabilitatea ca la orice moment de timp discret, o 2 LC traiectorie aleasă aleator prin diagrama trellis, să aibă o altă traiectorie, care diverge de la ea la acel moment de timp discret şi converge la ea mai tîrziu, astfel încît distanţa euclidiană dintre ele să fie egală cu dmi n 1=ψ dacă pentru orice traiectorie prin diagrama trellis, la orice moment de timp discret, există o altă traiectorie, care diverge de la ea la acel moment de timp discret şi 38 -Teza de doctorat- converge la ea mai tîrziu, astfel încît distanţa euclidiană dintre ele să fie egală cu dmi n (în general această proprietate nu este adevărată) ψ este ca toate traiectoriile din diagrama O condiţie suficientă pentru ca 1= trellis să contribuie în mod egal la P(e) O condiţie suficientă pentru ca 1=ψ este ca G să aibă toate elementele egale O condiţie suficientă pentru ca 1=ψ este ca aibă toate elementele egale toate matricile G(ei) să Se notează: m′ = numărul de biţi folosiţi la codarea simbolurilor sursei, Λ = lungimea minimă a evenimentului eroare care produce dmi, n (dmi) = numărul de traiectorii de lungime Λ care au o traiectorie cu care se întîlnesc NΛ n la distanţa dmi n Cu aceste notaţii rezultă: m′Λ rul de traiectorii de lungime Λ care încep din orice stare, N2 = numă Λmin)d(N = este probabilitatea ca la orice moment de timp discret, o traiectorie ′ Λm 2N aleasă aleator prin diagrama trellis, să aibă o altă traiectorie, care diverge de la ea la acel moment de timp discret şi converge la ea mai tîrziu, astfel încît distanţa euclidiană fie Λ dintre ele să fie egală cu dmi ş ni lungimea lor să Astfel se obţine : Λ min)d(N ≥Ψ (3 34) ′ Λm 2N Concluzia 3 13-Probabilitatea erorii de secvenţă (de simbol) pentru receptorul TCM real este limitată inferior astfel : d)d(1N minminΛ (3 35) )(erfc)e(P≥ ′ N22mΛ 2N0 3 2 1 Limita inferioară a probabilităţii erorii de bit Concluzia 3 14-Probabilitatea erorii de bit este limitată inferior astfel: d)d(1d1N minminminΛ (3 36) ≥ )(erfc)(erfcP= bΨ ′ ′ N2m2N2m2mΛ 2N00′ 3 3 Calculul funcţiei de transfer scalară a diagramei de stare a erorii Funcţia de transfer scalară a diagramei de stare a erorii T(D,I) se calculează astfel: ♦ din diagrama de stare a erorii se calculează funcţia de transfer scalară T(D,I) prin metode bazate pe soluţia unui sistem de ecuaţii liniare sau prin tehnici de 39 -Teza de doctorat- reducere a grafurilor, ♦ se calculează matricile G(ei) a ponderilor vectorilor eroare ei (sau în cazul unei scheme TCM uniforme se calculează profilurile W(e) a ponderilor vectorilor eroare e), ♦ se înlocuiesc etichetele formale ale funcţiei de transfer T(D,I) cu matricile G(ei) a ponderilor vectorilor eroare ei sau cu profilurile W(e) a ponderilor vectorilor ε eroare e după ce în prealabil s-au înmulţit cu I 3 4 Calculul distanţei euclidiene minime 3 4 1 Calculul distanţei euclidiene minime utilizînd diagrama de stare a erorii Se consideră G(ei) matricile ponderilor vectorilor eroare ei În general elementele lor sunt sume de termeni generaţi de tranziţii paralele Aceste tranziţii paralele apar datorită prezenţei biţilor necodaţi convoluţional în transmiţătorul TCM Pentru calculul dmi, este suficient să m pentru elementele matricilor G(ei) a n consideră ponderilor vectorilor eroare ei, în locul sumelor de termeni, doar un singur termen ş i anume acela cu exponentul mai mic (întrucît ne interesează calculul lui dmi) n Funcţia de transfer scalară T(D) a diagramei de stare a erorii poate fi scrisă ca o serie de genul : 22 dd )DN(d)DN(dT(D)nextmin+ = (3 37) nextmin+ unde dnex este după toarea distanţă euclidiană dintre două traiectorii din t dmin urmă diagrama trellis care formează un eveniment eroare Se definesc funcţiile Φ i Φ(D) : 1(D) ş2  )eD(T =Φln)D( (3 38) 1  )D(T  )]D(Tln[ (3 39) )D(= 2Φ Dln Se arată că funcţiile Φ(D) şi Φ(D) au următoarele propietăţi : 12 2 ⇒ funcţia Φ te monoton către d cînd D→0, deci pentru orice D>0 1(D) descreş min 2 superioară pentru d, funcţia )D(Φ 1 va furniza o limită min 2 te monoton către d cînd D→0, deci pentru orice D>0 ⇒ funcţia Φ 2(D) creş min 2 inferioară pentru d funcţia )D(Φ 2 va furniza o limită min Concluzia 3 15-Dacă considerăm două secvenţe de valori pozitive, descrescătoare i (D)Φ furnizează două către zero ale nedeterminatei D, atunci funcţiile )D(Φ 1 ş2 2 secvenţe de valori care converg către d min 40 -Teza de doctorat- În cazul schemelor TCM uniforme, funcţia de transfer scalară T(D) a diagramei de stare a erorii poate fi calculată utilizînd etichete de ramură scalare în diagrama de stare a erorii Aceste etichete de ramură scalare sunt egale cu suma elementelor de pe orice linie sau coloană a matricilor G(ei) a ponderilor vectorilor eroare ei Pentru calculul dmi, este suficient să m pentru etichetele de ramură scalare din n consideră diagrama de stare a erorii, în locul sumelor de termeni, doar un singur termen şi anume acela cu exponentul mai mic (întrucît ne interesează calculul lui dmi) n Se defineşte w(e) ponderea euclidiană a vectorului eroare e : 2 )(f)(fmin)(wecce⊕−= (3 40) Cc∈ 0 Concluzia 3 16-În cazul schemelor TCM uniforme, funcţia de transfer scalară T(D) a diagramei de stare a erorii poate fi calculată utilizînd etichete de ramură scalare în w(e) diagrama de stare a erorii pentru care se consideră valorile D 3 4 2 Algoritmul Saxona-Mulligan-Wilson Pentru fiecare pereche de ramuri dintr-o secţiune a diagramei trellis se defineşte o distanţă euclidiană între vectorii semnal care etichetează ramurile respective Dacă există tranziţii paralele atunci fiecare ramură va fi asociată cu o subconstelaţie de semnale În acest caz va fi utilizată distanţa euclidiană minimă dintre oricare doi vectori semnal extraşi din cele două subconstelaţii de semnale Pătratul distanţei euclidiene dintre două secvenţe de vectori semnal asociate cu două traiectorii din diagrama trellis se obţine prin suma pătratelor distanţelor euclidiene individuale (vezi Fig 2 2) Algoritmul Saxona-Mulligan-Wilson se bazează pe actualizarea )n()n()n( elementelor matricii D, unde termenul general ][δ D = minimul pqpq= pătratelor distanţelor euclidiene dintre toate perechile de traiectorii din diagrama i care ajung în stările p şi q la acelaşi trellis, care diverg dintr-o stare iniţială σ 0 ş moment de timp discret (Fig 3 1) )n( şi elementele ei de pe diagonala Concluzia 3 17-Matricea D este simetrică )n( principală sunt egale cu δ tratelor distanţelor euclidiene dintre toate pp=minimul pă i care perechile de traiectorii din diagrama trellis care diverg dintr-o stare iniţială σ 0 ş ajung în aceeaşi stare p la acelaşi moment de timp discret 41 -Teza de doctorat- Fig 3 1 Două perechi de traiectorii din diagrama trellis care diverg dintr-o stare iniţială σ i care ajung în stările p şi q la acelaşi moment de timp discret 0 ş Algoritmul Saxona-Mulligan-Wilson se desfăşoară astfel: Etapa 1-Pentru fiecare stare p, se găsesc cele M stări predecesoare din care sunt posibile tranziţii către starea p şi se memorează aceste M stări predecesoare într-un )0( tabel Se setează 1− δ pentru toate p şi q≥p pq= Etapa 2-Pentru fiecare pereche de stări (p, q), se găseşte minimul pătratelor distanţelor euclidiene dintre toate perechile de traiectorii din diagrama trellis care diverg dintr-o i care ajung în stările p şi q într-o unitate timp discret şi se memorează stare iniţială σ 0 ş )1( ca δ pq (Fig 3 2) Etapa 3-Pentru fiecare pereche de stări (p, q) se găsesc în tabelele memorate la Etapa 2 i q1, q2,…,qM (Fig 3 3) În general există M 1 cele M stări predecesoare: p1, p2,…,pM ş traiectorii posibile la momentul de timp discret n-1 care trec prin stările p şi q la momentul de timp discret n Aceste traiectorii trec prin perechile de stări: … ,(p ),q ,(p ,(p,),q )q  M12111  … ,(p ),q ,(p ,(p,),q )q  M22212 (3 41)  ……………………………   … ,(p ),q ,(p ,(p,),q )q  MM2M1M (n) Se calculează δ pq astfel: 42 -Teza de doctorat- 21)(n(n)− →→+= q)qp,(pdmin{δδ 11pq+ qp 11 21)(n− →+ q)qp,(pdδ+ 21→ qp 21 21)(n− →+ q)qp,(pdδ+ M1→ qp M1 (3 42) 21)(n− →+ q)qp,(pdδ+ 1M→ qp 1M 21)(n− →+ q)qp,(pdδ+ 2M→ qp 2M 21)(n− + →q)}qp,(pdδ MM→ qp MM )1n()1n(− − 2 sunt calculate din etapa anterioară şi )qq,pp(d→→ sunt δ unde , ,δ qpqpji MM11 pătratele distanţelor euclidiene dintre cei doi vectori semnal asociaţi cu tranziţiile pp→qq→ i, unde i, j = 1, 2, …,M işj )n()n( δ atunci se Etapa 4-Dacă există cel puţin o pereche de stări astfel încît minδ pppq R0 3) fiecare grup să fie invariant la o rotaţie de 0°, 90°, 180°, 270° Astfel dacă un punct de semnal 2-D dintr-un anumit grup este rotit cu 0°, 90°, 180°, 270° trebuie să se obţină un alt punct de semnal 2-D din acelaşi grup Această cerinţă permite demodulatorului să se caleze pe fază doar într-un interval de 90° ceea ce implică timpi de sincronizare mai mici De asemenea achiziţia de semnal poate fi făcută fără o 100 -Teza de doctorat- perioadă de antrenare (sau secvenţe de antrenare) Construcţia unei constelaţii de semnale 2N-D rectangulare (subset finit al latticei 2N Z) se face prin concatenarea a N puncte de semnal 2-D constituente ş i eliminînd acele puncte de semnal 2N-D rezultante care conţin mai mult decît un punct de semnal 2-D din OG Construcţia unei constelaţii de semnale 2N-D nerectangulare (subset finit al 2N2N2N lattticei D, E, DE) se face prin concatenarea a N puncte de semnal 2-D constituente şi eliminînd acele puncte de semnal 2N-D rezultante care conţin mai mult decît un punct de semnal 2-D din OG sau care nu sunt puncte valide ale latticii nerectangulare În continuare am prezentat construcţia unei constelaţii de semnale 2N-D rectangulare Se construiesc puncte de semnal 2N-D prin concatenarea a N puncte de semnal 2-D şi se ţine cont de următoarele observaţii : • numărul punctelor de semnal 2N-D care conţin N puncte de semnal 2-D din IG NQ este 2, • numărul punctelor de semnal 2N-D care conţin N-1 puncte de semnal 2-D din NQ IG şi un punct de semnal 2-D din OG în una din primele N poziţii este βN2 NQ+1NQNQ Deoarece în total sunt necesare 2=2+2 puncte de semnal 2N-D, se presupune că numărul punctelor de semnal 2N-D care conţin N-1 puncte de semnal NQ 2-D din IG şi un punct de semnal 2-D din OG este mai mare decît 2 care este numărul necesar de puncte de semnal 2-D : NQNQ1 2⇔β≥ βN2≥ (4 92) N Se presupune că fiecare punct de semnal 2-D al constelaţiei de semnale 2-D, prezintă o regiune Voronoi de arie unitară : S(Voronoi)=1 (4 93) Se presupune că în constelaţia de semnale 2-D există un număr suficient de mare de puncte de semnal 2-D astfel încît Eme a constelaţ poate fi dieiei de semnale 2-D să aproximată cu integrala pe suprafaţa cercului C2(O, R2) care o încadrează 4 4 1 2 1 Extinderea simplă a constelaţiei de semnale 2-D Am urmărit demonstrarea formulelor pentru numărul de puncte de semnal 2-D, energia medie şi parametrii CER şi PAR ai constelaţiei de semnale 2-D extinsă simplu Se construiesc puncte de semnal 2N-D prin concatenarea a N puncte de semnal 2-D şi se ţine cont de următoarele observaţii : • numărul punctelor de semnal 2N-D care conţin N puncte de semnal 2-D din IG NQ este 2, • numărul punctelor de semnal 2N-D care conţin N-1 puncte de semnal 2-D din 1 ii este IG şi un punct de semnal 2-D din OG în una din primele poziţ β 101 -Teza de doctorat- 1 NQNQ β2=2 β NQ+1NQNQ Rezultă în total 2=2+2 puncte de semnal 2N-D Lema 4 2-O constelaţie de semnale 2-D, încadrată de un cerc C(O, ρ) de rază ρ şi cu 2 2ρ puncte de semnal 2-D şi are energia medie Emedie= centrul în origine, conţine πρ 2 Fig 4 9 Constelaţia de semnale 2-D obţinută prin extindere simplă Lema 4 3-O constelaţie de semnale 2-D, încadrată de două cercuri C1(O, ρ i 1) ş 22 C2(O, ρ i cu centrul în origine, conţine )ρπ(ρ− puncte de semnal 2) de raze ρ1 0= Partiţionarea constelaţiei de semnale 2N-D în subseturi 2N-D am realizat-o diferit pentru cazul în care N este par, respectiv impar Dacă N este par atunci: 1 Am partiţionat fiecare din cele 2 constelaţii de semnale N-D constituente cu 2 în familii N-D, subfamilii N-D şi subseturi N-D cu MSED din ce în ce 0dMSED= mai mare Fiecare nivel de partiţionare a constelaţiei de semnale N-D dublează MSED (partiţionare binară) Am continuat cu partiţionarea pînă cînd am obţinut subseturi 22 N-D cu ddMSED> 0= După primul nivel de partiţionare a constelaţiei de semnale N-D am obţinut familii N-D iar după următoarele nivele de partiţionare (cu excepţia ultimului) am obţinut subfamilii N-D După ultimul nivel de partiţionare am obţinut subseturi N-D Cu cît este mai fină partiţionarea constelaţiei de semnale N-D cu atît rezultă mai multe subseturi N-D şi cu MSED mai mare 22 2 Am construit tipuri 2N-D cu ddMSED> prin concatenarea a cîte 2 subseturi 0= N-D 22s 3 Am construit subseturi 2N-D cu ddMSED> prin gruparea a cîte S=2 tipuri 0= 2N-D, unde S este numărul subseturilor N-D conţinute într-o subfamilie N-D Gruparea tipurilor 2N-D în subseturi 2N-D trebuie făcută astfel încît să fie îndeplinite două cerinţe : ¾ construcţia codului TCM 2N-D să fie cît mai simplă posibil (să se obţină un număr cît mai mic de posibil de subseturi 2N-D adică selecţia unui subset 2N-D să se facă cu un număr cît mai mic posibil de biţi codaţi de codorul convoluţional), ¾ fiecare subset 2N-D să fie invariant rotaţional la cît mai multe rotaţii posibile (ambiguităţi de fază ale constelaţiei de semnale 2N-D) Dacă nu este posibil ca fiecare subset 2N-D să fie invariant la toate rotaţiile, atunci ar trebui cel puţin ca orice rotaţie să transforme un subset 2N-D în alt subset 2N-D 4 Am grupat subseturile 2N-D în subfamilii 2N-D şi familii 2N-D cu MSED din ce în ce mai mică Fiecare nivel de grupare a subseturilor 2N-D înjumătăţeşte MSED Am 2 continuat cu gruparea pînă cînd am obţinut două familii 2N-D cu d2MSED= i 0 ş notate F0 ş i F1 112 -Teza de doctorat- Dacă N este impar atunci: 1 Am partiţionat fiecare din cele N constelaţii de semnale 2-D constituente cu 2 în familii 2-D, subfamilii 2-D şi subseturi 2-D cu MSED din ce în ce 0dMSED= mai mare Fiecare nivel de partiţionare a constelaţiei de semnale 2-D dublează MSED (partiţionare binară) Am continuat cu partiţionarea pînă cînd am obţinut subseturi 22 2-D cu ddMSED> 0= După primul nivel de partiţionare a constelaţiei de semnale 2-D am obţinut familii 2-D iar după următoarele nivele de partiţionare cu excepţia ultimului am obţinut subfamilii 2-D După ultimul nivel de partiţionare am obţinut subseturi 2-D Cu cît este mai fină partiţionarea constelaţiei de semnale 2-D cu atît rezultă mai multe subseturi 2-D şi cu MSED mai mare 22 2 Am construit tipuri 2N-D cu ddMSED> prin concatenarea a cîte N subseturi 0= 2-D 22sNu-m-1 3 Am construit subseturi 2N-D cu ddMSED> prin gruparea a cîte S=2=2 0= u tipuri 2N-D, unde U=2 este numă rul de subseturi 2-D ale constelaţiei de semnale 2-D constituentă Gruparea tipurilor 2N-D în subseturi 2N-D trebuie făcută astfel încît să fie îndeplinite cele două cerinţe de mai sus : ¾ construcţia codului TCM 2N-D să fie cît mai simplă posibil, ¾ fiecare subset 2N-D să fie invariant rotaţional la cît mai multe rotaţii posibile 4 Am grupat subseturile 2N-D în subfamilii 2N-D şi familii 2N-D cu MSED din ce în ce mai mică Fiecare nivel de grupare a subseturilor 2N-D înjumătăţeşte MSED Am 2 continuat cu gruparea pînă cînd am obţinut două familii 2N-D cu d2MSED= i 0 ş notate F0 ş i F1 4 4 3 Asignarea subseturilor 2N-D la tranziţiile de stare ale codului TCM 2N-D Asignarea subseturilor 2N-D la tranziţiile de stare ale codului TCM 2N-D se face astfel încît să fie îndeplinite trei cerinţe : 1 Subseturile 2N-D, asociate cu tranziţiile care pornesc dintr-o anumită stare, să difere unele de altele şi să aparţină aceleaşi familii 2N-D Acelaşi considerent se aplică şi pentru subseturile 2N-D asociate cu tranziţiile care ajung într-o anumită stare 2 MSED dintre două secvenţe de subseturi 2N-D, care corespund la două traiectorii 22 distincte prin diagrama trellis, să fie mai mare sau egală cu dd= , care este free MSED a fiecărui subset 2N-D (MSED de intersubset 2N-D să fie mai mare sau cel puţin egală cu MSED de intrasubset 2N-D) Această cerinţă garantează că MSED dintre două secvenţe de puncte de semnal 2 i elimină evenimentele eroare care diferă cu mai mult decît un punct de 2N-D este d ş semnal 2N-D faţă de secvenţa corectă de puncte de semnal 2N-D 113 -Teza de doctorat- Dacă se neglijează efectul de margine, atunci coeficientul de eroare Nfree al codului TCM 2N-D este egal cu numărul vecinilor cei mai apropiaţi ai unui punct de semnal 2N-D şi aflaţi într-un subset 2N-D cu un număr finit de puncte de semnal 2N-D Dacă nu se neglijează efectul de margine, atunci coeficientul de eroare al codului TCM 2N-D este mai mic decît cel din definiţia anterioară 3 Dacă X este subsetul 2N-D asociat cu tranziţiile din starea curentă σ n în starea următoare σ şi dacă X1, X2 ,X3 sunt subseturile 2N-D obţinute prin rotirea lui X Nn+ între cu 90°, 180°, 270° atunci se pot defini trei funcţii F1, F2, F3 de corespondenţă stările codului TCM 2N-D astfel încît : fie asociat cu tranziţiile din starea curentă )(σ în starea următoare ⇒ X1 să F n1 )(σF, Nn1+ fie asociat cu tranziţiile din starea curentă )(σ în starea următoare ⇒ X2 să F n2 )(σF, Nn2+ fie asociat cu tranziţiile din starea curentă )(σ în starea următoare ⇒ X3 să F n3 )(σF Nn3+ Acestă cerinţă garantează că vom obţine un cod TCM 2N-D transparent la toate ambiguităţile de fază ale constelaţiei de semnale 2N-D, cu alte cuvinte un cod TCM 2N-D invariant rotaţional la o rotaţie de 0°, 90°, 180°, 270° 4 4 4 Realizarea corespondenţei dintre cei NQ+1 biţi şi constelaţia de semnale 2N-D Corespondenţa complicată dintre cei NQ+1 biţi (primii m+1 biţi codaţi cu un codor convoluţional şi ceilalţi NQ-m biţi necodaţi) şi constelaţia de semnale 2N-D am convertit-o în N corespondenţe simple ale constelaţiilor de semnale 2-D constituente cu ajutorul unui convertor de bit şi a unui codor bloc Convertorul de bit şi codorul bloc transformă cei NQ+1 biţi în N grupuri de cîte Q+1 biţi fiecare, adică în NQ+N biţi Fiecare grup de Q+1 biţi adresează un bloc de corespondenţă 2-D şi selectează un punct de semnal 2-D Astfel sunt selectate N puncte de semnale 2-D şi prin concatenarea lor se obţine un punct de semnal 2N-D corespunzător celor NQ+1 biţi 4 5 Generalizarea construcţiei Wei a codorului TCM 2N-D Schema bloc codorului TCM 2N-D (construcţia propusă) am prezentat-o în Fig 4 11 NQ+1 Am utilizat o constelaţie de semnale 2N-D cu 2 puncte de semnal 2N-D Construcţia constelaţiei de semnale 2N-D şi partiţionarea ei în subseturi 2N-D le-am realizat ca şi în paragrafele anterioare 4 4 1 şi 4 4 2 Primii m biţi de la intrare sunt codaţi cu un codor convoluţional cu rata de m v codare ş ri i cu V=2 stă 1m+ 114 -Teza de doctorat- Am notat starea curentă şi starea următoare a codorului convoluţional astfel: W1,p, W2,p ,…, Wv,p (4 126) unde: p=n este starea curentă a codorului convoluţional, p=n+N este starea următoare a codorului convoluţional Valoarea zecimală a stării (W1,p, W2,p ,…, Wv,p) a codorului convoluţ ional este egală cu: 11v− +++= (4 127) pv,p1,vp1,pWW2 W2σ− unde: p=n este starea curentă codorului convoluţional, p=n+N este starea următoare codorului convoluţional Am presupus că numărul de subseturi 2N-D dintr-o familie 2N-D este egal cu F, unde este necesar ca F≤V Setul stărilor curente {} a codorului convoluţional le-am partiţionat în F nσ subseturi: jFi−≤≤−≤≤+ 1,Fi0 1Fj0 (4 128) {} a codorului convoluţional le-am partiţionat în F Setul stărilor următoare {} Nnσ + subseturi: kFj−≤≤−≤≤+ 1,Fj0 1Fk0 (4 129) {} Conectivitatea diagramei trellis a codorului TCM 2N-D am realizat-o astfel: pentru orice 0≤j≤F-1 există tranziţii de stare între fiecare stare curentă Fi+j, 0≤i≤F-1 şi fiecare stare următoare Fj+k, 0≤k≤F-1 Am determinat rotaţia subseturilor 2N-D la rotaţii sucesive cu 90° a subseturilor 2-D şi rotaţia stărilor codorului convoluţional la rotaţii sucesive cu 90° a subseturilor 2N-D Din rotaţia stărilor codorului convoluţional am dedus cele trei funcţii de i F3 corespondenţă F1, F2 ş Asignarea subseturilor 2N-D la tranziţiile de stare ale codului TCM 2N-D am realizat-o ca şi în paragraful 4 4 3 Astfel: ⇒ la tranziţiile care pornesc din stări pare (Y0,n=Wv,n=0) sunt asignate subseturi 2N-D din familia 2N-D notată F0, ⇒ la tranziţiile care pornesc din stări impare (Y0,n=Wv,n=1) sunt asignate subseturi 2N-D din familia 2N-D notată F1 Am folosit un codor convoluţional cu reacţie cu conexiune înapoi (systematic feedback convolutional encoder) Biţii de intrare I1,n, I2,n,…,Im, sunt conectaţ ni direct la ieşire Bitul Y0,n este codat convoluţ ional şi depinde de biţii de intrare I1,n, I2,n,…,Im, n şi de biţii de stare curentă W1,n, W2,n ,…, Wv,n 115 -Teza de doctorat- Din diagrama trellis a codorului TCM 2N-D am dedus relaţiile pentru biţii de stare următoare W1,n+N, W2,n+N,…, Wv,n+N în funcţ ie de biţii de intrare I1,n, I2,n,…,Im, şi n de biţii de stare curentă W1,n, W2,n ,…, Wv, n Cei m+1 biţi de la ieşirea codorului convoluţional selectează subsetul 2N-D care va fi utilizat Ceilalţi NQ-m biţi de la intrare, rămaşi necodaţi convoluţional, selectează cu ajutorul unui codor bloc şi a unui convertor de bit, punctul de semnal 2N-D care va fi transmis din subsetul 2N-D selectat anterior u Am notat cu U=2 numă rul de subseturi 2-D ale constelaţiei de semnale 2-D, constituentă a constelaţiei de semnale 2N-D Dacă N=par, la intrarea convertorului de bit: ⇒ primii m+1 biţi selectează subsetul 2N-D utilizat, ⇒ următorii s biţi selectează tipul 2N-D în cadrul subsetului 2N-D utilizat, ⇒ următorii Nu-m-1-s biţi selectează tipul N-D, pentru fiecare din cele 2 subseturi N-D constituente ale tipului 2N-D şi aşa mai departe pînă la selecţia celor N subseturi 2-D constituente Dacă N=impar la intrarea convertorului de bit: ⇒ primii m+1 biţi selectează subsetul 2N-D utilizat, ⇒ următorii s=Nu-m-1 biţi selectează tipul 2N-D, în cadrul subsetului 2N-D utilizat şi implicit cele N subseturi 2-D constituente Am notat cu r = restul împărţirii lui Nu la Q şi cu c = cîtul împărţirii lui Nu la Q, te cu IQ,n+c-1 unde este necesar ca r > 1 Dacă r ≤ 1 atunci Ir-1,n+c se înlocuieş Convertorul de bit converteşte cei Nu biţi de la intrarea sa Y0,n, I1,n, I2,n,…, Ir-1,n+c în N grupuri de cîte u biţ i fiecare: Z0,p, Z1,p ,…, Zu-1,p, p=n, n+1, …, n+N-1 (4 130) Aceste N grupuri de cîte u biţi fiecare selectează N subseturi 2-D, corespunzătoare tipului 2N-D selectat în cadrul subsetului 2N-D utilizat Pentru a obţine un cod TCM 2N-D transparent la orice ambiguitate de fază a constelaţiei de semnale 2N-D trebuie sa fie îndeplinită cerinţa detailată în continuare Dacă X este tipul 2N-D asociat cu structura de biţi ′′ atunci se notează cu X1, X2, ,I,Y++++− I, ,I , ,I ,I cn1,-rcn2,-rbna,bn1,an2,n1,n0,I X3 tipurile 2N-D obţ inute prin rotirea lui X cu 90°, 180°, 270°, unde 1ra2−≤≤ şi cb0≤≤ ′′′′′ Dacă I(′ )I , I(, )I I( )I sunt 1bna,bn1,a++−2bna,bn1,a++−3bna,bn1,a++− ′ perechile de biţi obţinute prin translaţia perechii de biţi )II(′ cu una, două bna,bn1,a++− sau trei poziţii în secvenţa circulară 00→10→01→11→00, atunci tipurile 2N-D asociate cu structurile de biţi: ′′ , ,I,Y++++− I(, ,I , ,)I ,I cn1,-rcn2,-r1bna,bn1,an2,n1,n0,I ′′ , ,I,Y++++− I(, ,I , ,)I ,I cn1,-rcn2,-r2bna,bn1,an2,n1,n0,I ′′ , ,I,Y++++− I(, ,I , ,)I ,I cn1,-rcn2,-r3bna,bn1,an2,n1,n0,I sunt X1, X2 ,X3 116 -Teza de doctorat- Fig 4 11 Schema bloc a codorului TCM 2N-D (construcţia propusă) 117 -Teza de doctorat- Pentru a suprima toate cele 4 ambiguităţi de fază a constelaţiei de semnale 2N-D (0°, 90°, 180°, 270°) am folosit un codor diferenţial pe 2 biţi de forma: ′′′′ = (4 131) I(+ I()I II )mod4I bna,bn1,aN-bna,N-bn1,abna,bn1,a++−++−++− Tabelul de adevăr al codorului diferenţial care rezultă din relaţia (4 131) este următorul: Tabelul 4 14 Tabelul de adevăr al codorului diferenţial BiţBiţii de ieşire anteriori ii de ieşire ′′′′ I I I++−bna,bn,1a++− N-bna,N-bn,1aI 00 01 10 11 Biţii de intrare 00 00 01 10 11 I I 01 01 10 11 00 bna,bn1,-a++ 10 10 11 00 01 11 11 00 01 10 Din tabelul de adevăr de mai sus am dedus următoarele relaţii între ieşirile şi intrările codorului diferenţial (Fig 4 12): ′′′ I )I(I)I(I ⊗⊕⊕=  +++++N-bna,bna,N-bn1,-abn1,-abn1,-a  (4 132)  ′′ I )I(I ⊕=  +++N-bna,bna,bna,  Fig 4 12 Codorul diferenţial pe 2 biţi (construcţia propusă) 118 -Teza de doctorat- Ultimii NQ-Nu+1 biţi de la intrare sunt codaţi cu un codor bloc, la ieşirea căruia se obţin N grupuri de cîte Q-u+1 biţi fiecare: Zu,p, Zu+1,p ,…, ZQ,, p=n, n+1, …, n+N-1 (4 133) p Aceste N grupuri, de cîte Q-u+1 biţi fiecare, selectează cele N puncte de semnal 2-D, astfel încît cel mult unul să fie situat în grupul extern OG a constelaţiei de semnale 2-D constituentă a constelaţiei de semnale 2N-D Cîştigul de codare asimptotic al codului TCM 2N-D faţă de o transmisie necodată este egal cu: 2   d   2     1dE medie   10log10logγ = =   1010 [dB] (4 134)   CER22 d d   00     E IG   unde: 22 dd==MSED a subseturilor 2N-D obţinute prin partiţionarea constelaţiei de free semnale 2N-D, 2 iei de semnale 2-D utilizată la o transmisie necodată, 0d =MSED a constelaţ Emedie=energia medie a constelaţ iei de semnale 2-D constituentă a constelaţiei de semnale 2N-D, EIG=energia medie a grupului intern IG a constelaţ iei de semnale 2-D constituentă a constelaţiei de semnale 2N-D, CER=rata de extindere a constelaţiei de semnale 2-D constituentă a constelaţiei de semnale 2N-D (Constellation Extension Rate) Întrucît formulele stabilite de Wei şi Sterian pentru cîştigul de codare asimptotic γ al codului TCM 2N-D faţă de o transmisie necodată, nu ţin cont de parametrii PAR şi CER ai constelaţiei de semnale 2-D constituentă a constelaţiei de semnale 2N-D, am dedus o nouă formulă care ia în calcul parametrii respectivi (relaţia 4 134) În relaţiile folosite de Wei şi Sterian (relaţiile 4 2 şi 4 33), energiile medii ale grupurilor se calculează dificil prin însumarea pătratelor distanţelor dintre originea axelor de coordonate a planului complex şi punctele de semnal 2-D situate în grupurile respective şi prin împărţirea la numărul de puncte de semnal 2-D ale grupurilor Spre deosebire de relaţiile folosite de Wei şi Sterian, relaţia propusă (4 134) permite calculul imediat a cîştigului de codare asimptotic γ al codului TCM 2N-D, după calculul parametrului CER a constelaţiei de semnale 2-D constituentă a constelaţiei de semnale 2N-D 119 -Teza de doctorat- 4 6 Generalizarea construcţiei Wei a decodorului TCM 2N-D Decodorul TCM 2N-D utilizează un algoritm de decodare de probabilitate condiţionată maximă (MLD=Maximum Likelihood Decoding) numit algoritmul Viterbi Decodorul Viterbi utilizează o tehnică de decodare cu decizie hard (hard- decision decoding technique) şi găseşte secvenţa de puncte de semnal 2N-D cele mai apropiate de secvenţa de puncte de semnal 2N-D recepţionate Schema propusă pentru decodorul Viterbi conţine două blocuri: ⇒ blocul de calcul a celor mai apropiate puncte de semnal 2N-D faţă de punctul de semnal 2N-D recepţionat în fiecare subset 2N-D şi a metricilor de subset 2N-D asociate, ⇒ decodorul Viterbi propriu-zis, care analizează traiectoriile din diagrama trellis şi alege traiectoria cu metrica totală minimă Astfel decodorul Viterbi propriu- zis, generează o secvenţă de puncte de semnal 2N-D care sunt cele mai apropiate de secvenţa de puncte de semnal 2N-D recepţionate Dacă N nu este o putere întreagă a lui 2 atunci am folosit diagrama de decodare Viterbi I iar dacă N este o întreagă putere a lui 2 atunci am folosit diagrama de decodare Viterbi II Pentru a decide dacă N este sau nu o putere întreagă a lui 2 am propus organigrama din Fig 4 13 Am citit valoarea lui N şi am atribuit la i valoarea iniţială 1 Dacă N este impar atunci am ales diagrama de decodare Viterbi I NN Dacă N este par atunci am efectuat operaţia Dacă rezultatul operaţiei ii 22 este impar atunci apar două cazuri: ⇒ dacă rezultatul operaţiei este diferit de 1, atunci am ales diagrama de decodare Viterbi I, ⇒ dacă rezultatul operaţiei este egal cu 1, atunci am ales diagrama de decodare Viterbi II NN i am reluat operaţia Dacă rezultatul operaţiei este par atunci am incrementat i ş ii 22 120 -Teza de doctorat- Fig 4 13 Organigrama de alegere a diagramei de decodare Viterbi (construcţia propusă) Pentru diagrama de decodare Viterbi I (Fig 4 14) am executat următoarele operaţii: ¾ am descompus punctul de semnal 2N-D recepţionat în N puncte de semnal 2-D constituente, ¾ am găsit în fiecare subset 2-D, cel mai apropiat punct de semnal 2-D faţă de fiecare din cele N puncte de semnal 2-D recepţionate şi am calculat metrica de subset 2-D asociată, ¾ pentru fiecare tip 2N-D, am calculat metrica de tip 2N-D ca suma a N metrici de subset 2-D, ¾ pentru fiecare subset 2N-D, am comparat cele S metrici de tip 2N-D şi metrica cea mai mică devine metrica de subset 2N-D asociată cu subsetul 2N-D respectiv Punctul de semnal 2N-D corespunzător tipului 2N-D cu cea mai mică metrică de tip 2N-D devine cel mai apropiat punct de semnal 2N-D faţă de punctul de semnal 2N-D recepţionat, ¾ am extins traiectoriile din diagrama trellis şi am generat decizia finală asupra punctului de semnal 2N-D recepţionat 121 -Teza de doctorat- Pentru diagrama de decodare Viterbi II (Fig 4 15) am executat următoarele operaţii: ¾ am descompus punctul de semnal 2N-D recepţionat în N puncte de semnal 2-D constituente, ¾ am găsit în fiecare subset 2-D, cel mai apropiat punct de semnal 2-D faţă de fiecare din cele N puncte de semnal 2-D recepţionate şi am calculat metrica de subset 2-D asociată, ¾ pentru fiecare tip 4-D, am calculat metrica de tip 4-D ca suma a două metrici de subset 2-D, ¾ pentru fiecare subset 4-D, am comparat cele două metrici de tip 4-D şi metrica cea mai mică devine metrica de subset 4-D asociată cu subsetul 4-D respectiv Punctul de semnal 4-D corespunzător tipului 4-D cu cea mai mică metrică de tip 4-D devine cel mai apropiat punct de semnal 4-D faţă de primul punct de semnal 4-D recepţionat (primul şi al doilea punct de semnal 2-D recepţionate) N Am repetat procedeul pentru toate cele puncte de semnal 4-D recepţ ionate, 2 ¾ pentru fiecare tip 2N-D, am calculat metrica de tip 2N-D ca suma a două metrici de subset N-D, ¾ pentru fiecare subset 2N-D, am comparat cele S metrici de tip 2N-D şi metrica cea mai mică devine metrica de subset 2N-D asociată cu subsetul 2N-D respectiv Punctul de semnal 2N-D corespunzător tipului 2N-D cu cea mai mică metrică de tip 2N-D devine cel mai apropiat punct de semnal 2N-D faţă de punctul de semnal 2N-D recepţionat, ¾ am extins traiectoriile din diagrama trellis şi am generat decizia finală asupra punctului de semnal 2N-D recepţionat Dacă ţinem cont de efectul de margine al unei constelaţii de semnale 2N-D cu un număr finit de puncte de semnal 2N-D, este necesar ca decizia finală pentru punctul de semnal 2N-D recepţionat, să corespundă unui punct de semnal valid al constelaţiei de semnale 2N-D, adică să conţină cel mult un punct de semnal 2-D din grupul extern OG 122 -Teza de doctorat- Fig 4 14 Diagrama de decodare Viterbi I (N nu este o putere întreagă a lui 2) 123 -Teza de doctorat- Fig 4 15 Diagrama de decodare Viterbi II (N este o putere întreagă a lui 2) 124 -Teza de doctorat- În continuare, decodorul TCM 2N-D conţine N blocuri de corespondenţă 2-D, un convertor de bit, un decodor diferenţial şi un decodor bloc care realizează operaţiile inverse de la codare Schema bloc a decodorului TCM 2N-D (construcţia propusă) am prezentat-o în Fig 4 16 Am prelucrat relaţiile (4 132) şi am obţinut următoarele relaţii între ieşirile şi intrările decodorului diferenţial (Fig 4 17):  ′′′′′ ⊗⊕⊕⊕= Iˆ ]Iˆ)IˆIˆ([)IˆIˆ( ++++++N-bna,N-bna,bna,N-bn1,-abn1,-abn1,-a   ′′ ⊕= Iˆ )IˆIˆ(  N-bna,bna,bna, +++  (4 135) Tabelul de adevăr al decodorului diferenţial care rezultă din relaţia (4 135) este următorul: Tabelul 4 15 Tabelul de adevăr al decodorului diferenţial ii de ieşire BiţBiţii de intrare anteriori ˆˆˆˆ ′′ I I I++bna,bn,1a++− N-bna,N-bn1,-aI 00 01 10 11 Biţii de intrare 00 00 01 10 11 ˆˆ′01 01 10 11 00 I′ I bna,bn,1a++− 10 10 11 00 01 11 11 00 01 10 Fig 4 17 Decodorul diferenţial pe 2 biţi (construcţia propusă) 125 -Teza de doctorat- Fig 4 16 Schema bloc a decodorului TCM 2N-D (construcţia propusă) 126 -Teza de doctorat- 4 7 Exemple pentru construcţia Wei generalizată a codurilor TCM 2N-D 4 7 1 Codul TCM 4-D, cu rata 2/3 şi cu 16 stări Am considerat transmisia a Q=7 biţi/interval de semnalizare 2-D Schema bloc a codorului TCM 4-D cu rata 2/3 şi 16 stări am prezentat-o în Fig 4 18 În acest caz N=2 este o putere întreagă a lui 2 deci am obţinut conform relaţiilor (4 124): 11  β ==  2N  0α =   1111 (4 136)  37500,1)1(2NCERCER ==++== 0  8N2N  24β)2(1 +  18182,2PAR ≅==  11CER 0  Deosebirea faţă de construcţia Wei a codului TCM 4-D, cu rata 2/3 şi cu 16 stări prezentată în paragraful 4 2, constă doar în calculul parametrilor CER şi PAR ai constelaţiei de semnale 2-D constituentă a constelaţiei de semnale 4-D Întrucît în acest caz N=2 este o putere întreagă a lui 2, am observat că valorile obţinute pentru parametrii CER şi PAR ai constelaţiei de semnale 2-D constituentă a constelaţiei de semnale 4-D sunt optimale deşi Wei nu a calculat aceste valori Constelaţia de semnale 2-D constituentă a constelaţiei de semnale 4-D este 2 un subset finit al latticei (2Z+1) (Z = mulţ imea numerelor întregi) şi este formată din 192 puncte de semnal 2-D: 7  2=128 puncte de semnal 2-D în grupul IG, 7 2 puncte de semnal 2-D în subgrupul OGL,  64= 2  0 puncte de semnal 2-D în subgrupul OGH u Am partiţionat constelaţia de semnale 2-D (Fig 4 11) constituentă în U=2=4 subseturi 2-D: 2 ⇒ A=(4Z+1), 2 ⇒ B=(4Z+3), ⇒ C=(4Z+1)(4Z+3), ⇒ D=(4Z+3)(4Z+1) 2 Dacă constelaţia de semnale 2-D constituentă are MSED=d atunci subseturile 0 22 2-D au MSED de intrasubset egală cu d4d= 0 127 -Teza de doctorat- Fig 4 18 Schema bloc a codorului TCM 4-D cu rata 2/3 şi 16 stări (construcţia Wei) 128 -Teza de doctorat- Numărul dedesubtul fiecărui punct de semnal 2-D este Z2,p, Z3,p, Z4,p, Z5,p, Z6,p, Z7,p, p=n, n+1 Fig 4 19 Constelaţia de semnale 2-D constituentă cu 192 puncte de semnal 2-D (construcţia Wei) 129 -Teza de doctorat- Dacă se rotesc succesiv subseturile 2-D în sens orar cu cîte 90° au loc următoarele transformări ale subseturilor 2-D (vezi relaţia 4 3) : A→C→B→D→A Prin concatenarea a 2 constelaţii de semnale 2-D constituente şi prin eliminarea punctelor de semnal 4-D care conţin mai mult decît un singur punct de semnal 2-D din 2x7+115 grupul extern OG se obţine constelaţia de semnale 4-D cu 2=2 puncte de semnal 4-D 22 Se construiesc 16 tipuri 4-D cu MSED=d4d= 0 prin concatenarea a cîte 2 subseturi 2-D şi notate astfel (vezi relaţia 4 4) : (A, A), (A, B), (A, C), (A, D),…,(D, D) Cele 16 tipuri 4-D sunt grupate în 8 subseturi 4-D (deci S=2) astfel încît MSED 22 a fiecărui subset 4-D obţinut să fie egală cu d4d= ia 4 5) : 0 (vezi relaţ S0=(A, A) ∪ (B, B) (D, D) S1=(C, C) ∪ (B, A) S2=(A, B) ∪ (D, C) S3=(C, D) ∪ (B, D) S4=(A, C) ∪ (D, A) S5=(C, B) ∪ (B, C) S6=(A, D) ∪ (D, B) S7=(C, A) ∪ Dacă se rotesc succesiv subseturile 2-D în sens orar cu cîte 90° au loc următoarele transformări ale subseturilor 4-D (vezi relaţia 4 6) : S1→S0 S0→ S3→S2 S2→ S5→S4 S4→ S7→S6 S6→ 2 Subseturile 4-D sunt grupate în 2 familii 4-D cu MSED=d2 astfel încît fiecare 0 familie 4-D să conţină F=4 subseturi 4-D (vezi relaţia 4 7) : F0=S0 ∪ S1 ∪ S2 ∪ S3 S5 ∪ S6 ∪ S7 F1=S4 ∪ În acest caz avem m=2 biţi de la intrare care sunt codaţi cu un codor 4 convoluţional cu rata de codare 2/3 şi cu V=2=16 stă ri De asemenea, restul şi cîtul împarţirii lui Nu=4 la Q=7 sunt r=4 şi respectiv c=0 130 -Teza de doctorat- Setul stărilor curente {} a codorului convoluţional se partiţionează în F=4 nσ subseturi (vezi relaţia 4 8) : j4i≤≤≤≤+ ,3i0 3j0 {} a codorului convoluţional se partiţionează în Setul stărilor următoare {} 2nσ + F=4 subseturi (vezi relaţia 4 9) : k4j≤≤≤≤+ ,3j0 3k0 {} Conectivitatea diagramei trellis se realizează astfel: pentru orice 0≤j≤3 există tranziţii de stare între fiecare stare curentă 4i+j, 0≤i≤3 şi fiecare stare următoare 4j+k, 0≤k≤3 Asignarea subseturilor 4-D la tranziţiile de stare ale codului TCM 4-D este realizată astfel încît să fie îndeplinite trei cerinţe : 1 Subseturile 4-D, asociate cu tranziţiile care pornesc dintr-o anumită stare, să difere unele de altele şi să aparţină aceleaşi familii 4-D Acelaşi considerent se aplică şi pentru subseturile 4-D asociate cu tranziţiile care ajung într-o anumită stare Conform acestei cerinţe: ⇒ la tranziţiile care pornesc din stări pare (Y0,n=W4,n=0) sunt asignate subseturi 4-D din familia 4-D notată F0, ⇒ la tranziţiile care pornesc din stări impare (Y0,n=W4,n=1) sunt asignate subseturi 4-D din familia 4-D notată F1 2 MSED dintre două secvenţe de subseturi 4-D, care corespund la două traiectorii 222 distincte prin diagrama trellis, să fie mai mare sau egală cu 4ddd= , care este 0= free MSED a fiecărui subset 4-D (MSED de intersubset 4-D să fie mai mare sau cel puţin egală cu MSED de intrasubset 4-D) Această cerinţă garantează că MSED dintre două secvenţe de puncte de semnal 22 i elimină evenimentele eroare care diferă cu mai mult decît un 4-D este d4d= 0 ş punct de semnal 4-D faţă de secvenţa corectă de puncte de semnal 4-D Dacă se neglijează efectul de margine, atunci coeficientul de eroare Nfree al codului TCM 4-D este egal cu numărul vecinilor cei mai apropiaţi ai unui punct de semnal 4-D şi aflaţi într-un subset 4-D cu un număr finit de puncte de semnal 4-D În acest caz Nfree=24 per punct de semnal 4-D (Nfree=12 per punct de semnal 2-D) Dacă nu se neglijează efectul de margine, atunci coeficientul de eroare al codului TCM 4-D este mai mic decît 24 3 Dacă X este subsetul 4-D asociat cu tranziţiile din starea curentă σ n în starea şi dacă X1 este subsetul 4-D obţinut prin rotirea lui X cu 90° atunci se următoare σ 2n+ poate defini o funcţie F1 de corespondenţă între stările codului TCM 4-D astfel încît X1 să fie asociat cu tranziţiile din starea curentă )(σ F în starea urmă)(σF n1toare 2n1+ 131 -Teza de doctorat- Pentru o rotaţie de 90° în sens orar se defineşte funcţia F1 astfel (vezi relaţia 4 10) : n4,n3,n2,n1,n4,n3,n2,n1,WWWWWWWW→ ii de unde bara deasupra înseamnă inversare binară şi W1,n, W2,n, W3,n, W4,n sunt biţ stare curentă a codorului convoluţional Acestă cerinţă garantează că vom obţine un cod TCM 4-D transparent la toate ambiguităţile de fază ale constelaţiei de semnale 4-D, cu alte cuvinte un cod TCM 4-D invariant rotaţional la o rotaţie de 0°, 90°, 180°, 270° Pentru a suprima toate cele 4 ambiguităţi de fază a constelaţiei de semnale 4-D (0°, 90°, 180°, 270°) se foloseşte un codor diferenţial pe 2 biţi (vezi relaţia 4 11) : ′′′′ )mod4IIII()II(+ n3,n2,2-n3,2-n2,n3,n2,= deci în acest caz avem a=3 şi b=0 Convertorul de bit converteşte cei 4 biţi de la intrarea sa Y0,n, I1,n, I′ 2,n, I′3,n în 2 grupuri de cîte 2 biţi fiecare (vezi relaţia 4 12) : Z0,p, Z1,p, p=n, n+1 care selectează 2 subseturi 2-D, corespunzătoare tipului 4-D selectat în cadrul subsetului 4-D utilizat Pentru a obţine un cod TCM 4-D transparent la orice ambiguitate de fază a constelaţiei de semnale 4-D trebuie sa fie îndeplinită cerinţa prezentată în continuare ′ Dacă X este tipul 4-D asociat cu structura de biţi ,Y′ ,I ,I I n3,n2,n1,n,0 atunci se notează cu X1, X2, X3 tipurile 4-D obţ inute prin rotirea lui X cu 90°, 180°, 270° ′′′′′ i obţinute Dacă ,I(′ )I 1n3,n2,, ,I(, 2n,3n2,)I ,I( 3n3,n2,)I sunt perechile de biţ ′ prin translaţia perechii de biţi ,I(′ )I sau trei poziţii în secvenţa n3,n2, cu una, două circulară 00→10→01→11→00, atunci tipurile 4-D asociate cu structurile de biţi: ′′ ,Y, ,I(I 1n3,n2,n1,n,0)I ′′ ,Y, ,I(I 2n3,n2,n1,n,0)I ′′ ,Y, ,I(I 3n3,n2,n1,n,0)I sunt X1, X2 ,X3 Corespondenţa dintre cei 2 biţi Z0,p, Z1,p, p=n, n+1 ş i subseturile 2-D este prezentată în Tabelul 4 1 iar funcţionarea convertorului de bit este prezentată în Tabelul 4 2 Tabelul 4 1 Corespondenţa dintre cei 2 biţi Z0,p, Z1,p ş i subseturile 2-D Subset 2-D Z0,p Z1,p A 0 0 B 0 1 C 1 0 D 1 1 132 -Teza de doctorat- Tabelul 4 2 Convertorul de bit Familie 4-D Subset 4-D Y0,n I1,n I′ 2,n I′3,n Tip 4-D Z0,n Z1,n Z0,n+1 Z1,n+1 0 0 0 0 0 0 (A, A) 0 0 0 0 0 0 0 1 (B, B) 0 1 0 1 1 0 0 1 0 (C, C) 1 0 1 0 0 0 1 1 (D, D) 1 1 1 1 2 0 1 0 0 (A, B) 0 0 0 1 0 1 0 1 (B, A) 0 1 0 0 3 0 1 1 0 (C, D) 1 0 1 1 0 1 1 1 (D, C) 1 1 1 0 1 4 1 0 0 0 (A, C) 0 0 1 0 1 0 0 1 (B, D) 0 1 1 1 5 1 0 1 0 (C, B) 1 0 0 1 1 0 1 1 (D, A) 1 1 0 0 6 1 1 0 0 (A, D) 0 0 1 1 1 1 0 1 (B, C) 0 1 1 0 7 1 1 1 0 (C, A) 1 0 0 0 1 1 1 1 (D, B) 1 1 0 1 Am notat cu Si subseturile 4-D ş i astfel indexul i este furnizat de relaţia: ′ += , i = 0, 1,…,7 (4 137) n2,n1,n0,I2IY4i+ Codorul bloc foloseşte cei 3 biţi I4,n, I5,n, I6,n pentru a genera 2 grupuri de cîte 2 biţi fiecare (vezi relaţia 4 13) : Z2,p, Z3,p, p=n, n+1 care selectează pentru cele 2 subseturi 2-D utilizate, prima jumătate a grupului intern IG1, a doua jumă tate a grupului intern IG2 sau subgrupul extern de energie mică OGL Corespondenţa dintre cei 2 biţi Z2,p, Z3,p, p=n, n+1 ş i grupurile constelaţiei de semnale 2-D este prezentată în Tabelul 4 16 iar funcţionarea codorului bloc este prezentată în Tabelul 4 4 Cele 2 grupuri de cîte 4 biţi fiecare (vezi relaţia 4 14) : Z4,p, Z5,p, Z6,p, Z7,p, p=n, n+1 selectează pentru fiecare din cele 2 subseturi 2-D utilizate, cîte un punct de semnal 2-D din grupul constelaţiei de semnale 2-D selectat anterior 133 -Teza de doctorat- Tabelul 4 16 Corespondenţa dintre cei 2 biţi Z2,p, Z3,p ş i grupurile constelaţiei de semnale 2-D Z2,p Z3,p Grup constelaţ ie de semnale 2-D 0 0 IG1 0 1 IG2 1 0 OGL Tabelul 4 17 Codorul bloc I4,n I5,n I6,n Z2,n Z3,n Z2,n+1 Z3,n+1 0 0 0 0 0 0 0 0 0 1 0 0 0 1 0 1 0 0 0 1 0 0 1 1 0 1 1 0 1 0 0 1 0 0 0 1 0 1 1 0 0 1 1 1 0 0 1 0 0 1 1 1 0 1 0 1 Cîştigul de codare asimptotic al codului TCM 4-D faţă de o transmisie necodată este conform relaţiei (4 134):     321    = [dB] (4 138) 63757,4log104log10γ≅=  1010 1111      8  4 7 2 Codul TCM 6-D, cu rata 3/4 şi cu 64 stări Am considerat transmisia a Q=7 biţi/interval de semnalizare 2-D Schema bloc a codorului TCM 6-D cu rata 3/4 şi 64 stări am prezentat-o în Fig 4 20 În acest caz N=3 nu este o putere întreagă a lui 2 deci am obţinut conform relaţiilor (4 125): 134 -Teza de doctorat- ′  11)I(4K 0 β ==  32Q 2   14 α ==  ′ 11 )I(4K  0  ′ 1KP ==  0  2P-N =  ′ Q   4K2 + ′ 0  44)I(4K == 0 N   1111 22222,1)1(2NCER ≅=++=  0 9N2N   )(1 11α − ][β-[CERCER +=1,22233]-β ≅ 0  )N(12 Nα −   β)2(1 + 19867,2PAR ≅=  11α)(1 (4 139) − ]-β][β-[CER +  0  Nα)N(12 −  Constelaţia de semnale 2-D constituentă a constelaţiei de semnale 6-D este 2 imea numerelor întregi) şi este formată din un subset finit al latticei (2Z+1) (Z = mulţ 172 puncte de semnal 2-D: 7  2=128 puncte de semnal 2-D în grupul IG,  40 puncte de semnal 2-D în subgrupul OGL,  4 puncte de semnal 2-D în subgrupul OGH u Am partiţionat constelaţia de semnale 2-D (Fig 4 21) constituentă în U=2=4 subseturi 2-D: 2 ⇒ A=(4Z+1), 2 ⇒ B=(4Z+3), ⇒ C=(4Z+1)(4Z+3), ⇒ D=(4Z+3)(4Z+1) 2 Dacă constelaţia de semnale 2-D constituentă are MSED=d atunci subseturile 0 22 2-D au MSED de intrasubset egală cu d4d= 0 Dacă se rotesc succesiv subseturile 2-D în sens orar cu cîte 90° au loc următoarele transformări ale subseturilor 2-D (vezi relaţia 4 3) : A→C→B→D→A Prin concatenarea a 3 constelaţii de semnale 2-D constituente şi prin eliminarea punctelor de semnal 6-D care conţin mai mult decît un singur punct de semnal 2-D din 3x7+122 grupul extern OG se obţine constelaţia de semnale 6-D cu 2=2 puncte de semnal 6-D 135 -Teza de doctorat- Fig 4 20 Schema bloc a codorului TCM 6-D cu rata 2/3 şi 16 stări (construcţia propusă) 136 -Teza de doctorat- Numărul dedesubtul fiecărui punct de semnal 2-D este Z2,p, Z3,p, Z4,p, Z5,p, Z6,p, Z7,p, p=n, n+1, n+2 Fig 4 21 Constelaţia de semnale 2-D constituentă cu 172 puncte de semnal 2-D (construcţia propusă) 137 -Teza de doctorat- 22 Se construiesc 64 tipuri 6-D cu MSED=d4d= prin concatenarea a cîte 3 0 subseturi 2-D şi notate astfel (vezi relaţia 4 34) : (A, A, A), (A, A, B), (A, A, C), (A, A, D),…,(D, D, D) Cele 64 tipuri 6-D sunt grupate în 16 subseturi 6-D (deci S=4) astfel încît 22 MSED a fiecărui subset 6-D obţinut să fie egală cu d4d= ia 4 35) : 0 (vezi relaţ S0=(A, A, A) ∪ (A, B, B) ∪ (B, A, B) ∪ (B, B, A) (C, B, A) ∪ (D, A, A) ∪ (D, B, B) S1=(C, A, B) ∪ (A, D, D) ∪ (B, C, D) ∪ (B, D, C) S2=(A, C, C) ∪ (C, D, D) ∪ (D, C, D) ∪ (D, D, C) S3=(C, C, C) ∪ (A, B, A) ∪ (B, A, A) ∪ (B, B, B) S4=(A, A, B) ∪ (C, B, B) ∪ (D, A, B) ∪ (D, B, A) S5=(C, A, A) ∪ (A, D, C) ∪ (B, C, C) ∪ (B, D, D) S6=(A, C, D) ∪ (C, D, C) ∪ (D, C, C) ∪ (D, D, D) S7=(C, C, D) ∪ (C, B, D) ∪ (D, A, D) ∪ (D, B, C) S8=(C, A, C) ∪ (A, B, C) ∪ (B, A, C) ∪ (B, B, D) S9=(A, A, D) ∪ (C, D, B) ∪ (D, C, B) ∪ (D, D, A) S10=(C, C, A) ∪ (A, D, B) ∪ (B, C, B) ∪ (B, D, A) S11=(A, C, A) ∪ (C, B, C) ∪ (D, A, C) ∪ (D, B, D) S12=(C, A, D) ∪ (A, B, D) ∪ (B, A, D) ∪ (B, B, C) S13=(A, A, C) ∪ (C, D, A) ∪ (D, C, A) ∪ (D, D, B) S14=(C, C, B) ∪ (A, D, A) ∪ (B, C, A) ∪ (B, D, B) S15=(A, C, B) ∪ Dacă se rotesc succesiv subseturile 2-D în sens orar cu cîte 90° au loc următoarele transformări ale subseturilor 6-D (vezi relaţia 4 36) : S3→S4→S7→S0 S0→ S2→S5→S6→S1 S1→ S8→S11→S12→S15→S8 S10→S13→S14→S9 S9→ 2 Subseturile 6-D sunt grupate în 2 familii 6-D cu MSED=d2 astfel încît fiecare 0 familie 6-D să conţină F=8 subseturi 6-D (vezi relaţia 4 37) : F0=S0 ∪ S2 ∪ S4 ∪ S6 ∪ S8 ∪ S10 ∪ S12 ∪ S14 S3 ∪ S5 ∪ S7 ∪ S9 ∪ S11 ∪ S13 ∪ S15 F1=S1 ∪ În acest caz avem m=3 biţi de la intrare care sunt codaţi cu un codor 6 convoluţional cu rata de codare 3/4 şi cu V=2=64 stă ri De asemenea, restul şi cîtul împărţirii lui Nu=6 la Q=7 sunt r=6 şi respectiv c=0 138 -Teza de doctorat- Setul stărilor curente {} a codorului convoluţional se partiţionează în F=8 nσ subseturi (vezi relaţia 4 38) : j8i≤≤≤≤+ ,7i0 7j0 {} a codorului convoluţional se partiţionează în Setul stărilor următoare {} 3nσ + F=8 subseturi (vezi relaţia 4 39) : k8j≤≤≤≤+ ,7j0 7k0 {} Conectivitatea diagramei trellis se realizează astfel: pentru orice 0≤j≤7 există tranziţii de stare între fiecare stare curentă 8i+j, 0≤i≤7 şi fiecare stare următoare 8j+k, 0≤k≤7 Asignarea subseturilor 6-D la tranziţiile de stare ale codului TCM 6-D este realizată astfel încît să fie îndeplinite trei cerinţe : 1 Subseturile 6-D, asociate cu tranziţiile care pornesc dintr-o anumită stare, să difere unele de altele şi să aparţină aceleaşi familii 6-D Acelaşi considerent se aplică şi pentru subseturile 6-D asociate cu tranziţiile care ajung într-o anumită stare Conform acestei cerinţe: ⇒ la tranziţiile care pornesc din stări pare (Y0,n=W6,n=0) sunt asignate subseturi 6-D din familia 6-D notată F0, ⇒ la tranziţiile care pornesc din stări impare (Y0,n=W6,n=1) sunt asignate subseturi 6-D din familia 6-D notată F1 2 MSED dintre două secvenţe de subseturi 6-D, care corespund la două traiectorii 222 distincte prin diagrama trellis, să fie mai mare sau egală cu 4ddd= , care este 0= free MSED a fiecărui subset 6-D (MSED de intersubset 6-D să fie mai mare sau cel puţin egală cu MSED de intrasubset 6-D) Această cerinţă garantează că MSED dintre două secvenţe de puncte de semnal 22 i elimină evenimentele eroare care diferă cu mai mult decît un 6-D este d4d= 0 ş punct de semnal 6-D faţă de secvenţa corectă de puncte de semnal 6-D Dacă se neglijează efectul de margine, atunci coeficientul de eroare Nfree al codului TCM 6-D este egal cu numărul vecinilor cei mai apropiaţi ai unui punct de semnal 6-D şi aflaţi într-un subset 6-D cu un număr finit de puncte de semnal 6-D În acest caz Nfree=60 per punct de semnal 6-D (Nfree=20 per punct de semnal 2-D) Dacă nu se neglijează efectul de margine, atunci coeficientul de eroare al codului TCM 6-D este mai mic decît 60 3 Dacă X este subsetul 6-D asociat cu tranziţiile din starea curentă σ n în starea următoare σ şi dacă X1, X2, X3 sunt subseturile 6-D obţinute prin rotirea lui X cu 3n+ între stările 90°, 180°, 270° atunci se pot defini trei funcţii F1, F2, F3 de corespondenţă fie asociat cu tranziţiile din starea curentă )(σ codului TCM 6-D astfel încît X1 să F n1 în starea următoare )(σ fie asociat cu tranziţiile din starea curentă )(σ F, X2 săF 3n1+n2 în starea următoare )(σ fie asociat cu tranziţiile din starea curentă )(σ F, X3 săF 3n2+n3 139 -Teza de doctorat- în starea următoare )(σ F 3n3+ ia Pentru o rotaţie de 90° în sens orar se defineşte funcţia F1 astfel (vezi relaţ 4 40) : n6,n5,n4,n3,n2,n1,n6,n5,n4,n3,n2,n1,VVWVVWWWWWWW→ ii de stare curentă a codorului convoluţional şi (vezi unde W1,n, W2,n ,…, W6,n sunt biţ relaţia 4 41) : 01)mod4W(WVV+ n3,n2,n3,n2,= 4mod)01W(WVV+ n6,n5,n6,n5,= ia Pentru o rotaţie de 180° în sens orar se defineşte funcţia F2 astfel (vezi relaţ 4 42) : n6,n5,n4,n3,n2,n1,n6,n5,n4,n3,n2,n1,WWWWWWWWWWWW→ unde bara deasupra înseamnă inversare binară ia 4 43) Pentru o rotaţie de 270° în sens orar definim funcţia F3 astfel (vezi relaţ : n6,n5,n4,n3,n2,n1,n6,n5,n4,n3,n2,n1,VVWVVWWWWWWW→ unde (vezi relaţia 4 44) : 1)mod41W(WVV+ n3,n2,n3,n2,= 4mod)11W(WVV+ n6,n5,n6,n5,= Acestă cerinţă garantează că vom obţine un cod TCM 6-D transparent la toate ambiguităţile de fază ale constelaţiei de semnale 6-D, cu alte cuvinte un cod TCM 6-D invariant rotaţional la o rotaţie de 0°, 90°, 180°, 270° Pentru a suprima toate cele 4 ambiguităţi de fază a constelaţiei de semnale 6-D (0°, 90°, 180°, 270°) se foloseşte un codor diferenţial pe 2 biţi (vezi relaţia 4 45) : ′′′′ )mod4IIII()II(+ n2,n1,3-n2,3-n1,n2,n1,= deci în acest caz avem a=2 şi b=0 Convertorul de bit converteşte cei 6 biţi de la intrarea sa Y0,n, I′ 1,n, I′2,n, I3,n, I4,n, I5,n în 3 grupuri de cîte 2 biţ i fiecare (vezi relaţia 4 46) : Z0,p, Z1,p, p=n, n+1, n+2 care selectează 3 subseturi 2-D, corespunzătoare tipului 6-D selectat în cadrul 140 -Teza de doctorat- subsetului 6-D utilizat Pentru a obţine un cod TCM 6-D transparent la orice ambiguitate de fază a constelaţiei de semnale 6-D trebuie sa fie îndeplinită cerinţa prezentată în continuare ′ Dacă X este tipul 6-D asociat cu structura de biţi ,Y′ ,I ,I II,I n5,n4,n3,n2,n1,n,0 atunci se notează cu X1, X2, X3 tipurile 6-D obţ inute prin rotirea lui X cu 90°, 180°, 270° ′′′′′ i obţinute prin Dacă ,I(′ )I 1n2,n1,, ,I(, 2n,2n1,)I ,I( 3n2,n1,)I sunt perechile de biţ ′ translaţia perechii de biţi )I,I(′ sau trei poziţii în secvenţa circulară n2,n1, cu una, două 00→10→01→11→00, atunci tipurile 6-D asociate cu structurile de biţi: ′′ ,Y, ,)I,I( n5,n4,n3,1n2,n1,n,0II,I ′′ ,Y, ,)I,I( n5,n4,n3,2n2,n1,n,0II,I ′′ ,Y, ,)I,I( n5,n4,n3,3n2,n1,n,0II,I sunt X1, X2 ,X3 Corespondenţa dintre cei 2 biţi Z0,p, Z1,p, p=n, n+1, n+2 ş i subseturile 2-D este prezentată în Tabelul 4 1 iar funcţionarea convertorului de bit este prezentată în Tabelul 4 8 Am notat cu Si subseturile 6-D ş i astfel indexul i este furnizat de relaţia (vezi relaţia 4 47) : ′′ ++=, i = 0, 1,…,15 n3,n2,n1,n0,II2I48Yi+ Codorul bloc foloseşte cei 16 biţi I6,, I7,n, I1,n+1, I2,n+1, I3,n, I4,n+1 I5,n, I6,n+1, n+1+1 I7,n+1, I1, I2,n+2, I3,, I4,n+2 I5,n+2, I6,n+2, I7,n+2 pentru a genera 3 grupuri de cîte 6 biţ ,n+2n+2i fiecare: Z2,p, Z3,p, Z4,p, Z5,p, Z6,p, Z7,p, p=n, n+1, n+2 (4 140) care selectează pentru fiecare din cele 3 subseturi 2-D utilizate, punctul de semnal 2-D care va fi transmis Funcţionarea codorului bloc este prezentată în Tabelul 4 18 Cîştigul de codare asimptotic al codului TCM 6-D faţă de o transmisie necodată este conform relaţiei (4 134): 4  14871,5log10γ≅ = [dB] (4 141)  10 222,1   Tabelul 4 8 Convertorul de bit Subset 6-D Y0,n I′ I4,n I5,n Tip 6-D Z0,n Z1,n Z0,n+1 Z1,n+1 Z0,n+2 Z1,n+2 1,n I′2,n I3,n 0 0 0 0 0 0 0 (A, A, A) 0 0 0 0 0 0 0 0 0 0 0 1 (A, B, B) 0 0 0 1 0 1 0 0 0 0 1 0 (B, A, B) 0 1 0 0 0 1 0 0 0 0 1 1 (B, B, A) 0 1 0 1 0 0 … 15 1 1 1 1 0 0 (A, C, B) 0 0 1 0 0 1 1 1 1 1 0 1 (A, D, A) 0 0 1 1 0 0 1 1 1 1 1 0 (B, C, A) 0 1 1 1 0 0 1 1 1 1 1 1 (B, D, B) 0 1 1 1 0 1 141 -Teza de doctorat- Tabelul 4 18 Codorul bloc 7,n+2 Z7,n+27,n+27,n+27,n+27,n+24,n+14,n+17,n+27,n+2 I IIIIIIII 6,n+2 Z6,n+26,n+26,n+26,n+26,n+23,n+16,n+26,n+2 I IIIII0 I IZ5,n+25,n+25,n+25,n+25,n+25,n+22,n+15,n+25,n+2 I IIIII0 I IZ4,n+24,n+24,n+24,n+24,n+24,n+24,n+24,n+2 I IIII0 I Z3,n+23,n+23,n+23,n+23,n+23,n+23,n+23,n+21 I IIII0 I I Z2,n+20 0 I 0 7,n+12,n+22,n+22,n+24,n+14,n+17,n+27,n+22,n+2 Z I IIIIIII1 0 6,n+1 Z1,n+21,n+21,n+23,n+16,n+26,n+21,n+2 I 0 III0 II1 5,n+1 7,n+17,n+17,n+12,n+15,n+25,n+27,n+1 0 IZ I III0 II1 1 4,n+1 IZ6,n+16,n+16,n+14,n+24,n+26,n+1 I II0 II1 1 Z3,n+15,n+15,n+15,n+13,n+23,n+25,n+11 III00 II1 I 2,n+1 0 Z0 0 I 0 7,n4,n+14,n+14,n+12,n+22,n+22,n+22,n+22,n+2 Z IIIIIII1 0 IZ6,n3,n+13,n+11,n+21,n+21,n+21,n+21,n+2 1 II0 III1 5,n 2,n+12,n+17,n+17,n+17,n+17,n+17,n+1 1 IIZ II0 III1 0 4,n IIZ1,n+16,n+16,n+16,n+16,n+16,n+1 I0 III1 0 1 IZ3,n7,n5,n+15,n+15,n+15,n+15,n+1 0 III1 I 0 I 2,n 1 IZ 0 0 I4,n+1 1 I X 1 X I3,n+1 X 0 X X I2,n+1 X 0 X 0 X 0 I1,n+1 X X X 0 0 X 7,n 0 I 1 X 1 X X 1 X 6,n I 0 0 0 0 1 0 0 X 1 1 1 1 1 1 0 0 1 142 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -Teza de doctorat- 4 7 3 Codul TCM 8-D, cu rata 3/4 şi cu 64 stări Am considerat transmisia a Q=7 biţi/interval de semnalizare 2-D Schema bloc a codorului TCM 8-D cu rata 3/4 şi 64 stări am prezentat-o în Fig 4 22 În acest caz N=4 este o putere întreagă a lui 2 deci am obţinut conform relaţiilor (4 124): 11  β ==  4N  0α =   3711 (4 142)  15625,1)1(2NCERCER ≅=++== 0  32N2N  80β)2(1 +  16216,2PAR ≅==  37CER 0  Deosebirea faţă de construcţia Wei a codului TCM 8-D, cu rata 3/4 şi cu 64 stări prezentată în paragraful 4 2, constă doar în calculul parametrilor CER şi PAR ai constelaţiei de semnale 2-D constituentă a constelaţiei de semnale 8-D Întrucît în acest caz N=2 este o putere întreagă a lui 2, am observat că valorile obţinute pentru parametrii CER şi PAR ai constelaţiei de semnale 2-D constituentă a constelaţiei de semnale 8-D sunt optimale deşi Wei nu a calculat aceste valori Constelaţia de semnale 2-D constituentă a constelaţiei de semnale 8-D este 2 un subset finit al latticei (2Z+1) (Z = mulţ imea numerelor întregi) şi este formată din 160 puncte de semnal 2-D: 7  2=128 puncte de semnal 2-D în grupul IG, 7 2 puncte de semnal 2-D în subgrupul OGL,  32= 4  0 puncte de semnal 2-D în subgrupul OGH u Am partiţionat constelaţia de semnale 2-D (Fig 4 23) constituentă în U=2=4 subseturi 2-D: 2 ⇒ A=(4Z+1), 2 ⇒ B=(4Z+3), ⇒ C=(4Z+1)(4Z+3), ⇒ D=(4Z+3)(4Z+1) 2 Dacă constelaţia de semnale 2-D constituentă are MSED=d atunci subseturile 0 22 2-D au MSED de intrasubset egală cu d4d= 0 Dacă se rotesc succesiv subseturile 2-D în sens orar cu cîte 90° au loc următoarele transformări ale subseturilor 2-D (vezi relaţia 4 3) : A→C→B→D→A (4 3) 143 -Teza de doctorat- Fig 4 22 Schema bloc a codorului TCM 8-D cu rata 3/4 şi 64 stări (construcţia Wei) 144 -Teza de doctorat- Numărul dedesubtul fiecărui punct de semnal 2-D este Z2,p, Z3,p, Z4,p, Z5,p, Z6,p, Z7,p, p=n, n+1, n+2, n+3 Fig 4 23 Constelaţia de semnale 2-D constituentă cu 160 puncte de semnal 2-D (construcţia Wei) 145 -Teza de doctorat- Prin concatenarea a 4 constelaţii de semnale 2-D constituente şi prin eliminarea punctelor de semnal 8-D care conţin mai mult decît un singur punct de semnal 2-D din 4x7+129 grupul extern OG se obţine constelaţia de semnale 8-D cu 2=2 puncte de semnal 8-D 22 Se construiesc 64 tipuri 8-D cu MSED=d4d= 0 prin concatenarea a cîte 2 subseturi 4-D şi notate astfel (vezi relaţia 4 16) : (0, 0), (0, 1),…,(7, 7) unde i=0, 1,…,7 reprezintă indicele subseturilor 4-D notate Si în paragraful 4 7 2 Cele 64 tipuri 8-D sunt grupate în 16 subseturi 8-D (deci S=4) astfel încît 22 MSED a fiecărui subset 8-D obţinut să fie egală cu d4d= ia 4 17) : 0 (vezi relaţ S0=(0, 0) ∪ (1, 1) ∪ (2, 2) ∪ (3, 3) (1, 0) ∪ (2, 3) ∪ (3, 2) S1=(0, 1) ∪ (1, 3) ∪ (2, 0) ∪ (3, 1) S2=(0, 2) ∪ (1, 2) ∪ (2, 1) ∪ (3, 0) S3=(0, 3) ∪ (5, 5) ∪ (6, 6) ∪ (7, 7) S4=(4, 4) ∪ (5, 4) ∪ (6, 7) ∪ (7, 6) S5=(4, 5) ∪ (5, 7) ∪ (6, 4) ∪ (7, 5) S6=(4, 6) ∪ (5, 6) ∪ (6, 5) ∪ (7, 4) S7=(4, 7) ∪ (1, 5) ∪ (2, 6) ∪ (3, 7) S8=(0, 4) ∪ (1, 4) ∪ (2, 7) ∪ (3, 6) S9=(0, 5) ∪ (1, 7) ∪ (2, 4) ∪ (3, 5) S10=(0, 6) ∪ (1, 6) ∪ (2, 5) ∪ (3, 4) S11=(0, 7) ∪ (5, 1) ∪ (6, 2) ∪ (7, 3) S12=(4, 0) ∪ (5, 0) ∪ (6, 3) ∪ (7, 2) S13=(4, 1) ∪ (5, 3) ∪ (6, 0) ∪ (7, 1) S14=(4, 2) ∪ (5, 2) ∪ (6, 1) ∪ (7, 0) S15=(4, 3) ∪ Dacă se rotesc succesiv subseturile 2-D în sens orar cu cîte 90° nu au loc transformări ale subseturilor 8-D, cu alte cuvinte subseturile 8-D sunt invariante la aceste rotaţii 2 Subseturile 8-D sunt grupate în 2 familii 8-D cu MSED=d2 astfel încît fiecare 0 familie 8-D să conţină F=8 subseturi 8-D (vezi relaţia 4 18) : F0=S0 ∪ S1 ∪ S2 ∪ S3 ∪ S4 ∪ S5 ∪ S6 ∪ S7 S9 ∪ S10 ∪ S11 ∪ S12 ∪ S13 ∪ S14 ∪ S1 F1=S8 ∪ În acest caz avem m=3 biţi de la intrare care sunt codaţi cu un codor 6 convoluţional cu rata de codare 3/4 şi cu V=2=64 stă ri De asemenea, restul şi cîtul împărţirii lui Nu=8 la Q=7 sunt r=1 şi respectiv c=1 146 -Teza de doctorat- Setul stărilor curente {} a codorului convoluţional se partiţionează în F=8 nσ subseturi (vezi relaţia 4 19) : j8i≤≤≤≤+ ,7i0 7j0 {} a codorului convoluţional se partiţionează în Setul stărilor următoare {} 4nσ + F=8 subseturi (vezi relaţia 4 20) : k8j≤≤≤≤+ ,7j0 7k0 {} Conectivitatea diagramei trellis se realizează astfel: pentru orice 0≤j≤7 există tranziţii de stare între fiecare stare curentă 8i+j, 0≤i≤7 şi fiecare stare următoare 8j+k, 0≤k≤7 Asignarea subseturilor 8-D la tranziţiile de stare ale codului TCM 8-D este realizată astfel încît să fie îndeplinite trei cerinţe : 1 Subseturile 8-D, asociate cu tranziţiile care pornesc dintr-o anumită stare, să difere unele de altele şi să aparţină aceleaşi familii 8-D Acelaşi considerent se aplică şi pentru subseturile 8-D asociate cu tranziţiile care ajung într-o anumită stare Conform acestei cerinţe: ⇒ la tranziţiile care pornesc din stări pare (Y0,n=W6,n=0) sunt asignate subseturi 8-D din familia 8-D notată F0, ⇒ la tranziţiile care pornesc din stări impare (Y0,n=W6,n=1) sunt asignate subseturi 8-D din familia 8-D notată F1 2 MSED dintre două secvenţe de subseturi 8-D, care corespund la două traiectorii 222 distincte prin diagrama trellis, să fie mai mare sau egală decît 4ddd= , care 0= free este MSED a fiecărui subset 8-D (MSED de intersubset 8-D să fie mai mare sau cel puţin egală cu MSED de intrasubset 8-D) Această cerinţă garantează că MSED dintre două secvenţe de puncte de semnal 22 i elimină evenimentele eroare care diferă cu mai mult decît un 8-D este d4d= 0 ş punct de semnal 8-D faţă de secvenţa corectă de puncte de semnal 8-D Dacă se neglijează efectul de margine, atunci coeficientul de eroare Nfree al codului TCM 8-D este egal cu numărul vecinilor cei mai apropiaţi ai unui punct de semnal 8-D şi aflaţi într-un subset 8-D cu un număr finit de puncte de semnal 8-D În acest caz Nfree=240 per punct de semnal 8-D (Nfree=60 per punct de semnal 2-D) Dacă nu se neglijează efectul de margine, atunci coeficientul de eroare al codului TCM 8-D este mai mic decît 240 3 Dacă X este subsetul 8-D asociat cu tranziţiile din starea curentă σ n în starea următoare σ şi dacă X1, X2, X3 sunt subseturile 8-D obţinute prin rotirea lui X cu 4n+ 90°, 180°, 270° atunci X=X1=X2=X3 Acestă cerinţă garantează că vom obţine un cod TCM 4-D transparent la toate ambiguităţile de fază ale constelaţiei de semnale 8-D, cu alte cuvinte un cod TCM 8-D invariant rotaţional la o rotaţie de 0°, 90°, 180°, 270° 147 -Teza de doctorat- Pentru a suprima toate cele 4 ambiguităţi de fază a constelaţiei de semnale 8-D (0°, 90°, 180°, 270°) se foloseşte un codor diferenţial pe 2 biţi (vezi relaţia 4 21) : ′′′′ )mod4IIII()II(+ n7,n6,4-n7,4-n6,n7,n6,= deci în acest caz avem a=7 şi b=0 Convertorul de bit converteşte cei 8 biţi de la intrarea sa Y0,n, I1,n, I2,, I3,n, I4,, nn i fiecare (vezi relaţia 4 22) : I5,n, I′ 6,n, I′7,n în 4 grupuri de cîte 2 biţ Z0,p, Z1,p, p=n, n+1, n+2, n+3 care selectează 4 subseturi 2-D, corespunzătoare tipului 8-D selectat în cadrul subsetului 8-D utilizat Corespondenţa dintre cei 2 biţi Z0,p, Z1,p, p=n, n+1, n+2, n+3 ş i subseturile 2-D este prezentată în Tabelul 4 1 iar funcţionarea convertorului de bit este prezentată în Tabelul 4 5 Am notat cu Si subseturile 8-D ş i astfel indexul i este furnizat de relaţia: ++=, i = 0, 1,…,15 (4 143) n3,n2,n1,n0,II2I4Y8i+ Codorul bloc foloseşte cei 9 biţi I1, I2,n+1, I3,n, I4,n+1, I5,n, I6,n+1, I7,n+1, I1,n+2, ,n+1+1+1 I2,n+2 pentru a genera 4 grupuri de cîte 3 biţ i fiecare (vezi relaţia 4 23) : Z2,p, Z3,p, Z4,p, p=n, n+1, n+2, n+3 care selectează pentru cele 4 subseturi 2-D utilizate primul sfert al grupului intern IG1, al doilea sfert al grupului intern IG2, al treilea sfert al grupului intern IG3, al patrulea sfert al grupului intern IG4 sau subgrupul extern de energie mică OGL Corespondenţa dintre cei 3 biţi Z2,p, Z3,p, Z4,p, p=n, n+1, n+2, n+3 ş i grupurile constelaţiei de semnale 2-D este prezentată în Tabelul 4 19 iar funcţionarea codorului bloc este prezentată în Tabelul 4 7 Cele 4 grupuri de cîte 3 biţi fiecare (vezi relaţia 4 24) : Z5,p, Z6,p, Z7,p, p=n, n+1, n+2, n+3 selectează pentru fiecare din cele 4 subseturi 2-D utilizate, cîte un punct de semnal 2-D din grupul constelaţiei de semnale 2-D selectat anterior Cîştigul de codare asimptotic al codului TCM 8-D faţă de o transmisie necodată este conform relaţiei (4 134):     1281    = [dB] (4 144) 39008,5log104log10γ≅=  1010 3737      32  148 -Teza de doctorat- Tabelul 4 1 Corespondenţa dintre cei 2 biţi Z0,p, Z1,p ş i subseturile 2-D Subset 2-D Z0,p Z1,p A 0 0 B 0 1 C 1 0 D 1 1 Tabelul 4 5 Convertorul de bit Subset 8-D Y0,n I1,n I2,n I3,n I4,n I5,n Tip 8-D I′ Z1,n Z0,n+ Z1,n+ Z0,n+ Z1,n+ Z0,n+ Z1,n+ 6,n I′7,n Subset 2-D Z0,n112233 0 0 0 0 0 0 0 (0, 0) 0 0 (A, A)∪0 0 0 0 0 0 0 0 (A, A) 0 0 0 0 0 0 0 1 (A, A)∪0 0 0 0 0 1 0 1 (B, B) 0 0 0 0 0 0 1 0 (B, B)∪0 1 0 1 0 0 0 0 (A, A) 0 0 0 0 0 0 1 1 (B, B)∪0 1 0 1 0 1 0 1 (B, B) 0 0 0 0 0 1 (1, 1) 0 0 (C, C)∪1 0 1 0 1 0 1 0 (C, C) 0 0 0 0 0 1 0 1 (C, C)∪1 0 1 0 1 1 1 1 (D, D) 0 0 0 0 0 1 1 0 (D, D)∪1 1 1 1 1 0 1 0 (C, C) 0 0 0 0 0 1 1 1 (D, D)∪1 1 1 1 1 1 1 1 (D, D) 0 0 0 0 1 0 (2, 2) 0 0 (A, B)∪0 0 0 1 0 0 0 1 (A, B) 0 0 0 0 1 0 0 1 (A, B)∪0 0 0 1 0 1 0 0 (B, A) 0 0 0 0 1 0 1 0 (B, A)∪0 1 0 0 0 0 0 1 (A, B) 0 0 0 0 1 0 1 1 (B, A)∪0 1 0 0 0 1 0 0 (B, A) 0 0 0 0 1 1 (3, 3) 0 0 (C, D)∪1 0 1 1 1 0 1 1 (C, D) 0 0 0 0 1 1 0 1 (C, D)∪1 0 1 1 1 1 1 0 (D, C) 0 0 0 0 1 1 1 0 (D, C)∪1 1 1 0 1 0 1 1 (C, D) 0 0 0 0 1 1 1 1 (D, C)∪ (D, C) 1 1 1 0 1 1 1 0 … 15 1 1 1 1 0 0 (4, 3) 0 0 (A, C)∪0 0 1 0 1 0 1 1 (C, D) 1 1 1 1 0 0 0 1 (A, C)∪0 0 1 0 1 1 1 0 (D, C) 1 1 1 1 0 0 1 0 (B, D)∪0 1 1 1 1 0 1 1 (C, D) 1 1 1 1 0 0 1 1 (B, D)∪0 1 1 1 1 1 1 0 (D, C) 1 1 1 1 0 1 (5, 2) 0 0 (C, B)∪1 0 0 1 0 0 0 1 (A, B) 1 1 1 1 0 1 0 1 (C, B)∪1 0 0 1 0 1 0 0 (B, A) 1 1 1 1 0 1 1 0 (D, A)∪1 1 1 1 0 0 0 1 (A, B) 1 1 1 1 0 1 1 1 (D, A)∪1 1 1 1 0 1 0 0 (B, A) 1 1 1 1 1 0 (6, 1) 0 0 (A, D)∪0 0 1 1 1 0 1 0 (C, C) 1 1 1 1 1 0 0 1 (A, D)∪0 0 1 1 1 1 1 1 (D, D) 1 1 1 1 1 0 1 0 (B, C)∪0 1 1 0 1 0 1 0 (C, C) 1 1 1 1 1 0 1 1 (B, C)∪0 1 1 0 1 1 1 1 (D, D) 1 1 1 1 1 1 (7, 0) 0 0 (C, A)∪1 0 0 0 0 0 0 0 (A, A) 1 1 1 1 1 1 0 1 (C, A)∪1 0 0 0 0 1 0 1 (B, B) 1 1 1 1 1 1 1 0 (D, B)∪1 1 0 1 0 0 0 0 (A, A) 1 1 1 1 1 1 1 1 (D, B)∪1 1 0 1 0 1 0 1 (B, B) Tabelul 4 19 Corespondenţa dintre cei 3 biţi Z2,p, Z3,p Z4,p ş i grupurile constelaţiei de semnale 2-D Z2,p Z3,p Z4,p Grup constelaţ ie de semnale 2-D 0 0 0 IG1 0 0 1 IG2 0 1 0 IG3 0 1 1 IG4 1 0 0 OGL Tabelul 4 7 Codorul bloc I1,n+1 I2,n+1 I3,n+1 Z2,n Z3,n Z4,n Z2,n+1 Z3,n+1 Z4,n+1 Z2,n+2 Z3,n+2 Z4,n+2 Z2,n+3 Z3,n+3 Z4,n+3 0 X X 0 I2,n+1 I3,n+1 0 I4,n+1 I5,n+1 0 I6,n+1 I7,n+1 0 I1,n+2 I2,n+2 1 0 0 1 0 0 0 I4,n+1 I5,n+1 0 I6,n+1 I7,n+1 0 I1,n+2 I2,n+2 1 0 1 0 I4,n+1 I5,n+1 1 0 0 0 I6,n+1 I7,n+1 0 I1,n+2 I2,n+2 1 1 0 0 I4,n+1 I5,n+1 0 I6,n+1 I7,n+1 1 0 0 0 I1,n+2 I2,n+2 1 1 1 0 I4,n+1 I5,n+1 0 I6,n+1 I7,n+1 0 I1,n+2 I2,n+2 1 0 0 149 -Teza de doctorat- 4 7 4 Codul TCM 12-D, cu rata 4/5 şi cu 256 stări Am considerat transmisia a Q=7 biţi/interval de semnalizare 2-D Schema bloc a codorului TCM 12-D cu rata 4/5 şi 256 stări am prezentat-o în Fig 4 24 În acest caz N=6 nu este o putere întreagă a lui 2 deci am obţinut conform relaţiilor (4 125): ′  3)I(4K 0 β ==  16Q 2   14 α ==  ′ 6 )I(4K  0  ′ ==4KP  0  =2P-N  ′ Q   +4K2 ′ 0  ==24)I(4K 0 N   7911 ≅=++=09722,1)1(2NCER  0 72N2N   −11α)(1 ≅+=09733,1]-β][β-[CERCER 0  −Nα)N(12   +β)2(1 ≅=16434,2PAR  −11α)(1 (4 145) +]-β][β-[CER  0  −Nα)N(12  Constelaţia de semnale 2-D constituentă a constelaţiei de semnale 12-D este 2 imea numerelor întregi) şi este formată din un subset finit al latticei (2Z+1) (Z = mulţ 152 puncte de semnal 2-D: 7  2=128 puncte de semnal 2-D în grupul IG,  20 puncte de semnal 2-D în subgrupul OGL,  4 puncte de semnal 2-D în subgrupul OGH u Am partiţionat constelaţia de semnale 2-D (Fig 4 25) constituentă în U=2=4 subseturi 2-D: 2 ⇒ A=(4Z+1), 2 ⇒ B=(4Z+3), ⇒ C=(4Z+1)(4Z+3), ⇒ D=(4Z+3)(4Z+1) 2 Dacă constelaţia de semnale 2-D constituentă are MSED=d atunci subseturile 0 22 2-D au MSED de intrasubset egală cu d4d= se rotesc succesiv subseturile 0 Dacă 2-D în sens orar cu cîte 90° au loc următoarele transformări ale subseturilor 2-D (vezi relaţia 4 3) : A→C→B→D→A 150 -Teza de doctorat- Fig 4 24 Schema bloc a codorului TCM 12-D cu rata 4/5 şi 256 stări (construcţia propusă) 151 -Teza de doctorat- Numărul dedesubtul fiecărui punct de semnal 2-D este Z2,p, Z3,p, Z4,p, Z5,p, Z6,p, Z7,p, p=n, n+1, n+2, n+3, n+4, n+5 Fig 4 25 Constelaţia de semnale 2-D constituentă cu 152 puncte de semnal 2-D (construcţia propusă) Prin concatenarea a 6 constelaţii de semnale 2-D constituente şi prin eliminarea punctelor de semnal 12-D care conţin mai mult decît un singur punct de semnal 2-D 6x7+143 din grupul extern OG se obţine constelaţia de semnale 12-D cu 2=2 puncte de semnal 6-D 152 -Teza de doctorat- 22 Se construiesc 256 tipuri 12-D cu MSED=d4d= prin concatenarea a cîte 2 0 subseturi 6-D şi notate astfel (vezi relaţia 4 51) : (0, 0), (0, 1),…,(15, 15) unde i=0, 1,…,15 reprezintă indicele subseturilor 6-D notate Si în paragraful 4 7 3 Cele 256 tipuri 12-D sunt grupate în 32 subseturi 12-D (deci S=8) astfel încît 22 MSED a fiecărui subset 12-D obţinut să fie egală cu d4d= ia 4 52) : 0 (vezi relaţ S0=(0, 0) ∪ (2, 2) ∪ (4, 4) ∪ (6, 6) ∪ (8, 8) ∪ (10, 10) ∪ (12, 12) ∪ (14, 14) (2, 5) ∪ (4, 7) ∪ (6, 1) ∪ (8, 11) ∪ (10, 13) ∪ (12, 15) ∪ (14, 9) S1=(0, 3) ∪ (2, 6) ∪ (4, 0) ∪ (6, 2) ∪ (8, 12) ∪ (10, 14) ∪ (12, 8) ∪ (14, 10) S2=(0, 4) ∪ (2, 1) ∪ (4, 3) ∪ (6, 5) ∪ (8, 15) ∪ (10, 9) ∪ (12, 11) ∪ (14, 13) S3=(0, 7) ∪ (2, 0) ∪ (4, 6) ∪ (6, 4) ∪ (8, 10) ∪ (10, 8) ∪ (12, 14) ∪ (14, 12) S4=(0, 2) ∪ (2, 3) ∪ (4, 1) ∪ (6, 7) ∪ (8, 13) ∪ (10, 11) ∪ (12, 9) ∪ (14, 15) S5=(0, 5) ∪ (2, 4) ∪ (4, 2) ∪ (6, 0) ∪ (8, 14) ∪ (10, 12) ∪ (12, 10) ∪ (14, 8) S6=(0, 6) ∪ (2, 7) ∪ (4, 5) ∪ (6, 3) ∪ (8, 9) ∪ (10, 15) ∪ (12, 13) ∪ (14, 11) S7=(0, 1) ∪ (2, 10) ∪ (4, 12) ∪ (6, 14) ∪ (8, 0) ∪ (10, 2) ∪ (12, 4) ∪ (14, 6) S8=(0, 8) ∪ (2, 13) ∪ (4, 15) ∪ (6, 9) ∪ (8, 3) ∪ (10, 5) ∪ (12, 7) ∪ (14, 1) S9=(0, 11) ∪ (2, 14) ∪ (4, 8) ∪ (6, 10) ∪ (8, 4) ∪ (10, 6) ∪ (12, 0) ∪ (14, 2) S10=(0, 12) ∪ (2, 9) ∪ (4, 11) ∪ (6, 13) ∪ (8, 7) ∪ (10, 1) ∪ (12, 3) ∪ (14, 5) S11=(0, 15) ∪ (2, 8) ∪ (4, 14) ∪ (6, 12) ∪ (8, 2) ∪ (10, 0) ∪ (12, 6) ∪ (14, 4) S12=(0, 10) ∪ (2, 11) ∪ (4, 9) ∪ (6, 15) ∪ (8, 5) ∪ (10, 3) ∪ (12, 1) ∪ (14, 7) S13=(0, 13) ∪ (2, 12) ∪ (4, 10) ∪ (6, 8) ∪ (8, 6) ∪ (10, 4) ∪ (12, 2) ∪ (14, 0) S14=(0, 14) ∪ (2, 15) ∪ (4, 13) ∪ (6, 11) ∪ (8, 1) ∪ (10, 7) ∪ (12, 5) ∪ (14, 3) S15=(0, 9) ∪ (3, 3) ∪ (5, 5) ∪(7, 7) ∪ (9, 9) ∪ (11, 11) ∪ (13, 13) ∪ (15, 5) S16=(1, 1) ∪ (3, 4) ∪ (5, 6) ∪(7, 0) ∪ (9, 10) ∪ (11, 12) ∪ (13, 14) ∪ (15, 8) S17=(1, 2) ∪ (3, 7) ∪ (5, 1) ∪(7, 3) ∪ (9, 13) ∪ (11, 15) ∪ (13, 9) ∪ (15, 11) S18=(1, 5) ∪ (3, 0) ∪ (5, 2) ∪(7, 4) ∪ (9, 14) ∪ (11, 8) ∪ (13, 10) ∪ (15, 12) S19=(1, 6) ∪ (3, 5) ∪ (5, 3) ∪ (7, 1) ∪ (9, 15) ∪ (11, 13) ∪ (13, 11) ∪ (15, 9) S20=(1, 7) ∪ (3, 6) ∪ (5, 4) ∪ (7, 2) ∪ (9, 8) ∪ (11, 14) ∪ (13, 12) ∪ (15, 10) S21=(1, 0) ∪ (3, 1) ∪ (5, 7) ∪ (7, 5) ∪ (9, 11) ∪ (11, 9) ∪ (13, 15) ∪ (15, 13) S22=(1, 3) ∪ (3, 2) ∪ (5, 0) ∪ (7, 6) ∪ (9, 12) ∪ (11, 10) ∪ (13, 8) ∪ (15, 14) S23=(1, 4) ∪ (3, 11) ∪ (5, 13) ∪ (7, 15) ∪ (9, 1) ∪ (11, 3) ∪ (13, 5) ∪ (15, 7) S24=(1, 9) ∪ (3, 12) ∪ (5, 14) ∪ (7, 8) ∪ (9, 2) ∪ (11, 4) ∪ (13, 6) ∪ (15, 0) S25=(1, 10) ∪ (3, 15) ∪ (5, 9) ∪ (7, 11) ∪ (9, 5) ∪ (11, 7) ∪ (13, 1) ∪ (15, 3) S26=(1, 13) ∪ (3, 8) ∪ (5, 10) ∪ (7, 12) ∪ (9, 6) ∪ (11, 0) ∪ (13, 2) ∪ (15, 4) S27=(1, 14) ∪ (3, 13) ∪ (5, 11) ∪ (7, 9) ∪ (9, 7) ∪ (11, 5) ∪ (13, 3) ∪ (15, 1) S28=(1, 15) ∪ (3, 14) ∪ (5, 12) ∪ (7, 10) ∪ (9, 0) ∪ (11, 6) ∪ (13, 4) ∪ (15, 2) S29=(1, 8) ∪ (3, 9) ∪ (5, 15) ∪ (7, 13) ∪ (9, 3) ∪ (11, 1) ∪ (13, 7) ∪ (15, 5) S30=(1, 11) ∪ (3, 10) ∪ (5, 8) ∪ (7, 14) ∪ (9, 4) ∪ (11, 2) ∪ (13, 0) ∪ (15, 6) S31=(1, 12) ∪ 153 -Teza de doctorat- Dacă se rotesc succesiv subseturile 2-D în sens orar cu cîte 90° au loc următoarele transformări ale subseturilor 12-D (vezi relaţia 4 53) : Si+16→Si, i= 0, 1,…,15 Si→ 2 Subseturile 12-D sunt grupate în 2 familii 12-D cu MSED=d2 astfel încît 0 fiecare familie 12-D să conţină F=16 subseturi 12-D (vezi relaţia 4 54) : F0=S0 ∪ S2 ∪ S4 ∪ S6 ∪ S8 ∪ S10 ∪ S12 ∪ S14 ∪ S16 ∪ S18 ∪ S20 ∪ S22 ∪ S24 ∪ S26 ∪ S30 ∪ S28 ∪ S3 ∪ S5 ∪ S7 ∪ S9 ∪ S11 ∪ S13 ∪ S15 ∪ S17 ∪ S19 ∪ S21 ∪ S23 ∪ S25 ∪ S27 ∪ F1=S1 ∪ S31 ∪ S29 ∪ În acest caz avem m=4 biţi de la intrare care sunt codaţi cu un codor 8 convoluţional cu rata de codare 4/5 şi V=2=256 stă ri De asemenea, restul şi cîtul împărţirii lui Nu=12 la Q=7 sunt r=5 şi respectiv c=1 Setul stărilor curente {} a codorului convoluţional se partiţionează în F=16 nσ subseturi (vezi relaţia 4 55) : j16i≤≤≤≤+ ,15i0 15j0 {} a codorului convoluţional se partiţionează în Setul stărilor următoare {} 6nσ + F=16 subseturi (vezi relaţia 4 56) : k16j≤≤≤≤+ ,15j0 15k0 {} Conectivitatea diagramei trellis se realizează astfel: pentru orice 0≤j≤15 există tranziţii de stare între fiecare stare curentă 16i+j, 0≤i≤15 şi fiecare stare următoare 16j+k, 0≤k≤15 Asignarea subseturilor 12-D la tranziţiile de stare ale codului TCM 12-D este realizată astfel încît să fie îndeplinite trei cerinţe : 1 Subseturile 12-D, asociate cu tranziţiile care pornesc dintr-o anumită stare, să difere unele de altele şi să aparţină aceleaşi familii 12-D Acelaşi considerent se aplică şi pentru subseturile 12-D asociate cu tranziţiile care ajung într-o anumită stare Conform acestei cerinţe: ⇒ la tranziţiile care pornesc din stări pare (Y0,n=W8,n=0) sunt asignate subseturi 12-D din familia 12-D notată F0, ⇒ la tranziţiile care pornesc din stări impare (Y0,n=W8,n=1) sunt asignate subseturi 12-D din familia 12-D notată F1 2 MSED dintre două secvenţe de subseturi 12-D, care corespund la două traiectorii 222 distincte prin diagrama trellis, să fie mai mare sau egală cu 4ddd= , care este 0= free MSED a fiecărui subset 12-D (MSED de intersubset 12-D să fie mai mare sau cel 154 -Teza de doctorat- puţin egală cu MSED de intrasubset 12-D) Această cerinţă garantează că MSED dintre două secvenţe de puncte de semnal 22 12-D este d4d= i elimină evenimentele eroare care diferă cu mai mult decît un 0 ş punct de semnal 12-D faţă de secvenţa corectă de puncte de semnal 12-D Dacă se neglijează efectul de margine, atunci coeficientul de eroare Nfree al codului TCM 12-D este egal cu numărul vecinilor cei mai apropiaţi ai unui punct de semnal 12-D şi aflaţi într-un subset 12-D cu un număr finit de puncte de semnal 12-D În acest caz Nfree=6408 per punct de semnal 12-D (Nfree=1068 per punct de semnal 2-D) Dacă nu se neglijează efectul de margine, atunci coeficientul de eroare al codului TCM 12-D este mai mic decît 6408 3 Dacă X este subsetul 12-D asociat cu tranziţiile din starea curentă σ n în starea şi dacă X1 este subsetul 12-D obţinut prin rotirea lui X cu 90° atunci următoare σ 6n+ se poate defini o funcţie F1 de corespondenţă între stările codului TCM 12-D astfel fie asociat cu tranziţiile din starea curentă )(σ în starea următoare încît X1 să F n1 )(σF 6n1+ ia Pentru o rotaţie de 90° în sens orar se defineşte funcţia F1 astfel (vezi relaţ 4 57) : n8,n7,n6,n5,n4,n3,n2,n1,n8,n7,n6,n5,n4,n3,n2,n1,WWWWWWWWWWWWWWWW→ unde bara deasupra înseamnă inversare binară şi W1,n, W2,n ,…, W8,n sunt biţ ii de stare curentă a codorului convoluţional Acestă cerinţă garantează că vom obţine un cod TCM 12-D transparent la toate ambiguităţile de fază ale constelaţiei de semnale 12-D, cu alte cuvinte un cod TCM 12-D invariant rotaţional la o rotaţie de 0°, 90°, 180°, 270° Pentru a suprima toate cele 4 ambiguităţi de fază a constelaţiei de semnale 12-D (0°, 90°, 180°, 270°) se foloseşte un codor diferenţial pe 2 biţi (vezi relaţia 4 58) : ′′′′ )mod4IIII()II(+ n5,n4,6-n5,6-n4,n5,n4,= deci în acest caz avem a=5 şi b=0 Convertorul de bit converteşte cei 12 biţi de la intrarea sa Y0,, I1,n, I2,n, I3,n, I′ n4,n, I5,n, I6,, I7,n, I1,n+1, I2,n, I3,n, I4 în 6 grupuri de cîte 2 biţ ia 4 59) n+1+1,n+1i fiecare (vezi relaţ : Z0,p, Z1,p, p=n, n+1, n+2, n+3, n+4, n+5 care selectează 6 subseturi 2-D corespunzătoare tipului 12-D selectat în cadrul subsetului 12-D utilizat 155 -Teza de doctorat- Pentru a obţine un cod TCM 12-D transparent la orice ambiguitate de fază a constelaţiei de semnale 12-D trebuie sa fie îndeplinită cerinţa prezentată în continuare Dacă X este tipul 12-D asociat cu structura de biţi ′ ,Y′ ,I ,I ,I,II,I ,I ,I ,I ,I I atunci se 1n4,1n3,1n2,1n1,n7,n6,n5,n4,n3,n2,n1,n,0++++ notează cu X1, X2, X3 tipurile 12-D obţ inute prin rotirea lui X cu 90°, 180°, 270° ′′′′′ i obţinute Dacă ,I(′ )I 1n5,n4,, ,I(, 2n5,n4,)I ,I( 2n5,n4,)I sunt perechile de biţ ′ prin translaţia perechii de biţi ,I(′ )I sau trei poziţii în secvenţa n5,n4, cu una, două circulară 00→10→01→11→00, atunci tipurile 12-D asociate cu structurile de biţi: ′ ,Y′ ,I ,I ,I,)II(,I ,I ,I ,I ,I I , 1n4,1n3,1n2,1n1,n7,n6,1n5,n4,n3,n2,n1,n,0++++ ′ ,Y′ ,I ,I ,I,)II(,I ,I ,I ,I ,I I , 1n4,1n3,1n2,1n1,n7,n6,2n5,n4,n3,n2,n1,n,0++++ ′ ,Y′ ,I ,I ,I,)II(,I ,I ,I ,I ,I I , 1n4,1n3,1n2,1n1,n7,n6,3n5,n4,n3,n2,n1,n,0++++ sunt X1, X2, X3 Corespondenţa dintre cei 2 biţi Z0,p, Z1,p, p=n, n+1, n+2, n+3, n+4, n+5 ş i subseturile 2-D este prezentată în Tabelul 4 1 iar funcţionarea convertorului de bit este prezentată în Tabelul 4 11 Am notat cu Si subseturile 12-D ş i astfel indexul i este furnizat de relaţia (vezi relaţia 4 60) : ′ +++=, i = 0, 1,…,31 n4,n3,n2,n1,n0,II2I4I8Y16i+ Codorul bloc foloseşte cei 31 biţi I5,n, I6,n+1, I7,n+1, I1,n+2, I2,n, I3,n I4, I5,n, +1+2+2,n+2+2 I6,n+2, I7, I1,n+3, I2,n, I3,n+3 I4,n+3, I5,n+3, I6,n+3, I7,n, I1,n+4, I2, I3,n+4 I4,, I5,n+4, I6, ,n+2+3+3,n+4n+4,n+4 I7,n+4, I1, I2,n+5, I3, I4,n, I5,n+5, I6,n+5, I7,n+5 pentru a genera 6 grupuri de cîte 6 biţ ,n+5n+5+5i fiecare: Z2,p, Z3,p, Z4,p, Z5,p, Z6,p, Z7,p, p=n, n+1, n+2, n+3, n+4, n+5 (4 146) care selectează pentru fiecare din cele 6 subseturi 2-D utilizate, punctul de semnal 2-D care va fi transmis Funcţionarea codorului bloc este prezentată în Tabelul 4 20 Cîştigul de codare asimptotic al codului TCM 12-D faţă de o transmisie necodată este conform relaţiei (4 134): 4  61722,5log10γ≅ = [dB] (4 147)  10 098,1   Tabelul 4 1 Corespondenţa dintre cei 2 biţi Z0,p, Z1,p ş i subseturile 2-D Subset 2-D Z0,p Z1,p A 0 0 B 0 1 C 1 0 D 1 1 156 -Teza de doctorat- Tabelul 4 11 Convertorul de bit 5 Z1,n+ 0 1 0,n+5 Z 1 0 0 0,n+ 0 4 Z 0 0 0 1 0 0,n+ 1 4 Z 1 0 0 0 0 0 1 1,n+ 0 3 0 Z 0 0 0 10 0 1 0 0,n+ 3 10 Z 0 0 0 010 0 1 0 010 Z1,n+2 0 0 0 1000 0 1 0 1000 0,n+2 0 Z 0 0 0100 00 1 0 0100 00 1,n+1 0 Z 0 0001 00 10 0 1001 00 00 0,n+1 0 Z 1101 10 10 00 0101 10 00 00 1,n 0001 Z 10 11 00 1001 00 0 01 10 1101 Z0,n 10 0 11 10 0101 10 0 01 0 10 0000 4,n+1 10 0 I 10 0 10 1000 00 0 0 00 0 10 1100 10 1 0 3,n+1 10 0 I 10 0100 10 0 0 00 0 0 10 0000 10 1 0 11 0 0 2,n+1 10 I 1000 00 0 1 01 1 0 10 0 1100 11 1 1 11 1 0 10 0 1,n+1 0100 I 11 0 1 01 0 0 10 0 0 0 1011 1 1 10 0 0 10 0 0 0 0011 0 1 Tip 00 1 0 10 1 12-D 1 0 0 1110 11 10 11 1 1 0 0 (0,0)7,n 0110 00 10 I 11 1 0 0 0 1 1 1111 101 0 0 0 0 0 1 0 6,n 0111 010 0 0 I 1 0 0 1 0 1 1 0 111 100 00 0 1 0 5,n 1 0 0 011 001 00 I′ 0 0 0 1 1 0 0 0 0 01 01 0 0 0 0 1 0 0 4,n 0 1 0 01 11 0 0 I′ 0 0 1 1 0 0 1 0 1 1 0 0 0 1 0 0 1 0 0 0 0 0 1 1 0 3,n 0 I 0 1 0 0 1 1 0 1 0 0 1 1 0 1 00 1 0 10 0 1 1 0 1 0 0 1 00 2,n 0 1 10 0 I 1 1 0 1 0 0 1 00 1 00 10 000 1 0 0 0 00 1,n 0 10 10 000 I 1 0 0 0 0 1 0 0 1 011 0 0000 0 0 1 0 0 0,n 0 111 0 0000 Y 0 0 1 0 0 (2,2) 0 0 0 0 0 1 0 0 00 0 1 0 0 1 0 0 0 1 1 0 0 00 0 0 0 0 0 0 0 Subset 1 1 0 0 0 00 0 0 12-D 1 0 00 1 1 0 0 0 00 0 0 0 0 00 0 1 0 0 0 0 00 1 0 1 0 0 10 0 0 0 00 1 0 0 0 111 0 0 100 157 0 1 100 1 0 0 00 1 0 100 1 0 0 00 1 0 0 1 0 0 1000 1 0 0 0 0 1000 1 0 0 0 1 0 00 100 0 0 0 00 100 0 1000 0 0 0 0 1000 0 0 0 0 1000 00 0 … 0 1000 0 0 0 0 000 0 0 0 0 0 0 000 0 0 0 0 00 0 0 0 0 00 0 0 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -Teza de doctorat- Tabelul 4 20 Codorul bloc 7,n+5 Z7,n+57,n+57,n+57,n+57,n+57,n+57,n+57,n+57,n+57,n+57,n+57,n+57,n+57,n+5 I IIIIIIIIIII0 I I6,n+5 Z6,n+56,n+56,n+56,n+56,n+56,n+56,n+56,n+56,n+56,n+56,n+56,n+56,n+56,n+5 I IIIIIIIIIII0 I I5,n+5 Z5,n+55,n+55,n+55,n+55,n+55,n+55,n+55,n+55,n+55,n+55,n+55,n+55,n+5 I IIIIIIIIII0 I 4,n+51 Z4,n+54,n+54,n+54,n+54,n+54,n+54,n+54,n+54,n+54,n+54,n+54,n+54,n+5 I IIIIIIIIII0 I I 3,n+50 Z3,n+53,n+53,n+53,n+53,n+53,n+53,n+53,n+53,n+53,n+53,n+53,n+53,n+5 I IIIIIIIIII0 I I 2,n+50 Z 0 I 0 7,n+4 Z2,n+52,n+52,n+52,n+52,n+52,n+52,n+52,n+52,n+52,n+55,n+57,n+52,n+52,n+5 I IIIIIIIII0 III 0 I6,n+4 1,n+51,n+51,n+51,n+51,n+51,n+51,n+51,n+51,n+51,n+54,n+56,n+51,n+51,n+5 0 Z I IIIIIIIII0 III 0 I5,n+4 Z7,n+47,n+47,n+47,n+47,n+47,n+47,n+47,n+47,n+43,n+55,n+57,n+47,n+4 I IIIIIIII0 III 0 4,n+41 0 Z6,n+46,n+46,n+46,n+46,n+46,n+46,n+46,n+46,n+42,n+54,n+56,n+46,n+4 I IIIIIIII0 III I 0 3,n+40 Z5,n+45,n+45,n+45,n+45,n+45,n+45,n+45,n+45,n+41,n+53,n+55,n+45,n+4 I 0 IIIIIIII0 III I 0 Z2,n+40 0 0 I 0 7,n+3 Z4,n+44,n+44,n+44,n+44,n+44,n+44,n+44,n+47,n+42,n+57,n+42,n+54,n+44,n+4 1IIIIIIII0 IIIII 0 1 I6,n+33,n+43,n+43,n+43,n+43,n+43,n+43,n+43,n+46,n+41,n+56,n+41,n+53,n+43,n+4 0 Z I IIIIIII0 IIIII 0 0 I5,n+3 Z2,n+42,n+42,n+42,n+42,n+42,n+42,n+45,n+47,n+45,n+47,n+42,n+42,n+4 0IIIIIII0 IIIII 0 4,n+31 0 Z1,n+41,n+41,n+41,n+41,n+41,n+41,n+44,n+46,n+44,n+46,n+41,n+41,n+4 I IIIIII0 IIIII I 0 3,n+30 Z7,n+37,n+37,n+37,n+37,n+37,n+37,n+33,n+45,n+43,n+45,n+47,n+37,n+3 I 0 IIIIII0 IIIII I 1 Z2,n+30 0 1 I 0 7,n+2 Z6,n+36,n+36,n+36,n+36,n+36,n+32,n+44,n+42,n+44,n+42,n+44,n+46,n+36,n+3 0IIIIII0 IIIIIII 0 0 I6,n+25,n+35,n+35,n+35,n+35,n+35,n+31,n+43,n+41,n+43,n+41,n+43,n+45,n+35,n+3 0 Z I IIIII0 IIIIIII 0 0 I5,n+2 Z4,n+34,n+34,n+34,n+34,n+37,n+32,n+47,n+32,n+47,n+32,n+44,n+34,n+3 0IIIII0 IIIIIII 0 1 0 158 I 1 1 0 0 0 0 0 0 -Teza de doctorat- 4,n+2 Z3,n+33,n+33,n+33,n+33,n+36,n+31,n+46,n+31,n+46,n+31,n+43,n+33,n+3 I IIII0 IIIIIII 3,n+20 Z2,n+32,n+32,n+32,n+32,n+35,n+37,n+35,n+37,n+35,n+37,n+32,n+32,n+3 I IIII0 IIIIIII I 2,n+20 Z 0 I 0 7,n+1 Z1,n+31,n+31,n+31,n+34,n+36,n+34,n+36,n+34,n+36,n+34,n+36,n+31,n+3 I III0 IIIIIIII1 0 I6,n+1 7,n+27,n+27,n+27,n+23,n+35,n+33,n+35,n+33,n+35,n+33,n+35,n+37,n+2 0 Z I III0 IIIIIIII1 0 I5,n+1 Z6,n+26,n+26,n+22,n+34,n+32,n+34,n+32,n+34,n+32,n+34,n+36,n+2 I II0 IIIIIIII1 4,n+11 1 Z5,n+25,n+25,n+21,n+33,n+31,n+33,n+31,n+33,n+31,n+33,n+35,n+2 I II0 IIIIIIII1 I 0 3,n+10 Z4,n+24,n+24,n+27,n+22,n+37,n+22,n+37,n+22,n+37,n+22,n+34,n+2 I 0 II0 IIIIIIII1 I 0 Z2,n+10 0 0 I 0 7,n Z3,n+23,n+26,n+21,n+36,n+21,n+36,n+21,n+36,n+21,n+36,n+21,n+31,n+3 0 I0 IIIIIIIII1 I 0 0 II6,n2,n+22,n+25,n+27,n+25,n+27,n+25,n+27,n+25,n+27,n+25,n+27,n+27,n+2 1 Z II0 IIIIIIIII1 0 1 II5,n Z1,n+24,n+26,n+24,n+26,n+24,n+26,n+24,n+26,n+24,n+26,n+26,n+2 0 0 IIIIIIIII1 0 I 1 I4,n 0 Z7,n+13,n+25,n+23,n+25,n+23,n+25,n+23,n+25,n+23,n+25,n+25,n+2 I0 IIIIIIIII1 I 0 0 I3,n Z6,n+12,n+24,n+22,n+24,n+22,n+24,n+22,n+24,n+22,n+24,n+24,n+2 0 I0 IIIIIIIII1 I 0 0 IZ2,n 0 0 I 1 I3,n+2 0 X 1 0 X 2,n+2 I X 0 0 0 X 0 1,n+2 I 1 XX 0 0 1 0 I7,n+1 XX 0 X 0 00 0 6,n+1 0 I 1 X X1 0 X 0 0 5,n+1 10 0 I 0 1 0 1 X1 0 1 X 0 0 11 0 1 0 1 1 X1 0 1 159 X 0 0 01 1 1 1 0 1 X1 0 1 0 1 0 0 0 1 1 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -Teza de doctorat- CAPITOLUL V Modemuri analogice de bandă vocală cu modulaţie codată trellis multidimensională 5 1 Schema bloc funcţională Schema bloc funcţională a unui modem analogic de bandă vocală cu TCM multidimensională 2N-D am prezentat-o în Fig 5 1 5 1 1 Cifratorul şi decifratorul Cifratorul (scrambler) are două funcţii: ¾ facilitarea regenerării semnalului de tact la recepţie, ¾ distribuirea uniformă a energiei semnalului modulat QAM (Quadrature Amplitude Modulation) în întreaga banda de frecvenţe a canalului de comunicaţii Cifratorul este un registru de deplasare liniar cu reacţie înapoi şi este definit de un polinom generator (generating polynomial) : -18-23 ⇒ pentru modemul chemător: GPC=1+x+x (Generating Polynomial Calling modem), -5-23 ⇒ pentru modemul chemat: GPA=1+x+x (Generating Polynomial Answering modem) Cifratorul aleatorizează secvenţa binară de date de intrare {bn} adică converteşte această secvenţă într-o secvenţă binară pseudoaleatoare {cn} cu perioada 23 2-1=8 388 607 biţ i O astfel de secvenţă binară pseudoaleatoare lungă nu necesită utilizarea unui polinom de gardă pentru a preveni apariţia la ieşirea cifratorului a unor structuri de bit repetate La transmisie, polinomul care reprezintă secvenţa de mesaj, se divide la polinomul generator GPC/GPA iar coeficienţii cîtului luaţi în ordine descrescătoare formează secvenţa de date transmisă (Fig 5 2) La recepţie, polinomul care reprezintă secvenţa de date recepţionată, se multiplică cu polinomul generator GPC/GPA iar coeficienţii produsului luaţi în ordine descrescătoare formează secvenţa de mesaj (Fig 5 3) 160 -Teza de doctorat- Fig 5 1 Schema bloc funcţională a unui modem analogic de bandă vocală cu TCM multidimensională 2N-D 161 -Teza de doctorat- Fig 5 2 Cifrator (scrambler) Fig 5 3 Decifrator (descrambler) 5 1 2 Convertorul serie/paralel şi convertorul paralel/serie La transmisie, convertorul serie/paralel (S/P), care este un registru de deplasare de lungime Q biţi, realizează conversia serie/paralel a secvenţei binare {cn} de la ieşirea cifratorului 162 -Teza de doctorat- Pentru fiecare interval de semnalizare 2-D (interval de simbol 2-D) dintr-un grup de N intervale de semnalizare 2-D, distribuitorul asigură transmisia celor Q biţi de la ieşirea convertorului S/P spre cele Q intrări corespunzătoare ale codorului TCM 2N-D (Fig 5 4) La recepţie, pentru fiecare interval de semnalizare 2-D dintr-un grup de N intervale de semnalizare 2-D, distribuitorul asigură transmisia celor Q biţi de la ieşirile corespunzătoare ale decodorului TCM 2N-D spre cele Q intrări ale convertorului paralel/serie (P/S) Convertorul P/S, care este un registru de deplasare de lungime Q biţi, realizează conversia paralel/serie a secvenţei binare de la ieşirea decodorului TCM 2N-D (Fig 5 5) Fig 5 4 Convertor S/P + distribuitor Fig 5 5 Convertor P/S + distribuitor 5 1 3 Codorul neliniar Distorsiunile neliniare sunt cauzate de cuantizarea de la codarea PCM (Pulse Code Modulation) a semnalelor analogice de bandă vocală la intrarea într-o porţiune digitală a PSTN şi/sau de unele componente analogice (bobine de încărcare, transformatoare de linie) prezente în PSTN Întrucît distorsiunile neliniare afectează punctele de semnal 2-D din OG mai mult decît punctele de semnal 2-D din IG, codarea neliniară (warping) măreşte imunitatea la zgomot a punctelor de semnal 2-D din OG şi micşorează imunitatea la 163 -Teza de doctorat- zgomot a punctelor de semnal 2-D din IG Cu alte cuvinte, codarea neliniară întinde porţiunea OG şi strînge porţiunea IG a constelaţiei de semnale 2-D (Fig 5 6) În consecinţă, se măreşte distanţa euclidiană medie dintre punctele de semnal 2-D din OG şi se micşorează distanţa euclidiană medie dintre punctele de semnal 2-D din IG Fig 5 6 Efectul codării neliniare asupra constelaţiei de semnale 2-D Secvenţa de intrare {x(n)} este codată neliniar conform relaţiei : x′(n)=Φ(n) x(n) (5 1) unde: 2 (n)ξξ(n) +=Φ este funcţia de proiecţie neliniară, 1(n)+ 1206 22 (n)]x(n)[x+ ir δξ(n)=, (5 2) 22 (n)]x(n)[x+ ir 22 este energia medie a secvenţei de intrare {x(n)}, (n)]x(n)[x+ ir δ este o constantă care se selectează în faza 4 a procedurii de negociere dintre modemul local şi modemul distant şi poate fi 0 sau 0,3125 5 1 4 Filtrul de preaccentuare În faza 2 a procedurii de negociere dintre modemul local şi modemul distant se determină : ⇒ rata de bit=Db [bit/s] ⇒ rata de simbol=vm [Baud] sau [simboluri/s]: 164 -Teza de doctorat- a 0,01%2400v± (5 3) m= c unde a şi c sunt constante (Tabelul 5 1) ⇒ frecvenţa semnalului purtător=fp [Hz]: d (5 4) mpvf= e unde d şi e sunt constante (Tabelul 5 2) Tabelul 5 1 vm [Baud] a c 2400 1 1 2743 8 7 2800 7 6 3000 5 4 3200 4 3 3429 10 7 Tabelul 5 2 vm [Baud] Joasă frecvenţă Înaltă frecvenţă fp [Hz] d e fp [Hz] d e 2400 1600 2 3 1800 3 4 2743 1646 3 5 1829 2 3 2800 1680 3 5 1867 2 3 3000 1800 3 5 2000 2 3 3200 1829 4 7 1920 3 5 3429 1959 4 7 1959 4 7 Filtrul de preaccentuare folosit este specificat printr-un index numeric i=0,1,2,…,10 care este furnizat de modemul distant în faza 2 a procedurii de negociere Amplitudinea spectrului de putere a semnalului transmis s′(t) trebuie să se încadreze între limitele precizate de graficele următoare pentru o frecvenţa normalizată ddf 0,45;[+−∈ 0,45 ] (Fig 5 7 şi Tabelul 5 3) eev m 165 -Teza de doctorat- Amplitudinea [dB] α 0 0 0 2 0 4 0 6 0 8 1 0 1 2 Frecvenţa normalizată f/vm Amplitudinea [dB] 2β β 0 0 0 2 0 4 0 6 0 8 1 0 1 2 Frecvenţa normalizată f/vm NOTĂ – Toleranţa pentru amplitudinea spectrului de putere a semnalului transmis este ± 1 dB Fig 5 7 Amplitudinea spectrului de putere a semnalului transmis 166 -Teza de doctorat- Tabelul 5 3 Indexul numeric al filtrului de preaccentuare Index numeric α [dB] β [dB] 0 0 - 1 2 - 2 4 - 3 6 - 4 8 - 5 10 - 6 - 0,5 7 - 1 8 - 1,5 9 - 2 10 - 2,5 5 1 5 Egalizarea adaptivă Egalizarea este procedura prin care se compensează distorsiunile de amplitudine şi de fază ale semnalului util introduse de canalul de comunicaţie Întrucît caracteristicile canalului de comunicaţie variază în timp se foloseşte o procedură de egalizare adaptivă care asigură o calitate constantă a transmisiei de date Se foloseşte o schemă de egalizare adaptivă neliniară cu precodare Tomlinson-Harashima (THP=Tomlinson Harashima Precoding) (Fig 5 8) [46, 47, 48] Ideea este de împărţi egalizarea adaptivă între transmiţător şi receptor Receptorul de date calculează coeficienţii optimali ai egalizorului adaptiv cu reacţie decizională (DFE=Decision Feedback Equalizer) în faza 2 a procedurii de negociere, după care îi transmite la transmiţător în faza 4 a procedurii de negociere Transmiţătorul de date utilizează aceşti coeficienţi de filtrare optimali pentru a precoda (preegaliza) semnalul util înainte de transmisia prin canalul de comunicaţie 167 -Teza de doctorat- Fig 5 8 Schema de egalizare adaptivă neliniară cu precodare La transmisie secvenţa de date {d(n)} de la ieşirea codorului TCM 2N-D, este precodată astfel încît la ieşirea precodorului se obţine secvenţa transmisă {x(n)} La fiecare moment de timp discret n, se calculează simbolul f(n) prin scăderea din simbolul de date curent d(n) a contribuţiei ISI a simbolurilor transmise anterioare {x(n-i), i≥1} : i)x(nhd(n)f(n)− = (5 5) ∑ i− 1i≥ -1-2 unde H(z)=1+h1z+h2z+… este funcţ ia de transfer a canalului de comunicaţie modelat drept un canal AWGN cu ISI Simbolul transmis curent x(n) se obţine prin aplicarea operaţiei modulo-2R [46, 47]: (5 6) −=−= a(n)i)x(nhd(n)a(n)f(n)x(n)− ∑ i− 1i≥ 2 i simbolul a(n)∈2RZ unde R este partea întreagă din R2 (R=[R2]) ş Dacă se aplică transformata z se obţine : X(z)=F(z)-A(z)=D(z)-A(z)-X(z)[H(z)-1]⇔X(z)H(z)=D(z)-A(z) (5 7) 168 -Teza de doctorat- Dar X(z)H(z)=Y(z) deci rezultă că : Y(z)=D(z)-A(z) (5 8) În absenţa zgomotului secvenţa de la ieşirea canalului de comunicaţie este egală cu secvenţa de date modificată : 2 y(n)=d(n)-a(n), a(n)∈2RZ (5 9) Filtrul cu funcţia de transfer H(z), care modelează canalul de comunicaţie, include efectele codorului neliniar, modulatorului QAM, filtrului de preaccentuare, supresorului de ecou, demodulatorului QAM, blocului AGC (Automatic Gain Control) şi filtrului de interpolare/decimare La recepţie, egalizorul adaptiv DFE îndepărtează ISI iar filtrul de înălbire a zgomotului (whitening noise filter) (care de fapt este un filtru de predicţie a erorii) reduce varianţa zgomotului de la ieşirea egalizorului adaptiv DFE prin cîştigul de predicţie: +π 11 (5 10) = pdfγ ∫ 2 2π H(f) −π unde H(f) este funcţia de transfer a filtrului care modelează canalul de comunicaţie Atenuarea introdusă la marginile benzii de frecvenţe este: H(π ) (5 11) = α H(0) La ieşirea filtrului de înălbire a zgomotului, se obţine secvenţa recepţionată {r(n)}, care conţine un zgomot aditiv, alb, gaussian: R(z)=Y(z)+Z(z) (5 12) Urmează blocul care realizează operaţia modulo-2R şi decodorul TCM 2N-D, pentru estimarea secvenţei de date (n)}dˆ{ ş i a biţilor de informaţie originali Egalizorul adaptiv DFE se bazează pe următorul principiu = odată ce un simbol recepţionat din canalul de comunicaţie fost estimat (detectat), contribuţia ISI a acestui simbol la simbolurile recepţionate ulterior, poate fi estimată şi îndepărtată înainte de estimarea (detecţia) simbolurilor recepţionate ulterior Schema bloc a egalizorului adaptiv DFE este prezentată în Fig 5 9 169 -Teza de doctorat- Fig 5 9 Schema bloc a egalizorului adaptiv DFE Am notat: ⇒ {u(n)}=secvenţa recepţionată din canalul de comunicaţie, ⇒ {v(n}=secvenţa egalizată, a de antrenare în ⇒ {w(n)}=secvenţa estimată (n)}vˆ{ în modul normal sau secvenţ modul antrenare, ⇒ {e(n)}=secvenţa eroare de estimare: e(n)=w(n)-v(n) (5 13) Înaintea de transmisia de date propriu-zisă nu se cunosc caracteristicile canalului de comunicaţie Din aceste motiv egalizorul adaptiv DFE trece în modul antrenare (training mode), în care se transmite o secvenţă de antrenare (training sequence), pentru a colecta informaţii despre caracteristicile canalului de comunicaţie Astfel se colectează o secvenţă eroare de estimare {e(n)}, care se foloseşte la ajustarea sincronă a coeficienţilor de filtrare ai filtrului cu reacţie cu conexiune înainte şi a filtrului cu reacţie cu coenxiune înapoi, în scopul minimizării erorii de estimare medii pătratice (MSEE=Mean Squared Estimation Error) Circuitul de decizie (Slicer) estimează, pentru fiecare simbol din secvenţa egalizată {v(n)}, care a fost simbolul cel mai probabil transmis dintr-un set disponibil de simboluri Filtrul cu reacţie cu conexiune înainte (Feedforward Filter=FFF) filtrează secvenţa recepţionată din canalul de comunicaţie {u(n)}, astfel încît să se obţină o eroare de estimare medie pătratică minimă, deci secvenţa estimată (n)}vˆ{ să fie corectă cu o probabilitate maximă Filtrul cu reacţie cu conexiune înapoi (Feedbackward Filter=FBF) înlătură contribuţia ISI a fiecărui simbol estimat la simbolurile recepţionate ulterior 170 -Teza de doctorat- 5 1 6 Modulatorul QAM Demodulatorul QAM Blocul de refacere a semnalului purtător Modulaţia de amplitudine în cuadratură (Quadrature Amplitude Modulation) este o tehnică de modulaţie care combină modulaţia de amplitudine şi modulaţia de fază, în sensul că informaţia pe care dorim să o transmitem prin canalul de comunicaţie se reprezintă prin variaţii de amplitudine şi schimbări de fază ale unui semnal purtător Schema bloc a modulatorului QAM este următoarea: Fig 5 10 Modulatorul QAM La momentul de timp discret n, simbolul de informaţie complex: ′′ ′ (n)jx(n)x(n)x+ (5 14) ir= ′′′′ reprezintă punctul de semnal 2-D cu coordonatele x(n),(x, (n)) unde x(n),x (n) irir sunt numere reale La intrarea filtrelor formatoare de impulsuri (Pulse Shaping Filter) avem : ′′ = (n)x(t)x=semnalul (t modulator în fază ∑ srr)nT-δ n ′′ = (n)x(t)x=semnalul (t modulator în cuadratură (5 15) ∑ sii)nT-δ n t1,=0  unde δ(t) este impulsul Dirac= iar Ts este durata simbolului de informaţ ie  t0,≠0  171 -Teza de doctorat- La ieşirea filtrelor formatoare de impulsuri (Pulse Shaping Filter) avem : ′ = ∑ (n)xI(t) =semnalul sr)nT-g(t modulator în fază n ′ = (5 16) ∑ (n)xQ(t) =semnalul si)nT-g(t modulator în cuadratură n unde g(t) este răspunsul la impuls a filtrelor formatoare de impulsuri (care sunt identice) Semnalul modulat QAM se obţine prin înmulţirea semnalelor modulatoare în fază, respectiv în cuadratură, de la ieşirea filtrelor formatoare de impulsuri, cu componentele în fază, respectiv în cuadratură ale semnalului purtător şi prin însumarea semnalelor rezultante : (5 17) I(t)(t)s(t)ss(t)+ cos Q(t)t)(ω t)(ωsin += ppQI= unde 2π f este pulsaţia semnalului purtător = ppω Componentele în fază şi în cuadratură ale semnalului purtător sunt generate de un oscilator local şi un defazor cu (-90°) Dacă răspunsul la impuls al filtrelor formatoare de impulsuri ar fi : T s 1,≤ t  2t TG(f))rect(g(t)=↔== )sinc(fT (5 18) ss TT ss 0,> t  2 atunci durata răspunsului la impuls g(t) este strict limitată la durata Ts a simbolului de informaţie dar spectrul de putere al semnalului modulat QAM ocupă o bandă de frecvenţe infinită 1 Întrucît funcţia sinc(fTs) scade lent cu creş terea frecvenţei f (cu rata ≅ ), f rezultă că este necesară o lăţime de bandă B foarte mare a canalului de comunicaţie pentru a conţine cea mai mare parte din puterea semnalului modulat QAM De exemplu, pentru ca banda de frecvenţe B să conţină 90% din puterea semnalului modulat QAM este necesar ca: 2 B≥ (5 19) sT Dacă răspunsul la impuls al filtrelor formatoare de impulsuri ar fi : 172 -Teza de doctorat- 1  ≤ 2π sf,T   s2Tt G(f))sinc(g(t) (5 20) ==↔=  2π T ss)rect(fT s1T  0, f >  s2T  atunci spectrul de putere al semnalului modulat QAM este strict limitat la lăţimea de 1 pe care o oferă canalul de comunicaţie, dar durata răspunsului la impuls bandă B= sT g(t) este infinită t ia valori nule la momentele de timp Întrucît răspunsul la impuls )sinc(g(t)= T s egale cu multipli întregi ai Ts, rezultă că pentru a obţine ISI=0 vom eşantiona semnalul Z modulat QAM la momente de timp t=kTs, unde k∈ Practic este imposibil de a realiza şi de a menţine o sincronizare perfectă la recepţie În consecinţă, fiecare eşantion ale semnalului demodulat QAM, va conţine urme ale eşantioanelor anterioare, sub forma unor sume de eşantioane ale cozilor din jurul impulsurilor principale şi deci ISI≠0 Deoarece sunt dificil de realizat practic filtre formatoare de impulsuri cu un t , se folosesc filtre formatoare de impulsuri de tip răspuns la impuls )sinc(g(t)= T s cosinus ridicat (Raised Cosine Pulse Shaping Filter) cu răspunsul la impuls : tπ tπ     )cos()sin(     TT ss     = (5 21) ,g(t) 0 atunci eş antionarea simbolului are loc prea tîrziu şi dacă en<0 atunci eşantionarea simbolului are loc prea devreme La convergenţa buclei de refacere a 0 şi eşantionarea simbolului are loc la momentele de timp tactului de simbol, en→ corecte Filtrul de buclă este identic cu cel folosit în bucla PLL din blocul de refacere a semnalului purtător şi este un filtru activ de ordinul doi cu funcţia de transfer : gsg+ ip H(s)= (5 29) s unde avem: =cîştigul proporţional, 2ξg= npω 2 g ==cîştigul integrator, niω ξ=factorul de amortizare, egal de obicei cu 0,707, BL=banda de frecvenţ e a buclei, egală cu 0,05 şi care poate fi redusă după convergenţa buclei de refacere a tactului de simbol, 8ξ LB = =frecvenţa naturală a buclei nω 2 4ξ1+ 177 -Teza de doctorat- Fig 5 14 Filtrul de buclă Filtrul de buclă este format din două ramuri (Fig 5 14): -1) ramura proporţională care multiplică eroarea de fază a tactului de simbol cu un cîştig proporţional gp ş i asigură urmărirea erorii de fază a tactului de simbol, -2) ramura integratoare care multiplică eroarea de fază a tactului de simbol cu un cîştig integrator gi ş i apoi integrează semnalul rezultat Se asigură urmărirea erorii de frecvenţă a tactului de simbol 5 1 8 Blocul de control automat a cîştigului Blocul de control automat al cîştigului (AGC=Automatic Gain Control) este un sistem adaptiv care are rolul de a menţine semnalul de la ieşirea sa la un nivel mediu de energie constant Acest lucru este absolut necesar pentru funcţionarea corectă a blocului de refacere a tactului de simbol Schema bloc a blocului de control automat a cîştigului este prezentată în Fig 5 15 Fiecare simbol de la intrare x′(n), este înmulţit cu un factor de cîştig α(n), iar rezultatul este furnizat la ieşire şi este folosit pentru actualizarea detectorul de energie totală din intervalul de baud Energia totală din intervalul de baud se calculează cu relaţia : ′ = ∑ totala(n)xE (5 30) intervalul de baud Energia totală din intervalul de baud este comparată cu un nivel de referinţă AGC şi se calculează o eroare de buclă AGC care este minimizată de blocul AGC 178 -Teza de doctorat- Fig 5 15 Blocul de control automat a cîştigului Energia semnalului modulat QAM pentru un punct de semnal 2-D, cu ′′ coordonatele (n)) x(n),(x, este egală ir cu : T 2π s 1 ′′ 22 dtt)](n)sin(ωxt)(n)cos(ω[x(t)dtsE(n)= == pipr+ ∫∫ 2π 00 2π2π t)cos(2ω1t)cos(2ω + 111pp− ′ 22′ = (5 31) ]dt[(n)][x]dt[(n)][x+ ir+ ∫∫ 22π 22π 00 ′′ 222π (n)][x(n)][x1+ ′ ir′ t)]dtω(n)][sin(2(n)x[x= pir+ ∫ 22π 0 Nivelul de referinţă AGC se calculează cu relaţia : )}n(E{EE= ref (5 32) unde E{ } este operatorul de mediere statistică În continuare se aplică o fereastră (fereastra 1) în jurul nivelului de referinţă AGC Astfel blocul AGC va aplica corecţii numai dacă eroarea de buclă AGC se află în exteriorul acestei ferestre De asemenea se aplică o pondere (în general egală cu 2) erorii de buclă AGC astfel încît să se realizeze o relaţie de paritate între erorile pozitive şi negative 179 -Teza de doctorat- Detectorul de mod de operare AGC funcţionează astfel: ⇒ atît timp cît eroarea de buclă AGC se află în interiorul ferestrei 2, blocul AGC operează în modul urmărire (tracking mode) şi compensează schimbările din nivelul mediu de energie al semnalului de intrare, prin înmulţirea acestuia cu un factor de cîştig, care este ajustat în mod corespunzător de către integratorul exponenţial, ⇒ în momentul în care eroarea de buclă AGC iese din interiorul ferestrei 2, blocul AGC trece în modul urmărire rapidă (fast tracking mode), iar în momentul în care eroarea de buclă AGC revine în interiorul ferestrei 2, blocul AGC trece din nou în modul urmărire (Fig 5 16) Integratorul exponenţial calculează factorul de cîştig conform relaţiei : α(n)1)α(n−=+ α(n)e(n)K (5 33) unde: e(n) este eroarea de buclă AGC, K este o constantă care determină viteza de convergenţă a buclei AGC În general, constanta K este egală cu : 1  =rapidã,K modulpentru urmãrire   2 (5 34)  1  <<urmãrire,K0 modulpentru   2 Fig 5 16 Modurile de operare a blocului AGC 180 -Teza de doctorat- 5 1 9 Supresorul de ecou Schema bloc a conexiunii telefonice prin PSTN dintre un modem local şi un modem distant este prezentată în Fig 5 17 Fig 5 17 Schema bloc a conexiunii telefonice prin PSTN dintre un modem local şi un modem distant Sistemul diferenţial 2F/4F realizează conversia de la o conexiune pe 2 fire (2F) în linia telefonică la o conexiune pe 4 fire (4F) în PSTN sau în modem Sistemele diferenţiale 2F/4F introduc neadaptări de impedanţă în punctele în care sunt utilizate Din acest motiv, dacă ne referim la modemul local, o parte din semnalul transmis de acesta s1(n), se reflectă înapoi sub forma unui ecou Apar două tipuri de ecou: ⇒ ecou apropiat (NE=Near end Echo) care este reflexia semnalului transmis de la sistemul diferenţial 2F/4F din capătul local al PSTN, ⇒ ecou îndepărtat (FE=Far end Echo) care reflexia semnalului transmis de la sistemul diferenţial 2F/4F din capătul distant al PSTN Pe linia telefonică asociată unui modem coexistă următoarele semnale:  semnalul transmis de modemul respectiv, care are un nivel de putere relativ mare,  semnalul recepţionat de la modemul corespondent, care are un nivel de putere relativ mic,  ecoul apropiat, care are un nivel de putere relativ mare,  ecoul distant, care are un nivel de putere relativ mic Se pune problema de a separa aceste semnale astfel încît modemul local să preia doar semnalul recepţionat de la modemul corespondent Semnalul recepţionat de modemul local este egal cu : ′ ISI)n(z)n(s)n(s)n(s)n(r+ ++= (5 35) FENE21+ 181 -Teza de doctorat- unde: ¾ s2(n) este semnalul transmis de modemul distant, ¾ sNE(n) este ecoul apropiat, ¾ sFE(n) este ecoul îndepă rtat, ¾ z(n) este zgomotul aditiv, alb, gaussian introdus de canalul de comunicaţie Suprimarea ecoului se face prin extragerea unui estimat al ecoului total din semnalul recepţionat de modemul local Astfel semnalul recepţionat de modemul local după suprimarea ecoului devine : ′ ISI)n(z)]n(sˆ)n(s[)]n(sˆ)n(s[)n(s)n(sˆ)n(r)n(r+ −+−+=−= FEFENENE2E11+ (5 36) unde estimatul ecoului total este: (5 37) )n(sˆ)n(sˆ)n(sˆ+ FENEE= Estimarea ecoului apropiat se face prin trecerea semnalului transmis de modemul local s1(n) printr-un filtru adaptiv liniar, transversal, FIR, de ordinul N1 cu funcţia de transfer HNE(z) Estimarea ecoului îndepărtat se face prin trecerea semnalului transmis de modemul local s1(n) printr-o linie de întîrzîiere cu N1 elemente ş i printr-un filtru adaptiv liniar, transversal, FIR, de ordinul N2 cu funcţ ia de transfer HFE(z) Algoritmul adaptiv utilizat pentru ajustarea coeficienţilor de filtrare a celor două filtre adaptive este LMS (Least Mean Squares) (Fig 5 18) Fig 5 18 Supresorul de ecou 182 -Teza de doctorat- Filtrul adaptiv, cu funcţia de transfer HNE(z) sau HFE(z), trebuie să fie suficient de lung astfel încît să fie capabil să aproximeze răspunsul la impuls şi întîrzîierea asociată funcţiei de transfer pe care încearcă să o adapteze O linie telefonică are un răspuns la impuls şi o întîrzîiere asociată funcţiei de transfer de apoximativ 13 milisecunde O conexiune telefonică prin PSTN are un răspuns la impuls şi o întîrzîiere Dc asociată funcţiei de transfer care depinde de lungimea acestei conexiuni: ¾ zeci de milisecunde pentru conexiuni locale, ¾ sute de milisecunde pentru conexiuni interurbane, ¾ cîteva secunde pentru conexiuni internaţionale Dacă frecvenţa de eşantionare a circuitului de codare/decodare (CODEC) este fe [eşantioane/s] atunci numărul de elemente de întărzîiere ale filtrelor adaptive cu funcţiile de transfer HNE(z), respectiv HFE(z) sunt : N1=0,013fe N2=Dc fe (5 38) unde Dc este întîrzîierea asociată funcţiei de transfer a conexiunii telefonice [s] De exemplu, pentru fe=9600 eş antioane/s şi pentru Dc=50 ms rezultă N1=125 elemente de întîrzîiere şi N2=480 elemente de întîrzîiere Antrenarea supresorului de ecou are loc în faza 3 a procedurii de negociere dintre modemul local şi modemul distant Secvenţa de antrenare se obţine prin aplicarea de biţi „1” la intrarea cifratorului 5 2 Implementare hard Schema bloc hard a uni modem analogic de bandă vocală este următoarea: Fig 5 19 Schema bloc hard a unui modem analogic de bandă vocală 183 -Teza de doctorat- 5 2 1 Procesorul digital de semnal (DSP) Procesorul digital de semnal (DSP=Digital Signal Processor) este un chip de tip microprocesor, dedicat pentru realizarea funcţiilor de procesare digitală a semnalelor Un microprocesor dintr-un PC sau un simplu microcontroller nu pot realiza funcţiile complexe pe care le realizează un procesor digital de semnal, pentru că procesarea digitală a semnalelor implică utilizarea extensivă a operaţiilor de multiplicare şi însumare În general, microprocesoarele convertesc operaţiile de multiplicare şi de divizare într-o serie de deplasări, însumări sau scăderi Din acest motiv un microprocesor nu poate realiza funcţiile complexe de procesare digitală a semnalelor la fel de rapid ca şi un procesor digital de semnal Funcţiile DSP-ului: ⇒ recepţia şi transmisia datelor, ⇒ formatarea şi executarea comenzilor în format AT (AttenTion command), ⇒ monitorizarea progresului apelului care constă în detecţia tonalităţilor (ton de disc, ton de ocupat, ton de fax,…), detecţia curentului de sonerie, detecţia răspunsului abonatului chemat, ⇒ monitorizarea procedurii de negociere şi luarea deciziilor cu privire la deplasarea în starea anterioară/următoare a unei maşini de stare, ⇒ tratarea protocoalelor de nivel legătură de date (V 42 pentru corecţia de erori şi V 42 bis pentru controlul erorilor) 5 2 2 Memoria SRAM Memoria SRAM (Static Random Access Memory) memorează datele statice, adică tabelele de date şi codul program Viteza memoriei SRAM depinde de timpul de acces al DSP-ului De exemplu, pentru un DSP Texas Instruments TMS320C51 care operează la o frecvenţă de ceas de 57 MHz (deci durata unui ciclu de citire/scriere este de 35 ns) este necesară o memorie SRAM cu timpul de acces de 20 ns Capacitatea memoriei SRAM depinde de schema bloc funcţională a modemului şi de algoritmii utilizaţi De exemplu, pentru modem anaolgic de bandă vocală construit conform recomandării ITU-T V 34, memoria SRAM necesară este de 64 Ko 5 2 3 Memoria ROM Memoria ROM (Read Only Memory) memorează codul program al DSP-ului Viteza memoriei ROM depinde de timpul de acces al DSP-ului De exemplu, pentru un DSP Texas Instruments TMS320C51 care operează la o frecvenţă de ceas de 57 MHz (deci durata unui ciclu de citire/scriere este de 35 ns) este necesară o memorie ROM cu timpul de acces de 20 ns Problema este că o astfel de memorie ROM nu există sau dacă ar exista ar fi foarte scumpă (costul ei ar fi prohibitiv) Pentru a rezolva 184 -Teza de doctorat- cerinţa de viteză a memoriei ROM DSP-ul conţine un generator de stări de aşteptare, programat la pornire Fiecare stare de aşteptare prelungeşte un ciclu de citire/scriere cu un ciclu maşină al DSP-ului Numărul stărilor de aşteptare se alege printr-un compromis între: ¾ costul memoriei ROM, ţinînd cont că un număr mic de stări de aşteptare înseamnă o memorie ROM rapidă şi scumpă, ¾ performanţele modemului, ţinînd cont că un număr mare de stări de aşteptare înseamnă o întîrzîiere totală mare introdusă de modem Capacitatea memoriei ROM depinde de schema bloc funcţională a modemului şi de algoritmii utilizaţi De exemplu, pentru modem anaolgic de bandă vocală construit conform recomandării ITU-T V 34, memoria ROM necesară este de 128 Ko 5 2 4 Interfaţa cu calculatorul gazdă Interfaţa cu calculatorul gazdă (Host Interface) este o interfaţă serială de tip RS-232 sau USB (Universal Serial Bus) pentru modemurile externe sau o interfaţă de magistrală de tip PCI (Peripheral Component Interconnect) sau ISA (Industry Standard Architecture) pentru modemurile interne În cazul modemurilor interne interfaţa cu calculatorul gazdă se implementează cu ajutorul unui circuit integrat I/O ASIC (Input/ Output Application Specific Integrated Circuit) care conţine: ¾ logica de interfaţă cu calculatorul gazdă, ¾ logica Plug&Play (PnP) care asigură identificarea modemului şi alocarea de resurse pentru modem de către sistemul de operare al calculatorului gazdă, ¾ logica UART 16550 (Universal Asynchronous Receiver/Transmitter), ¾ logica DMA (Direct Memory Access) care permite înregistrarea de voce în modemurile care suportă aplicaţii multimedia 5 2 5 Interfaţa analogică Schema bloc a interfeţei analogice (AFE=Analog Front End) este prezentată în Fig 5 20 Întrucît în canalul de comunicaţie se utilizează semnale analogice şi posibilităţile de prelucrare a semnalelor analogice sunt limitate, semnalele analogice sunt convertite în semnale digitale înainte de a fi prelucrate de DSP cu un circuit de codare/decodare (CODEC) care operează la o frecvenţă de eşantionare fe [eş antioane/s] sau [Hz] şi conform legii µ În general, circuitul de codare/decodare (CODEC) operează la o antioane/s frecvenţă de eşantionare fe=9600 eş 185 -Teza de doctorat- Fig 5 20 Schema bloc a interfeţei analogice Blocul de control automat al cîştigului extern (AGC extern) are rolul de a menţine semnalul analogic recepţionat la un nivel mediu de energie aproximativ constant Acest lucru este necesar pentru funcţionarea corectă a circuitului de codare/decodare (CODEC) şi a blocului de refacere a semnalului purtător În cazul în care nivelul de putere al semnalului recepţionat scade sub –28 dBm, DSP setează o linie de stare care comandă comutarea unei rezistenţe diferite în circuitul de reacţie al unui amplificator operaţional (AO) ceea ce introduce un cîştig de 12 dB Dacă nivelul de putere al semnalului recepţionat creşte peste –28 dBm, DSP resetează linia de stare Pentru a preveni activarea sau dezactivarea repetată a cîştigului de 12 dB cînd nivelul de putere al semnalului recepţionat variază repetat în jurul valorii de prag de –28 dBm, se introduce un histeresis de 4 dBm între activare şi dezactivare Astfel: ¾ dacă nivelul de putere a semnalului recepţionat scade sub –24 dBm se activează cîştigul, ¾ dacă nivelul de putere al semnalului recepţionat creşte peste –28 dBm se dezactivează cîştigul 5 2 6 Interfaţa cu linia telefonică Interfaţa cu linia telefonică (DAA=Data Arrangement Access) conţine: ¾ circuitele de protecţie la supratensiuni (varistoare), ¾ transformatorul de linie (TL), ¾ terminaţia de 600 Ω, 186 -Teza de doctorat- ¾ sistemul diferenţial 2F/4F, ¾ circuitul de limitare al curentului, ¾ circuitul de detecţie al curentului de sonerie, ¾ circuitul de detecţie a identităţii (numărului) abonatului chemător (Caller ID) Sistemul diferenţial 2F/4F funcţionează ca un supresor de ecou analogic Componentele R şi C ale sistemului diferenţial 2F/4F se aleg astfel încît să modeleze impedanţa liniei telefonice (Fig 5 21) Fig 5 21 Interfaţa cu linia telefonică Presupunem că avem o linie telefonică ideală cu impedanţa egală cu 600 Ω Fie Vx semnalul de la ieş irea amplificatorului operaţional AO1 Rezultă că semnalul de la V x bornele transformatorului de linie TL este 2 Cîştigul intrării invertoare (-) a amplificatorului operaţional AO2 este : − x2FVRR− (5 39) === xinVVV−− 1I2RR Cîştigul intrării neinvertoare (+) a amplificatorului operaţional AO2 este : 42FRRR =+=+= (5 40) xxinVV))(1()V1(V++ 431IRRRR+ 187 -Teza de doctorat- Semnalul VRx se anulează complet în traseul de recepţie analogică : (5 41) = 0VVV= -Rx+ + În cazul unei linii telefonice ideale cu impedanţa egală cu 600 Ω, sistemul diferenţial 2F/4F suprimă complet ecoul rezidual În practică impedanţa liniei telefonice variază atît cu lungimea liniei cît şi cu frecvenţa Din acest motiv în majoritatea cazurilor apare o neadaptare de impedanţă care cauzează un ecou care poate fi suprimat cu un supresor de ecou digital 5 3 Implementare soft Schema bloc soft a unui modem analogic de bandă vocală este prezentată în Fig 5 22 La nivelul microprocesorului din calculatorul gazdă se găsesc următoarele blocuri soft : ⇒ soft API (Application Programming Interface) care sunt o serie de funcţii ale sistemului de operare a calculatorului gazdă pe care modemul le utilizează pentru realizarea anumitor sarcini, ⇒ soft comezi AT (AttenTion command) care asigură formatarea şi executarea de către modem a comenzilor în format AT, ⇒ soft V 42 şi V 42 bis care asigură tratarea de către calculatorul gazdă a protocoalelor de nivel legătură de date (V 42 pentru corecţia de erori şi V 42 bis pentru controlul erorilor), ⇒ soft transmisie sincronă/asincronă care asigură transmisia datelor de către calculatorul gazdă şi modem prin tehnica sincronă sau asincronă, ⇒ soft protocol HDLC (High-level Data Link Control) care asigură tratarea de către calculatorul gazdă a protocolului de nivel legătură de date de tip HDLC La nivelul DSP-ului din modem se găsesc următoarele blocuri soft : ⇒ soft transmisie/recepţie de date (Data Pump) care asigură recepţia şi transmisia propriu-zisă a datelor Softul de transmisie/recepţie de date conţine toate blocurile funcţionale digitale ale modemului analogic de bandă vocală (cifrator/decifrator, codor/decodor TCM 2N-D, codor neliniar, egalizor adaptiv, blocul AGC, supresorul de ecou, modulator/demodulator QAM), ⇒ driver pentru CODEC care permite softului de transmisie/recepţie de date al DSP-ului să controleze şi să utilizeze circuitul de codare/decodare (CODEC) 188 -Teza de doctorat- Fig 5 22 Schema bloc soft a unui modem analogic de bandă vocală La nivelul interfeţei analogice (AFE) din modem se găsesc următoarele blocuri soft : ⇒ softul circuitului de codare/decodare (CODEC) care asigură funcţionarea corectă a acestuia, ⇒ driver pentru interfaţa cu linia telefonică care permite softului circuitului de codare/decodare al interfeţei analogice să controleze şi să utilizeze interfaţa cu linia telefonică La nivelul interfeţei cu linia telefonică din modem nu se găsesc blocuri soft 189 -Teza de doctorat- CAPITOLUL VI Concluzii şi contribuţii personale la modulaţia codată trellis multidimensională 6 1 Concluzii personale cu privire la modulaţia codată trellis multidimensională Concluzia 6 1-Construcţia propusă permite generalizarea construcţiei Wei a codurilor TCM 2N-D şi pentru cazul în care N nu este o putere întreagă a lui 2 şi permite obţinerea unor parametri CER, PAR şi γ mai buni decît construcţia 1 Sterian care devine astfel un caz particular al acestei construcţii (cazul α =) 2 Concluzia 6 2-Costul plătit pentru această îmbunătăţire a performanţelor codurilor TCM 2N-D în cazul în care N nu este o putere întreagă a lui 2, constă în creşterea complexităţii codorului bloc mai ales pentru cazurile în care N ia valori mari Acest cost nu este foarte mare întrucît în practică se construiesc coduri TCM 2N-D cu N≤8 şi în prezent costul memoriilor şi al DSP-urilor este relativ scăzut În Anexa 1 am calculat folosind funcţiile din Microsoft Excel 2002 şi am obţinut un tabel cu parametrii codurilor TCM 2N-D pentru dimensiunile codului N=2, 3, 4, 5, 6, 7, 8 şi pentru două valori Q=7 şi Q=10 biţi/interval de smenalizare 2-D ale ratei de transmisie a codorului TCM 2N-D Concluziile obţinute sunt prezentate în continuare Concluzia 6 3-În cazul în care N nu este o putere întreagă a lui 2, pentru un N fixat, cu cît Q este mai mare, cu atît parametrii CER şi PAR obţinuţi sunt mai mici şi cu atît parametrul γ este mai mare în comparaţie cu parametrii corespunzători ai construcţiei Sterian, în condiţiile menţinerii coeficientului de eroare Nfree Concluzia 6 4-În cazul în care N nu este o putere întreagă a lui 2, pentru un Q fixat, cu cît N este mai mare, cu atît îmbunătăţirea parametrilor CER, PAR şi γ obţinuţi relativ la construcţia Sterian este mai mică, dar numărul total de puncte de semnal 2-D a constelaţiilor de semnale 2-D constituente a constelaţiei de semnale 2N-D, este din ce în ce mai mic faţă de construcţia Sterian Concluzia 6 5-În care N este un număr întreg pozitiv, pentru un N şi un Q fixat, cu 2  d   = este mai mare, cu atît cîştigul de codare asimptotic cît cîştigul de distanţă ε   0d  γ relativ la o transmisie necodată este mai mare 190 -Teza de doctorat- Concluzia 6 6-La alegerea valorilor pentru 2N dimensiunea codului TCM, rata de transmisie Q a codorului TCM 2N-D, cîştigul de distanţă ε, numărul m de biţi de intrare codaţi convoluţional şi numărul V de stări ale codorului convoluţional este necesar să facem un compromis între complexitatea codorului TCM 2N-D şi valorile obţinute pentru parametrii CER, PAR şi γ În Anexa 2 am realizat 27 aplicaţii în MathWorks MATLAB-Simulink 6 0 şi anume: Codoare şi decodoare TCM: 01) codor TCM 4-D (construcţia Wei), 02) decodor TCM 4-D cu detecţie SSD (construcţia Wei), 03) decodor TCM 4-D cu detecţie MLD (construcţia Wei) 04) codor TCM 6-D (construcţia Sterian), 05) decodor TCM 6-D şi cu detecţie SSD (construcţia Sterian), 06) decodor TCM 6-D şi cu detecţie MLD (construcţia Sterian), 07) codor TCM 6-D (construcţia propusă), 08) decodor TCM 6-D şi cu detecţie SSD (construcţia propusă), 09) decodor TCM 6-D şi cu detecţie MLD (construcţia propusă), Sisteme de transmisii de date cu TCM şi fără modulaţie QAM: 10) sistem de transmisii de date cu TCM 4-D şi cu detecţie SSD (construcţia Wei), 11) sistem de transmisii de date cu TCM 4-D şi cu detecţie MLD (construcţia Wei) 12) sistem de transmisii de date cu TCM 6-D şi cu detecţie SSD (construcţia Sterian), 13) sistem de transmisii de date cu TCM 6-D şi cu detecţie MLD (construcţia Sterian), 14) sistem de transmisii de date cu TCM 6-D şi cu detecţie SSD (construcţia propusă), 15) sistem de transmisii de date cu TCM 6-D şi cu detecţie MLD (construcţia propusă), Sisteme de transmisii de date cu TCM şi cu modulaţie QAM în banda de bază: 16) sistem de transmisii de date cu TCM 4-D şi cu detecţie SSD (construcţia Wei), 17) sistem de transmisii de date cu TCM 4-D şi cu detecţie MLD (construcţia Wei) 18) sistem de transmisii de date cu TCM 6-D şi cu detecţie SSD (construcţia Sterian), 19) sistem de transmisii de date cu TCM 6-D şi cu detecţie MLD (construcţia Sterian), 20) sistem de transmisii de date cu TCM 6-D şi cu detecţie SSD (construcţia propusă), 21) sistem de transmisii de date cu TCM 6-D şi cu detecţie MLD (construcţia propusă), 191 -Teza de doctorat- Sisteme de transmisii de date cu TCM şi cu modulaţie QAM în banda de trecere: 22) sistem de transmisii de date cu TCM 4-D şi cu detecţie SSD (construcţia Wei), 23) sistem de transmisii de date cu TCM 4-D şi cu detecţie MLD (construcţia Wei) 24) sistem de transmisii de date cu TCM 6-D şi cu detecţie SSD (construcţia Sterian), 25) sistem de transmisii de date cu TCM 6-D şi cu detecţie MLD (construcţia Sterian), 26) sistem de transmisii de date cu TCM 6-D şi cu detecţie SSD (construcţia propusă), 27) sistem de transmisii de date cu TCM 6-D şi cu detecţie MLD (construcţia propusă) Pentru fiecare sistem de transmisii de date de mai sus, am rulat simulări în MathWorks MATLAB-Simulink 6 0 pe un PC Compaq cu sistemul de operare Microsoft Windows 2000 Professional şi avînd o memorie RAM de 128 Mo şi un procesor Intel Pentium III cu frecvenţa de ceas de 667 MHz (hardisk Quantum Fireball 20 Go, placă video Nvidia Vanta 8 Mo, placă audio Creative ES1371, CD-ROM Compaq CDR-8402B) Fiecare simulare am realizat-o pentru un număr total de 10000 de biţi transmişi şi pentru 12 valori ale raportului semnal/zgomot (SNR=Signal to Noise Ratio) din canalul de comunicaţie AWGN (în paşi de 1 dB) Am obţinut numărul de biţi recepţionaţi eronaţi şi rata erorilor de bit (BER=Bit Error Rate) şi am reprezentat grafic BER în funcţie de SNR (Anexa 3) Concluziile sunt prezentate în continuare Concluzia 6 7-Atît în cazul în care nu utilizăm modulaţia QAM cît şi în cazul în care utilizăm modulaţia QAM, sistemul de transmisii de date cu TCM 6-D (construcţia propusă) are rata erorilor de bit mai mică decît sistemul de transmisii de date cu TCM 6-D (construcţia Sterian) Astfel am obţinut un cîştig mediu de 0,2 dB pentru sistemul de transmisii de date cu TCM 6-D (construcţia propusă) relativ la sistemul de transmisii de date cu TCM 6-D (construcţia Sterian) Concluzia 6 8-Dacă utilizăm modulaţia QAM în banda de trecere, se obţin valori mai mici pentru rata erorilor de bit decît în cazul în care nu utilizăm modulaţia QAM sau utilizăm modulaţia QAM în banda de bază Raportul semnal/zgomot pentru care am obţinut o transmisie cu rata erorilor egală cu zero este 4 dB dacă utilizăm modulaţia QAM în banda de trecere şi 11 dB dacă nu utilizăm modulaţia QAM sau dacă utilizăm modulaţia QAM în banda de bază Concluzia 6 9-Dacă utilizăm modulaţia QAM, se obţin valori mai mici pentru rata erorilor de bit a sistemelor de transmisii de date cu detecţie MLD decît a celor cu detecţie SSD 192 -Teza de doctorat- 6 2 Contribuţii personale la modulaţia codată trellis multidimensională 1 Generalizarea construcţiei Wei a codurilor TCM 2N-D Am utilizat ideea extinderii unei constelaţii de semnale 2-D care a fost formulată de Dinh şi Hashimoto în articolul „A systematic approach to construction of bandwith- efficient multidimensional trellis codes”, IEEE Trans Comms , Vol 48, No 11, Nov 2000 Am demonstrat formulele pentru energia medie (Emedie) ş i parametrii PAR şi CER ai constelaţiei de semnale 2-D constituentă a constelaţiei de semnale 2N-D (paragraful 4 4 1 2) În paragraful 4 4 am propus o construcţie Wei generalizată pentru codurile TCM 2N-D Întrucît formulele stabilite de Wei şi Sterian pentru cîştigul de codare asimptotic γ al codului TCM 2N-D faţă de o transmisie necodată, nu ţin cont de parametrii PAR şi CER ai constelaţiei de semnale 2-D constituentă a constelaţiei de semnale 2N-D, am dedus o nouă formulă care ia în calcul parametrii respectivi (relaţia 4 134) În relaţiile folosite de Wei şi Sterian (relaţiile 4 2 şi 4 33), energiile medii ale grupurilor se calculează dificil prin însumarea pătratelor distanţelor dintre originea axelor de coordonate a planului complex şi punctele de semnal 2-D situate în grupurile respective şi prin împărţirea la numărul de puncte de semnal 2-D ale grupurilor Spre deosebire de relaţiile folosite de Wei şi Sterian, relaţia propusă (4 134) permite calculul imediat a cîştigului de codare asimptotic γ al codului TCM 2N-D, după calculul parametrului CER a constelaţiei de semnale 2-D constituentă a constelaţiei de semnale 2N-D 2 Generalizarea construcţiei Wei a codorului TCM 2N-D În anul 1987, Wei a prezentat trei codoare TCM multidimensionale pentru cazul în care N este o putere întreagă a lui 2:  un codor TCM 4-D cu 16 stări (paragraful 4 2 1),  un codor TCM 4-D cu 64 stări,  un codor TCM 8-D cu 64 stări (paragraful 4 2 2) În anul 1997, Sterian a găsit o metodă de extindere a construcţiei Wei pentru cazurile în care N nu este o putere întreagă a lui 2 şi a prezentat două codoare TCM multidimensionale:  un codor TCM 6-D cu 64 stări (paragraful 4 3 1),  un codor TCM 12-D cu 256 stări (paragraful 4 3 2) În paragraful 4 5 am propus un codor TCM 2N-D pentru cazul în care N este un număr întreg pozitiv Am propus un codor diferenţial pe 2 biţi (Fig 4 12) care asigură condiţiile de invarianţă rotaţională a constelaţiei de semnale 2N-D la o rotaţie de 0°, 90°, 180°, 270° şi am dedus relaţiile dintre intrările şi ieşirile sale (relaţiile 4 132) 193 -Teza de doctorat- 3 Generalizarea construcţiei Wei a decodorului TCM 2N-D În paragraful 4 6 am propus un decodor TCM 2N-D pentru cazul în care N este un număr întreg pozitiv Pentru decodorul Viterbi, am construit două diagrame de decodare Viterbi (Fig 4 14 şi Fig 4 15) şi o organigramă de alegere a diagramei Viterbi (Fig 4 13) Am propus un decodor diferenţial pe 2 biţi (Fig 4 17) care asigură condiţiile de invarianţă rotaţională a constelaţiei de semnale 2N-D la o rotaţie de 0°, 90°, 180°, 270° şi am dedus relaţiile dintre intrările şi ieşirile sale (relaţiile 4 135) 4 Construcţia unui modem analogic de bandă vocală cu TCM 2N-D În capitolul 5 am propus o aplicaţie a codorului/decodorului TCM 2N-D la construcţia unui modem analogic de bandă vocală 5 Realizarea în MathWorks MATLAB-Simulink 6 0 a 27 de aplicaţii Lista acestor aplicaţii este prezentată în paragraful 6 1 şi aplicaţiile se găsesc pe CD-ul ataşat în Anexa 2 Pîna în prezent am văzut doar realizări în MathWorks MATLAB-Simulink ale unui codor/decodor TCM 2-D În Anexa 3 am prezentat rezultatele simulărilor pe care le-am realizat în MathWorks MATLAB-Simulink 6 0 194 -Teza de doctorat- Bibliografie citată în teza de doctorat 1 M V EYUBOGLU -„Multimedia modems”-IEEE Comms Mag , Vol 34, No 12, pag 26-28, Dec 1996 2 M ALI -„Evolution from Voiceband to Broadband Internet Access”-DSPS R&D Center, Texas Instruments, 2000 3 GEORGIA TECH Studies - http://www gvu gatech edu/user surveys/survey-1998- 10/graphs/technology/q01 htm 4 B GOODMAN -„Internet telephony and modem delay”-IEEE Network, May/June, 1999 5 W CONWAY -„A management briefing on V 34 modems: History, technology and economics”, General DataComm, 1996 6 M NAFORNIŢĂ, C MUNTEAN -„Comunicaţii de date”-Ed „Gheorghe Asachi” Iaşi, 1996 7 E BIGLIERI, D DIVSALAR, P J McLANE şi M K SIMON -„Introduction to Trellis Coding Modulation with applications”-Macmillan Publishing Company New York, 1991 8 Z PAPIR, A SIMMONDS -„Competing for throughput in the local loop”-IEEE Comms Mag , Vol 37, No 5, pag 61-67, May 1999 9 J L MASSEY -„Coding and modulation in digital communications”, 1974 International Zurich Seminar on Digital Communications, Zurich, Switzerland, March, 1974 10 G UNGERBOECK, I CSAJKA -„On improving data-link performance by increasing the channel alphabet and introducing sequence coding”, 1976 International Symposium of Information Theory, Ronneby, Sweden, Jun 1976 11 H IMAI, S HIRAKAWA -„A new multilevel coding method using error- correcting codes”-IEEE Trans Inf Theory, Vol IT-23, No 3, pag 371-377, May 1977 12 G UNGERBOECK -„Channel coding with multilevels/phase signals”-IEEE Trans Inf Theory, Vol IT-28, No 1, pag 56-67, Jan 1982 13 G D FORNEY Jr , R G GALLAGER, G R LANG, F M LONGSTAFF, S U QURESHI -„Efficient modulation for bandlimited channels”-IEEE J Select Areas Comms, Vol SAC-2, No 5, pag 632-647, Sept 1984 14 L -F WEI -„Rotationally invariant convolutional channel coding with expanded signal space Part II: Nonlinear Coding”-IEEE J Select Areas Comms, Vol SAC-2, No 5, pag 672-686, Sept 1984 15 ITU-T „Recommendation V 32 for a family of 2-wire duplex modems operating at data signalling rates of up to 9600 bit/s for use on the General Switched Telephone Network and on leased telephone-type circuits”, May 1984 16 L -F Wei „Trellis-Coded Modulation with multidimensional constellations”-IEEE J Trans Inf Theory, Vol IT-33, No 4, pag 483-501, July 1987 17 ITU-T „Recommendation V 34 for a family of 2-wire duplex modems operating at data signalling rates of up to 33600 bit/s for use on the General Switched Telephone Network and on leased point-to-point 2-wire telephone-type circuits”, May 1994 18 R G C Williams -„A trellis code for V fast”, CCITT V fast Rapporteur Meeting, Bath , UK, Sept 1992 195 -Teza de doctorat- 19 L -F WEI -„A new 4-D, 64-state, rate-4/5 Trellis Code”-Cont D19, ITU-T, Study Group 14, Geneva, Switzerland, Sept 1993 20 ITU-T „Recommendation V 90 for a digital modem and analogue modem pair operating for use on the Public Switched Telephone Network (PSTN) at data signalling rates of up to 56000 bit/s downstream and up to 33600 bit/s upstream”, Sept 1998 21 ITU-T „Recommendation G 711 for Pulse Code Modulation of voice frequencies”, 1988 22 G UNGERBOECK -„TCM with redundant signals sets Part I: Introduction Part II: State of the art”-IEEE Comms Magazine, Vol 25, No 2, pag 5-21, Feb 1987 23 A R CALDERBANK, J E MAZO -„A new description of trellis codes”-IEEE Trans Inf Theory, Vol IT-30, pag 784-791, Nov 1984 24 O AGAZZI, D G MESSERSCHMITT şi D A HODGES -„Nonlinear echo cancellation of data signals”-IEEE Trans Comms , Vol COM-30, Nr 11, pag 2421- 2433, Nov 1982 25 G D FORNEY Jr -„The Viterbi Algorithm”-IEEE Proc , Vol 61, pag 268-278, 1973 26 BENEDETTO, M AJMONE MARSAN, G ALBETENGO, E GIACHIN - „Combined coding and modulation Theory and applications”-IEEE Trans Inf Theory, Vol IT-34, No 2, pag 223-236, Mar 1988 27 M M MULLIGAN, S G WILSON -„An improved algorithm for evaluating trellis phase codes”-IEEE Trans Inf Theory, Vol IT-30, No 6, pag 846-851, Nov 1984 28 E BIGILERI-„High-level modulation and coding for nonlinear satellite channels”- IEEE Trans Comms , Vol COM-32, No 5, pag 616-626, May 1984 29 G D FORNEY Jr -„Coset codes Part I: Introduction and geometrical clasification”-IEEE Trans Inf Theory, Vol IT-34, No 5, pag 1123-1151, Sept 1988 30 G J POTTIE, D P TAYLOR -„Sphere-packing upper bounds on the free distance of trellis codes”-IEEE Trans Inf Theory, Vol IT-34, No 2, pag 435-447, May 1988 31 A R CALDERBANK, J E MAZO, V K WEI -„Asymtotic upper bounds on the minimum distance of trellis codes”-IEEE Trans Comms , Vol COM-33, pag 305-309, Apr 1985 32 E BIGLIERI -„Ungerboeck codes do not shape the signal power spectrum”-IEEE Trans Inf Theory, Vol IT-32, No 4, pag 595-596, Jun 1986 33 G D FORNEY Jr -„Geometrically uniform codes”, IEEE Trans Inf Theory, Vol IT-37, No 9, pag 1241-1260, Sept 1991 34 A R Calderbank, N J A Sloane -„Four-dimensional modulation with eight-state trellis code”, AT&T Tech Journal, Vol 64, pag 1005-1018, May/June 1985 35 S S PIETROBON , D J COSTELLO -„Trellis coding modulation with multidimensional QAM signals sets”, IEEE Trans Inf Theory, Vol IT-39, No 3, pag 325-336, Mar 1993 36 S S PIETROBON, R H DENG, A LAFANACHERE, G UNGERBOECK, D J COSTELLO -„Trellis-coded multidimensional phase modulation”, IEEE Trans Inf Theory, Vol IT-36, No 3, pag 63-89, Jan 1990 37 F Q WANG, D J COSTELLO -„New rotationally invariant four-dimensional trellis codes”, IEEE Trans Inf Theory , Vol IT-42, pag 291-300, Jan 1996 196 -Teza de doctorat- 38 T C DINH, T HASHIMOTO -„A systematic approach to construction of bandwith-efficient multidimensional trellis codes”, IEEE Trans Comms , Vol 48, No 11, pag 1808-1817, Nov 2000 39 P PRANDONI -„Information theory in modem practice”, Laboratoire de Communication Audio Visuelle (LCAV) Tech-Report, École Polytechnique Fédérale de Lausanne, Feb 1997 40 C -E D STERIAN -„Extending Wei’trellis coded modulation technique with 2N- D rectangular constellation for N not a power of 2”, Telecomunicaţ, Anul XXI, Nr 4, ii 1994 41 C -E D STERIAN -„128-state, rate-4/5 rotationally invariant trellis code with 4- D rectangular signal constellation having a coding gain of 5,63 dB”, Telecomunicaţ, ii Anul XXIII, Nr 2, 1996 42 C -E D STERIAN -„Wei-type trellis-coded modulation with 2N-dimensional rectangular constellation for N not a power of 2”, IEEE Trans Inf Theory, Vol IT-43, No 3, pag 750-758, Mar 1997 43 G D FORNEY Jr -„Trellis Shaping”, IEEE Trans Inf Theory, Vol IT-38, No 3, pag 281-300, Mar 1992 44 A K KHANDANI, P KABAL -„Shaping Multidimensional Spaces - Part I: Optimum shaping, Shell Mapping”, IEEE Trans Inf Theory, Vol IT-39, No 11, pag 1799-1808, Nov 1993 45 Black Box Inc -„V 34”, Technology Overview, 1999 - http://www blackbox com 46 G D FORNEY Jr , M V EYUBOGLU -„Combined equalization and coding using precoding”-IEEE Comms Mag , Vol 29, No 12, pag 23-34, Dec 1991 47 G D FORNEY Jr , M V EYUBOGLU -„Trellis precoding: combined coding , precoding and shaping for intersymbol interference channels”, IEEE Trans Inf Theory, Vol IT-38, No 2, pag 301-313, Mar 1992 48 W H GERSTAKER, R F H FISCHER, J B HUBER -„A transmission scheme for twisted pair lines with coding, precoding and blind equalization”, Proceedings of IEEE Global Telecommunications Conference (Globecomm’97), pag 52-56, Phoenix, AZ, Nov 1997 49 F HARRIS –„Carrier synchronization techniques for DSP-based modems”, Communication Systems Design Magazine, July 2000 - http://www csdmag com/main/2000/07/0007feat2 htm 50 T MASSEY, R IYER - „DSP solutions for telephony and data/facsimile modems”, Texas Instruments - Application Book, Jan 1997 51 L LITWIN -„Matched filtering and timing recovery in digital receivers”, Thompson Multimedia Corporate Research, Sept 2001 52 J WENG -„Digital Transmissions Systems”, course 2002-2003, Department of Electrical and Computer Engineering, Concordia University , Montreal, Canada 53 A LOVRICH, R CHIRAYIL - „An all-digital automatic gain control”, Texas Instruments - Application Report, Jan 1997 54 F GAO -„Modem and fax standards and software”, GAO Research and Consulting Article , 2001 - http://www gaoresearch com/resources/articles/mfstds html 55 E KAREN -„ITU-T Recommendation V 34”, The University of Toronto, 2000 - http://www comm toronto edu/~karen/projects/38 ITUV34/index html 197 -Teza de doctorat- 56 F DĂRĂBAN -„Generalization of Wei's construction for the 2N-dimensional TCM codes“, Proceedings of the International Workshop “Trends and Recent Achivements in Information Technology”, pag 197-205, Cluj Napoca, Romania, May 2002 - http://193 226 6 174/IT2002/pdf/L35 pdf 57 F DĂRĂBAN -„Generalization of Wei's construction for the 2N-dimensional TCM encoder and decoder“, Proceedings of the International Workshop “Trends and Recent Achivements in Information Technology”, pag 205-213, Cluj Napoca, Romania, May 2002 - http://193 226 6 174/IT2002/pdf/L36 pdf 58 F DĂRĂBAN -„An analog voice modem with generalized Wei's type 2N- dimensional trellis coded modulation“, Proceedings of the International Workshop “Trends and Recent Achivements in Information Technology”, pag 213-223, Cluj Napoca, Romania, May 2002 - http://193 226 6 174/IT2002/pdf/L37 pdf 59 F DĂRĂBAN -„A Wei’s Type Multidimensional TCM Encoder“, Transactions on Electronics and Communications - Buletinul Stiinţific al Universităţii Politehnica Timişoara, pag 3-12, Tom 46(60), Fascicola 2, 2001 60 F DĂRĂBAN -„A Wei’s Type Multidimensional TCM Decoder“, Transactions on Electronics and Communications - Buletinul Stiinţific al Universităţii Politehnica Timişoara, pag 13-18, Tom 46(60), Fascicola 2, 2001 198 -Teza de doctorat- Bibliografie consultată pentru documentare 61 M E BORDA -„Teoria transmiterii informaţiei Teoria informaţiei şi codării Fundamente şi aplicaţii” -Ed Dacia, Cluj-Napoca, 1999 62 T STARR, J M CIOFFI, P J SILVERMAN -„Understanding Digital Subscriber Line Technology” -Prentice Hall PTR, 1999 63 S RAJPAL, D J RHEE, S LIN -„Multidimensional Trellis Coded Phase Modulation Using a Multilevel Concatenation Approach - Part I: Code Design” -IEEE Trans Comms , Vol COM-45, Nr 01, pag 64-72, Jan 1997 64 S RAJPAL, D J RHEE, S LIN -„Multidimensional Trellis Coded Phase Modulation Using a Multilevel Concatenation Approach - Part II: Codes for the AWGN and Fading Channels” -IEEE Trans Comms , Vol COM-45, Nr 02, pag 177- 186, Feb 1997 65 A MASOOMZADEH-FARD, S PASUPATHY -„Combined Equalization and Differential Detection Using Precoding” -IEEE Trans Comms , Vol COM-45, Nr 03, pag 274-278, Mar 1997 66 S -L SU, J -M WU -„Combination of BCM and TCM with 90º, 180º and 270º Phase Rotational Invariance” -IEEE Trans Comms , Vol COM-45, Nr 07, pag 800- 808, July 1997 67 J -Y WANG, M -C LIN -„On Constructing Trellis Codes with Large Free Distances and Low Decoding Complexities” -IEEE Trans Comms , Vol COM-45, Nr 09, pag 1017-1020, Sept 1997 68 E BACCARELLI, R CUSANI, G DI BLASIO -„Performance Bound and Trellis- Code Design Criterion for Discrete Memoryless Channels and Finite-Delay Symbol- by-Symbol Decoding” -IEEE Trans Comms , Vol COM-45, Nr 10, pag 1192-1199, Oct 1997 69 T JI, P AN, S C KWATRA -„Nonlinear Rotationally Invariant Trellis Codes for Multidimensional Modulation” -IEEE Trans Comms , Vol COM-45, Nr 10, pag 1231-1235, Oct 1997 70 E S ESTEVES, R SAMPAIO-NETO -„A Per-Survivor Phase Acquisition and Tracking Algorithm for Detection of TCM Signals with Phase Jitter and Frequency Error” -IEEE Trans Comms , Vol COM-45, Nr 11, pag 1381-1384, Nov 1997 71 E BACCARELLI, S GALLI -„New Results About Analysis and Design of TCM for ISI Channels and Combined Equalization/Decoding” -IEEE Trans Comms , Vol COM-46, Nr 04, pag 417-420, April 1998 199 -Teza de doctorat- 72 J K CAVERS, J -H KIM, P HO -„Exact Calculation of the Union Bound on Performance of Trellis-Coded Modulation in Fading Channels” -IEEE Trans Comms , Vol COM-46, Nr 05, pag 576-579, May 1998 73 K M MACKENTHUN JR -„Codes Based on a Trellis Cut-Set Transformation - Part I: Rotationally Invariant Codes” -IEEE Trans Comms , Vol COM-46, Nr 05, pag 610-620, May 1998 74 J E SMEE, N C BEAULIEU -„Error-Rate Evaluation of Linear Equalization and Decision Feedback Equalization with Error Propagation” -IEEE Trans Comms , Vol COM-46, Nr 05, pag 656-665, May 1998 75 E AYANOGLU, N R DAGDEVIREN, G D GOLDEN, J E MAZO -„An Equalizer Design Technique for the PCM Modem: A New Modem for the Digital Public Switched Network” -IEEE Trans Comms , Vol COM-46, Nr 06, pag 763-774, June 1998 76 J LABAT, O MACCHI, C LAOT -„Adaptive Decision Feedback Equalization: Can You Skip the Training Period?” -IEEE Trans Comms , Vol COM-46, Nr 07, pag 921-930, July 1998 77 T -T CHEN, J S LEHNERT -„TCM/SSMA Communication Systems with Cascaded Sequences and PAM/QAM Signal Sets” -IEEE Trans Comms , Vol COM- 46, Nr 07, pag 950-956, July 1998 78 D G LAINIOTIS, P PAPAPARASCHEVA -„Partitioned Adaptive Approach to Nonlinear Channel Equalization” -IEEE Trans Comms , Vol COM-46, Nr 10, pag 1325-1336, Oct 1998 79 A E VITYAEV, P H SIEGEL -„On Viterbi Detector Path Metric Differences” - IEEE Trans Comms , Vol COM-46, Nr 12, pag 1549-1554, Dec 1998 80 W LEE, K CHEUN -„Convergence Analysis of the Stop-and-Go Blind Equalization Algorithm” -IEEE Trans Comms , Vol COM-47, Nr 02, pag 177-180, Feb 1999 81 D HUANG, F GUSTAFSSON -„Sufficient Output Conditions for Identifiability in Blind Equalization” -IEEE Trans Comms , Vol COM-47, Nr 02, pag 191-194, Feb 1999 200 -Teza de doctorat- 82 H A RAJAB, M D YUCEL -„Efficient Performance Computations for Trellis- Coded Modulation” -IEEE Trans Comms , Vol COM-47, Nr 07, pag 796-798, June 1999 83 K M MACKENTHUN Jr -„Codes Based on a Trellis Cut Set Transformation - Part II: Codes for Noncoherent Detection” -IEEE Trans Comms , Vol COM-47, Nr 07, pag 998-1007, July 1999 84 K J KIM, R A ILTIS -„An Importance Sampling Technique for a Symbol-by- Symbol TCM Decoding/Equalization Algorithm” -IEEE Trans Comms , Vol COM- 48, Nr 07, pag 1141-1150, July 2000 85 M -C LIN, Y -L UENG, J -Y WANG -„Two Trellis Coding Schemes for Large Free Distances” -IEEE Trans Comms , Vol COM-48, Nr 08, pag 1286-1296, Aug 2000 86 Z DING, Z -Q LUO -„A Fast Linear Programming Algorithm for Blind Equalization” -IEEE Trans Comms , Vol COM-48, Nr 09, pag 1432-1436, Sept 2000 87 W H GERSTACKER, R R MULLER, J B HUBER -„Iterative Equalization with Adaptive Soft Feedback” -IEEE Trans Comms , Vol COM-48, Nr 09, pag 1462- 1467, Sept 2000 88 J X YU, Y LI, H MURATA, S YOSHIDA -„Hybrid-ARQ Scheme Using Different TCM for Retransmission” -IEEE Trans Comms , Vol COM-48, Nr 10, pag 1609-1613, Oct 2000 89 D RAPHAELI, T KAITZ -„A Reduced-Complexity Algorithm for Combined Equalization and Decoding” -IEEE Trans Comms , Vol COM-48, Nr 11, pag 1797- 1807, Nov 2000 90 E BACCARELLI -„On the Performance Limits of TCM in Fast-Fading Multipath Channels with Combined Equalization/Decoding” -IEEE Trans Comms , Vol COM- 48, Nr 12, pag 1957-1964, Dec 2000 91 J MANNERKOSKI, V KOIVUNEN, D P TAYLOR -„Performance Bounds for Multistep Prediction-Based Blind Equalization” -IEEE Trans Comms , Vol COM-49 Nr 01, pag 84-93, Jan 2001 201 -Teza de doctorat- 92 S A ALTEKAR, A E VITYAEV, J K WOLF -„Decision-Feedback Equalization via Separating Hyperplanes” -IEEE Trans Comms , Vol COM-49, Nr 03, pag 480-486, Mar 2001 93 L M GARTH -„A Dynamic Convergence Analysis of Blind Equalization Algorithms” -IEEE Trans Comms , Vol COM-49, Nr 04, pag 624-634, April 2001 94 K R NARAYANAN, G L STUBER -„Performance of Trellis-Coded CPM with Iterative Demodulation and Decoding” -IEEE Trans Comms , Vol COM-49, Nr 04, pag 676-687, April 2001 95 Q WANG, L WEI -„Iterative Viterbi Algorithm for Concatenated Multidimensional TCM” -IEEE Trans Comms , Vol COM-50, Nr 01, pag 12-15, Jan 2002 96 M E PELLENZ, J PORTUGHEIS -„ On the Analysis of the Performance of Coded Modulation Schemes for Unequal Error Protection”- IEEE Trans Comms , Vol COM-50, Nr 08, pag 1205-1208, Aug 2002 97 Y -L UENG, C -J ZEH, M -C LIN -„ On Trellis Codes with a Delay Processor and a Signal Mapper”- IEEE Trans Comms , Vol COM-50, Nr 12, pag 1906-1917, Dec 2002 98 W Y CHEN -„The Development and Standardization of Asymmetrical Digital Subscriber Line” -IEEE Comms Mag , Vol 37, Nr 05, pag 68-72, May 1999 99 J F HAYES -„The Viterbi Algorithm Applied to Digital Data Transmission” - th IEEE Comms Mag , 50 Anniversary Issue, pag 26-32, May 2002 st 100 E BIGLIERI -„Digital Transmission in the 21 Century” -IEEE Comms Mag , th 50 Anniversary Issue, pag 128-137, May 2002 101 N V TRAU -„Performance and Distance Spectrum Calculation of TCM over ISI Channels” -Master Thesis, The University of South Australia, 1998 102 S A AL-SEMARI -„Design and Performance of Trellis Codes for Wireless Channels” -PhD Thesis, The University of Maryland, 1995 103 R LAROIA -„An Optimal Shaping for Multidimensional Constellations-An Alternative Approach to Lattice-Bounded (Voronoi) Constellations” -Tech Report, The University of Maryland, 1992 202 -Teza de doctorat- 104 R LAROIA, S TRETTER, N FARVARDIN -„A Simple and Effective Precoding Scheme for Noise Whitening on ISI Channels” -Tech Report, The University of Maryland, 1992 105 S A AL-SAMARI, T E FUJA -„Performance Analysis of Coherent TCM Systems with Diversity Reception in Slow Rayleigh Fading” -Tech Report, The University of Maryland, 1996 203 [biti/i s ] 7777777QN 2 3 4 5 6 7 8 2 2 4 4 4 4 8N ' 4 4 4 4 4 4 4 0,00000 0,09091 0,00000 0,14286 0,16667 0,20000 0,00000 α 0,50000 0,34375 0,25000 0,21875 0,18750 0,15625 0,12500[biti/i s ] 10101010101010QN 1,37500 1,22222 1,15625 1,12000 1,09722 1,08163 1,07031 2 3 4 5 6 7 8 2 2 4 4 4 4 8N ' 1,37500 1,22233 1,15625 1,12012 1,09733 1,08175 1,07031 4 4 4 4 4 4 4 0,00000 0,01163 0,00000 0,01923 0,02326 0,02703 0,00000 2,18182 2,19867 2,16216 2,17611 2,16434 2,13774 2,10219 α 0,50000 0,33594 0,25000 0,20313 0,16797 0,14453 0,12500[biti/i s ] 1,38303 0,87150 0,63052 0,49218 0,40295 0,34080 0,29511 [dB]7777777QN 1,37500 1,22222 1,15625 1,12000 1,09722 1,08163 1,07031 2 3 4 5 6 7 8 1,38303 0,87189 0,63052 0,49263 0,40338 0,34128 0,29511 [dB] 2 2 4 4 4 4 8N ' 1,37500 1,22222 1,15625 1,12000 1,09722 1,08163 1,07031 8 8 8 8 8 8 8 3,38819 3,42160 3,34888 3,37681 3,35326 3,29954 3,22672 [dB] 0,00000 0,09091 0,00000 0,14286 0,16667 0,20000 0,00000 2,18182 2,18608 2,16216 2,14844 2,12895 2,11630 2,10219α 4,63757 5,14871 5,39008 5,52797 5,61722 5,67932 5,72549 [dB] 0,50000 0,34375 0,25000 0,21875 0,18750 0,15625 0,12500[biti/i s ] 1,38303 0,87150 0,63052 0,49218 0,40295 0,34080 0,2951110101010101010 [dB] totalpuncte 192 172 160 156 152 148 144 1,37500 1,22222 1,15625 1,12000 1,09722 1,08163 1,07031 2 3 4 5 6 7 8 1,38303 0,87151 0,63052 0,49218 0,40295 0,34081 0,29511 2 2 4 4 4 4 8 1,37500 1,15625 1,07031 [dB] W E - - - - I 1,37500 1,22233 1,15625 1,12012 1,09733 1,08175 1,07031 8 8 8 8 8 8 8 2,18182 2,16216 2,10219 [dB] 3,38819 3,39665 3,34888 3,32122 3,28166 3,25577 3,22672 W E 0,00000 0,01163 0,00000 0,01923 0,02326 0,02703 0,00000 I 2,18182 2,19867 2,16216 2,17611 2,16434 2,13774 2,10219 1,38303 0,63052 0,29511 [dB] 4,63757 5,14909 5,39008 5,52841 5,61765 5,67979 5,72549 W [dB] E 0,50000 0,33594 0,25000 0,20313 0,16797 0,14453 0,12500 I 1,38303 0,87150 0,63052 0,49218 0,40295 0,34080 0,29511 [dB] puncte 3,38819 3,34888 3,22672total 1536 1368 1280 1232 1196 1172 1152 W [dB] 1,37500 1,22222 1,15625 1,12000 1,09722 1,08163 1,07031 E I 1,38303 0,87189 0,63052 0,49263 0,40338 0,34128 0,29511 1,37500 1,15625 1,07031 [dB] W 4,63757 5,39008 5,72549 - - - - W [dB]EI 1,37500 1,22222 1,15625 1,12000 1,09722 1,08163 1,07031 E I 3,38819 3,42160 3,34888 3,37681 3,35326 3,29954 3,22672 2,18182 2,16216 2,10219 [dB] W total puncte 192 160 144E 2,18182 2,18608 2,16216 2,14844 2,12895 2,11630 2,10219 I 7,64787 8,15901 8,40038 8,53827 8,62752 8,68962 8,73579 1,38303 0,63052 0,29511 [dB] STERIAN W [dB] 1,229171,121881,098961,08259 E 1,38303 0,87150 0,63052 0,49218 0,40295 0,34080 0,29511 - - -I[dB] puncte 3,38819 3,34888 3,22672total 192 172 160 156 152 148 144 W [dB] STERIANE 1,38303 0,87151 0,63052 0,49218 0,40295 0,34081 0,29511 2,44068 2,22841 2,27488 2,30928I 1,37500 1,15625 1,07031 [dB] W 4,63757 5,39008 5,72549E - - - - W [dB]I E I 3,38819 3,39665 3,34888 3,32122 3,28166 3,25577 3,22672 STERIAN 2,18182 2,16216 2,10219 [dB] 0,89611 0,49944 0,40981 0,34464 W [dB] total puncte 1536 1280 1152EI 7,64787 8,15939 8,40038 8,53871 8,62795 8,69009 8,73579 1,38303 0,63052 0,29511 [dB] STERIANSTERIAN W [dB] 3,87510 3,47996 3,56959 3,634761,229171,121881,098961,08259E [dB]I - - - puncte 3,38819 3,34888 3,22672total 1536 1368 1280 1232 1196 1172 1152 W [dB] E STERIAN STERIANI 5,12449 5,52116 5,61079 5,67596 2,44068 2,22841 2,27488 2,30928 1,37500 1,15625 1,07031 [dB] W E - - - - Anexa 1-Tabel cu parametrii calcula W [dB] 7,64787 8,40038 8,73579I E Idd 0 puncte STERIAN W 2,18182 2,16216 2,10219 Nr total 192 160 160 160 0,89611 0,49944 0,40981 0,344642 QNβCER0γNr γNr γ  - [dB] puncteE  total 192 160 144I  CER 1,38303 0,63052 0,29511 N CER CERSTERIAN W [dB] ' STERIAN 3,87510 3,47996 3,56959 3,634761,229171,121881,098961,08259E [dB]I - - - P 3,38819 3,34888 3,22672 PA W [dB] ARE R P STERIANI A STERIAN 2,44068 2,22841 2,27488 2,30928 R 5,12449 5,52116 5,61079 5,67596 [dB] α 7,64787 8,40038 8,73579 CER0 CER W [dB] E Id d CER STERIAN0 puncte 0,89611 0,49944 0,40981 0,34464 βγγ Nr total 1536 1280 1280 1280 [dB]2 CER0Nr Nr γ  P - total puncte 1536 1280 1152 CERA R  CER P CER CER STERIANSTERIAN1,229171,121881,098961,08259 A 3,87510 3,47996 3,56959 3,63476 PR [dB] A - - - R P PA AR ţ i ai codurilor TCM 2N-D pentru N=2R PA STERIAN STERIAN 2,44068 2,22841 2,27488 2,30928 R [dB] 8,13479 8,53146 8,62109 8,68626 CER0 CER 0dd CER puncte STERIAN  βγ Nr γNr total 192 160 160 160 [dB] 0,89611 0,49944 0,40981 0,344642 CER0Nr γ P - CERA  R CER P CER CER STERIAN A 3,87510 3,47996 3,56959 3,63476 PR [dB] A R P PA AR R P STERIAN A 8,13479 8,53146 8,62109 8,68626 R [dB] CER0 CER 0dd CER puncte  βγ Nr γNr total 1536 1280 1280 12802 CER0Nr γ P - CERA  R CER P CER CER A PR A R ≤ P 8 PA AR R P A R CER0 CER CER P CERA204 R P A PR A R Anexa 2-Modelele MATLAB-Sim TCM 6-D (construc CD-ul at a ş uat link pentru sistemele d ţ ia Sterian ) ş i TCM 6-D (construc e t r ansmisii de date cu TCM 4-D (construc ţ ia propus 20 ă5 ) ţ ia W e i), Anexa 3a 206 TC M 6 - D ST TC EM R 6 I A- D P N R OP US Ă ) ER B it ( e b r d rilo o er ata R 0,0,0,0,0,1, 000000001010100000 011000000000 000000000000 000000 01 23 R a po r t ul s eDe m nt e 4alc /ţ zie gom S S otD s ( Ea debu / 5N 0t )e ac lţ ie c Ma n aL lD ul fui, 67ă der ă m c oodu m u nl iaţ c aie ţ Qi e (AA M W G N ) [d 89B ] 1 0 1 1 Anexa 3b 207 TC M 6 - D ST TC EM R I 6- A ND P R O P US Ă ) ER B it ( b e r d rilo o ata er R 0,0,0,0,0,1, 000010 000001010100000000 100000000000000000000000 0 1 2 3 R a p o r t u l s e m 45n a lDe / z got e mc ţ oie St ( ES bD c/ N 0u m) a lo cd au nala lţ uie Ql u deA67i M c Bao m us s ni B ca an ţ ied ( A W G N ) [ d B 8] 9 1 0 1 1 Anexa 3c 208 TC M 6 - D T S TC EM R I 6 A- ND P R O P US Ă ) ER (B it e b r d rilo o er ata R 000001 ,,,,,, 000000011000 000110000000 100000000000 000000000000 012 345 R a po r t ul s e m n aDe l / z gote c mţ oie Mt ( EL b/D cN 0u m) a lo cd au nl aa lţ uie Ql u deA6i M c Bao m us n s Ba i c an ţ id 789e (A W G N) [d B ] 1 0 1 1 Anexa 3d 209 TC M 6 - D T S TC EM R I 6 A- ND P R O P US Ă ) ER (B it e b r d rilo o er ata R 000001 ,,,,,, 000000011000 000110000000 100000000000 000000000000 -6 -5 -4 -3 R a po r t ul sDe e mt e -2nc aţ lie / Sz go S mD co tu ( mE b/o -1Nd 0ula ) ţ ai le QA c a n aM l u Pl u deas 0i s c B oan md u n(fp i c=1 a ţ 8 1ie0 (A0 H Wz ) G N) [d B 2] 3 4 Anexa 3e 210 TC M 6 - D ST TC EM R 6 I A- D PR N O P US Ă ) ER B it ( b e r d rilo o ata er R 0,000010000,000100000,001000000,010000000,100000001,00000000 -6 -5 -4 -3 Ra p o rt ul De s et me -2ncţ aie Ml / z g o L mD co tu (E m b/odul -1N 0 ) aa ţ i Qle c a nA aM l ul Pui da 0 s es cB oa (mnd uni f cp a= ţ Hi1800 (A 1e W Gz N) ) [dB] 2 3 4